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SUMMARY 


A simplified mathematical model is derived that is useful for studying the 
effects of vibration-dissociation coupling in fluid flows. The derivation is based 
on an energy-moment procedure for simplifying the master equations . The 
formalism yields a final set of rate equations that are similar in many respects 
to those of Marrone and Treanor. Instead of a single temperature for the vibra- 
tional energy mode, however, two vibrational temperatures are defined, one for 
the distribution in the lower vibrational states and one for the upper states. In 
this way the perturbing effect of dissociation and recombination on the vibrational 
energy states is approximately accounted for. Also included are the effects of 
molecular anharmonicity, both as regards the energy-level spacing and the 
collisional transition rates. Such effects require determination of the "rate of 
quantum transfer" occurring as a result of molecular collisions, which is approx- 
imated by means of linear relations. The rate equations that result contain the 
necessary modification of the Landau-Teller theory to account for anharmonic 
effects. The effects of dissociation and recombination are treated by a procedure 
similar to that used by Marrone and Treanor, except that only vibrationally ex- 
cited molecules are allowed to dissociate. The thermodynamic quantities asso- 
ciated with the vibrational energy mode are derived, and it is shown that the 
procedure for considering the effects of anharmonicity results in relations that 
agree favorably with the more accurate quantities obtained from spectroscopically 
determined values of the vibrational energies. 

The parameters of the model equations are evaluated by comparison of the 
pre-exponential temperature dependence of the effective dissociation rate with 
experimental results. For this comparison, analytical solutions are obtained 
that apply in the quasi-steady zone behind a normal shock wave. The effect of 
altering certain parameters is also discussed. 

Complete numerical solutions of the model equations are then given for the 
flow behind a steady-shock wave and for nozzle flow. The shock structure is 
discussed in detail including the effects of the various terms that comprise the 
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rate equations. The transient, quasi-steady, and final relaxation zones, that 
have been described by previous researchers, are displayed; and the relaxation 
effects peculiar to the different zones are pointed out. Values of the induction 
time, during which the effects of dissociation are not observed, are found and 
compared with experimentally obtained delay times. The values obtained from 
the calculations are an order of magnitude too large, and the readjustments of 
the assigned parameters that would improve the agreement are indicated. 

The nozzle solution is discussed in detail, and the effects of the various 
terms in the model rate equations are again described. It is shown that the 
characteristic vibrational relaxation time that results in the case of near equi- 
librium flow is about l/4 of that given by Landau-Teller theory and is in quali- 
tative agreement with recent experimental data. However, such reductions of 
the vibrational relaxation time, resulting from consideration of molecular 
anharmonicity effects, are realized only if freezing occurs near the throat when 
the translational and vibrational temperatures are large. 

The equations introduced in this paper contain the minimum structure 
allowing agreement with presently available experimental data. The formalism 
is sufficiently general that additional structure can be included as more refined 
experimental measurements become available. 

The numerical procedure used to integrate the system of rate and flow 


equations is also described. 
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INTRODUCTION 


The purpose of this work is to develop a relatively simple mathematical 
model that is useful for studying the effects of vibration-dissociation coupling in 
fluid flows. The work may be viewed as a logical continuation of the sequence of 
papers by Hammerling, Teare, and Kivel (36), Treanor and Marrone (90, 91, 93), 
and Marrone and Treanor (58). In each of these cases an effort was made to de- 
vise a model containing the smallest possible number of parameters that would 
yield results in agreement with existing experimental data. Such a model would 
be of value for extrapolation of rate processes beyond the temperatures and 
densities accessible to laboratory evaluation. Of equal importance, the resulting 

formalism would lend itself readily to fluid-flow calculations. The model pre- 
* 

sented here is analogous in many respects to those introduced in the aforemen- 
tioned papers, although more complex. The theory is applicable to polyatomic 
species (see, e.g., (41a)), but only diatomic species are considered here. The 
rotational mode is taken to be in equilibrium with the translational (kinetic- 
energy) mode and only the lowest molecular electronic states are considered. 
However, instead of the single vibrational temperature considered in the previous 
papers (36, 58, 90, 91, 93), two independent vibrational temperatures are used to 
allow for the perturbing effects of dissociation and recombination in the upper 
levels. Also, the reduced energy-level spacing of the upper levels is approxi- 
mately accounted for as well as the anharmonic oscillator rates associated with 
these levels. 

Historically, the subject of energy exchange between external degrees of 
freedom (translation) and "internal" degrees of freedom (rotation, vibration, 
dissociation, electronic excitation, ionization) has received considerable atten- 
tion. In addition to the intrinsic interest the subject commands, it is significant 
in the general field of chemical kinetics and has an important bearing on many 
gas-dynamical and combustion problems. Recently an increasing interest in gas 
lasers has added impetus to the need for understanding the mechanisms of 



molecular energy exhange (33, 42, 60). Theoretical treatment of the problem is 
impeded by the complexity arising from the large number of interacting parti- 
cles. During a collision a molecule in a given energy state (specified transla- 
tional, rotational, vibrational, and electronic coordinates) many undergo a 
transition to one of many possible final states. A satisfactory theoretical treat- 
ment must necessarily consider the likelihood and importance of all of these 
transitions (4, 5, 34, 56). Considerable research, both theoretical and experi- 
mental, has been carried out on this subject, involving many specialized studies 
on rate constants and collision or relaxation mechanisms. However, important 
as the prior research is, a comprehensive survey of all the significant papers 
that contribute to an understanding of the complete problem and that would be 
helpful in placing the present work in perspective will not be attempted. Instead 
the reader is referred to the books by Vincenti and Kruger (97), Clarke and 
McChesney (22) (an excellent source for a review of prior research), Stupochenko 
et al. (86), Bradley (10), and ZePdovich and Raizer (100). Herzberg and Litovitz 
(41), Cottrell and McCoubrey (23), and Nikitin (64) are also excellent sources for 
an understanding of rates and relaxation mechanisms. Although a fully compre- 
hensive survey will not be given, it is worthwhile to discuss briefly a few papers 
that are closely related and contribute to the underlying ideas of the specific 
approach adopted here. 

Previously developed schemes for finding first-energy-moment "summed" 
relations of the population equations and thus obtaining the simple relaxation 
equation 

dq°° (Ty) = q°°(T) - q°°(T V ) 

dt T (T,p) ( } 

are basic to the model to be introduced here. In this relation q°° is the average 
vibrational energy (per mole), a known function of the temperature T or Ty, and 
r(T,p) is a characteristic relaxation time that depends on temperature and pres- 
sure. The equation describes vibrational relaxation in the absence of dissocia- 
tion and is often referred to as the Landau-Teller or the Bethe-Teller relax- 
ation model. Landau and Teller (49) were the first to observe that, providing 



the vibrational transition rates vary linearly with vibrational quantum number, 
as do those for a harmonic oscillator (46, 49), a rate equation can be derived 
from the first-energy-moment relations that contains only the single character- 
istic relaxation time r. Later Bethe and Teller (8) identified the summation 
quantities that appear in the equation as the thermodynamic quantity q°°*. The 
basic equations used were the population or "master" equations which express 
the rate of change of the population of each energy state. Imbedded in the 
Landau-Teller or Bethe-Teller approach are the following assumptions: (1) the 
vibrational transitions occur step-wise from one vibrational energy level to the 
next adjacent level ("ladder climbing"); (2) vibrational transition rates vary 
linearly with vibrational quantum number; (3) the energy levels are equally 
spaced; and (4) the molecular population over the various vibrational energy 
states is Boltzmann and thus describable in terms of a single parameter, the 
vibrational temperature Ty. 



The Landau-Teller relaxation equation is a remarkably simple result con- 
sidering the internal processes that it describes. The question arises as to how 
accurately the actual physical processes are represented. In studies by Rubin 
and Shuler (75, 76) and Montroll and Shuler (62), assumptions (1), (2) and (3) were 
imposed and for the case of oscillators subjected to an impulsive change in heat- 
bath temperature, a solution of the complete set of master equations was obtained. 
It was shown that, for molecules having a Boltzmann distribution initially, relax- 
ation occurs through a continuous sequence of Boltzmann distributions. These 
results add validity to the Landau-Teller theory. Early experimental studies 
using shock-tubes (e.g., Camac (20), Cottrell and McCoubrey (23), Losev and 
Osipov (56) have also yielded results in excellent agreement with the simple 
theory. However, more exact theoretical studies including the effects of molec- 
ular anharmonicity (Bazley et al. (6), Osipov and Stupochenko (67), Bray (15, 16), 
and Nikitin (65)) suggest that at higher temperatures relaxation occurs with a 


*Landau and Teller (49) applied the relaxation equation to the acoustic problem 
of sound dispersion in gases. Later Bethe and Teller (8) used the equation for 
studying vibrational relaxation behind shock waves. 
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smaller characteristic relaxation time. These studies tend to support the 
validity of the vibrational temperature concept (that is, a Boltzmann distribution 
of vibrational energy). More recent experimental studies of flow in nozzles also 
indicate a departure from the Landau-Teller relaxation model (1, 2, 43, 44, 45, 

61, 78, 79, 80). Additional complication enters the problem when the molecular 
concentration can no longer be considered dilute (74). One must then include the 
effects of vibration-vibration (V-V) energy interchange (15, 16, 66, 74, 88, 92, 

92a, 101), which also lead to departures from simple relaxation theory (15, 16). 

On the basis of the foregoing studies it may be concluded that the Landau- 
Teller relaxation model (8, 49) accurately represents the processes occurring 
only when the upper vibrational states are not appreciably populated. When 
large populations of the upper states exist, modification of the theory is required 
to account for nonuniform energy-level spacing (4) and nonlinear vibrational 
transition rates (24, 37, 39, 65, 73). The major simplifying feature that the prior 
studies suggest might remain valid is the concept of a single vibrational temper- 
ature, that is, population according to a Boltzmann distribution (62, 81, 82, 83, 84). 

When the effects of dissociation and recombination can no longer be neglected, 
additional complication is introduced. The simple Landau-Teller theory for vi- 
brational relaxation no longer applies, since chemical effects by their very nature 
suggest a significant population of the upper states. These effects also have a 
perturbing influence on the population distribution, yielding other than Boltzmann 
distributions. If the characteristic time required to populate the upper energy 
states is comparable to or greater than the time required for dissociation, vibra- 
tional relaxation will also have an appreciable effect on the rates of dissociation 
(vibration-dissociation coupling) (Wray (99) and Rice (72a, 72b)). Finally, the 
possibility of coupling of both dissociation and vibrational excitation with rota- 
tional excitation should be considered (Bauer (3) and Bauer and Tsang (4)). 

The foregoing complications, not accounted for by the Landau-Teller theory, 
are summarized in the following list: 

(1) Molecular anharmonicity 

(a) nonuniformity of the vibrational energy-level spacing 



(b) nonlinear vibrational transition rates 

(c) multiple-level transitions 

(2) Non-Boltzmann population distributions 

(3) Vibration-rotation energy interchange 

(4) For nondilute molecular mixtures the effects of vibration-vibration 
energy interchange 

and, in addition, the chemical effects that occur according to the following 
mechanisms: 

(5) Dissociation via vibrational excitation 

(6) Dissociation via rotational excitation 

(7) Recombination occurring preferentially into the upper vibrational 
and rotational levels. 

Varied approaches have been adopted in theoretical studies of dissociation 
rates and of the transition probabilities that appear in the master equations. The 
variational theory of reaction rates utilizes the basic assumption that a reactive 
system can be described by the motion of a representative point in the phase 
space of the system (Keck (47) ). The procedure requires the computation of 
"trajectories" of these points to determine the net rate at which the represen- 
tative points pass through a particular surface in the phase space. By this 
method Snider (85)* obtained values of the effective dissociation and recombi- 
nation rates for the quasi-steady zone behind a normal shock wave (for a de- 
scription of this zone see Chapters III and IV), and showed that the ratio of these 
rates is the chemical equilibrium constant. In this work, Snider considered all 
effects in the above list except multiple-level transitions (lc). In an independent 
study, Treanor (94) obtained similar results. Montroll and Shuler (62a) con- 
sidered the problem of dissociation as a one-dimensional random walk (through 
the vibrational energy levels) with an absorption barrier (dissociation). The 
rate of activation is related to the mean passage time of the "walker." In other 
studies (Keck and Carrier (48) and Brau et al. (12)), the investigators introduced 



*The paper by Snider (85) also contains an excellent summary of prior research 
results, in particular, as applied to dissociation and recombination studies. 
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r a diffusion analogy for the master equations (effects (la, b) , (2) , (5) and (7) were 

accounted for). Vibrational excitation with eventual dissociation is thereby ex- 
plained in terms of energy diffusion through the vibrational energy states. The 
conclusions arrived at by this method are similar to those of the studies de- 
scribed above. These investigations relate primarily to the processes asso- 
ciated with dissociation. Benson and Fueno (7) using a deactivation theory, 
obtained results applicable for the inverse process of recombination. 

The above investigations have contributed to an understanding of the observed 
temperature dependence of experimentally determined rates, as well as indicating 
the types of population distributions to be expected. It is significant that deviations 
from a Boltzmann distribution are found primarily in the uppermost vibrational 
energy levels, where the influence of dissociation or recombination is greatest 
(see, e.g., 11, 12, 13, 48, 84, 94). These studies have been confined, for the most 
part, to the quasi-steady zone where the fractional molecular population in each 
vibrational energy state is constant and as a result solutions are relatively easy 
to obtain. In this situation integration of the complete set of rate equations is not 
required. 

Few solutions are available for the more difficult problem of integration of 
the complete set of master equations. Difficulties encountered in this task are 
well known (18, 29, 77). The importance of such solutions, aside from providing 
a description of the entire relaxation process, is in obtaining values for the 
"incubation" or "induction" time behind a normal-shock wave (11, 12, 13, 58), 
which is the time required to populate the upper vibrational energy levels, ob- 
served experimentally as a delay prior to the onset of dissociation (99). To 
explain this phenomenon, a manifestation of coupling of vibration and dissociation, 
Marrone and Treanor (58) introduced a relatively simple set of relations that were 
an outgrowth of equations used by previous researchers (36). The resulting 
formalism may be considered an extension of the Landau-Teller model, since 
the basic structure of "vibrational relaxation" follows Landau-Teller theory and 
additional terms are included to account for the effects of dissociation and 
recombination. 
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Despite its simplicity the Marrone-Treanor model contains most of the 
essential elements of vibration-dissociation coupling and has been widely used 
(27, 31, 32, 35, 41a, 57, 87, 89, also see discussion in 86). Their model is based 
on the assumption that the vibrational degrees of freedom may be represented 
by a single vibrational temperature. Dissociation preferentially from the upper 
vibrational levels is presumed, the dissociation probability being weighted in 
such a manner that it decreases exponentially for the lower energy levels (see 
Appendix C). The weighting factor contains a parameter that can be varied. As 
previously mentioned, the term that describes vibrational relaxation in the 
equations is formally identical to that of the Landau-Teller or Bethe-Teller 
model but involves different variables. In one version a truncated harmonic 
oscillator was used as a basis for defining the temperature dependence of the 
quantity q°° (93). Alternatively, q°° was based on the actual vibrational energies 
obtained from spectra (58). The "effective rate of dissociation" that results 
from their model contains a "vibrational coupling factor" that depends on the 
vibrational as well as the translational temperature. This factor is responsible 
for a delay in dissociation behind a normal shock wave. Thus the solutions for 
such flows do exhibit an incubation time prior to the onset of appreciable disso- 
ciation. Treanor and Marrone (90) suggest that their model is not completely 
satisfactory in that a strong temperature dependence occurs for the pre- 
exponential factor in the effective dissociation rate and this is not observed 
experimentally. Also, allowance is not made for other complications given in 
the earlier listing. In particular, nonlinear transition probabilities, deviations 
from a Boltzmann distribution, and rotational effects were not considered, 
although it is possible that an implicit inclusion of these effects may occur 
depending on the choice of parameters. 




In the present work the Marrone-Treanor model is generalized to include 
effects (1) and (2), that is, the effects of molecular anharmonicity and a non- 
Boltzmann distribution. This is done formally by starting with the master 
equations and using a procedure similar to that of Bethe and Teller (8). Rota- 
tional effects are not explicitly included. The formalism is sufficiently general, 
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however, that later inclusion of these effects may be possible when the associated 
rate parameters are better understood. The method follows, in principle, a 
suggestion by Carrington (21): "in systems of many levels, it may be useful to 
combine some levels into groups, and calculate average transition probabilities 
between groups. It may be that the difficulties of interpretation connected with 
the significance of the averages will be less than the difficulties caused by in- 
sufficient quantity or accuracy of the data in the system with the full number of 
levels." In the spirit of these remarks an effort has been made to hold to a 
minimum the number of parameters that are to be evaluated from experimental 
data. 

In Chapter II the model equations are derived by means of an energy-moment 
procedure operating on the master equations. It is shown that by separating the 
complete set of vibrational energy levels into two sub-groups and assigning a 
separate vibrational temperature to each sub-group, the complete set of master 
equations is reduced to a set of only four equations. The "reduced" equations 
are a complete set and are still general in that they include the effects of multi- 
level vibrational transitions. However, as in previous investigations (7, 8, 12, 

48, 49, 73, 75), the generalization is not continued throughout the paper because 
explicit representations for the various rate parameters are not available. 

In the later sections of Chapter II it is shown that further simplification of 
the equations can be made. An additional independent separation of the vibra- 
tional energy levels is introduced, below which the energy-level spacings are 
large and the transition rates small relative to the corresponding quantities 
above the separation. The procedure for evaluating these quantities is discussed 
in Appendix B. Essentially, the procedure involves consideration of the "rate of 
quantum transfer" resulting from molecular collisions. This quantity is evalu- 
ated and approximated by linear segments. Such a procedure allows an extension 
of the Landau-Teller theory in a systematic manner to include the effects of 
anharmonicity. The thermodynamic quantities associated with the vibrational 
mode are evaluated in Appendix A and it is shown that, with the method used to 
approximate the effects of anharmonicity, accurate evaluations of the 
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thermodynamic quantities result (partition function, energy, and specific heat). 

In the final section of Chapter II the dissociation and recombination terms are 
introduced. These terms are treated in a manner similar to that of Marrone 
and Treanor in that preferential dissociation is presumed. Only vibrationally 
excited molecules that belong to the "upper group" of vibrational energy states 
are allowed to dissociate. The model equations that result are listed at the end 
of the chapter. These equations are not greatly different in their outward 
appearance from those used by Marrone and Treanor (58). The major difference, 
aside from the appearance of two vibrational temperatures and the inclusion of 
anharmonic effects, is in the vibrational relaxation term. It is shown that the 
proper inclusion of the effects of truncation results in a term that is the same 
as the Landau-Teller relaxation equation except for a truncation factor. This 
factor is important when the vibrational temperatures are large. An interesting 
feature of the model equations is that they reduce to the Landau-Teller equation 
when the vibrational temperatures are not large. 



At this stage the model equations contain parameters that must be evaluated. 
In Chapter III tentative values of the parameters are assigned and solutions ob- 
tained for the quasi-steady zone behind a normal-shock wave. An iterative 
numerical procedure is developed for this purpose. A unique determination of 
all parameters does not result from this investigation but several general state- 
ments can be made. It is shown that when the level that separates the regions 
of different vibrational temperatures is placed about midway up the vibrational- 
energy well, a highly preferential dissociation is required to achieve the proper 
pre-exponential temperature dependence of the effective dissociation rate. Sep- 
arations lying above the midpoint energy yield results in disagreement with 
experiment. 

In Chapter IV numerical solutions of the model equations coupled with the 
equations of flow are presented. For flow behind a normal-shock wave, it is 
shown that the midpoint separation of vibrational energy states yields induction 
times that are about a factor of ten too large compared with experimental re- 
sults. It is concluded from this investigation that readjustments of the tentatively 
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chosen point of separation as well as other parameters of the model are required 
for agreement. In later sections of Chapter IV the model equations are investi- 
gated as applied to nozzle flow. It is shown that the characteristic vibrational 
relaxation time that the theory gives for such flows is qualitatively consistent 
with recent experimental results in that it is smaller by a factor of about 1/4 
as compared with results from Landau-Teller theory. Such reductions in the 
characteristic vibrational relaxation time, however, are realized only when 
freezing occurs at high vibrational temperatures. 

It should be emphasized that consideration of the complete set of master 
equations is complex and numerical solutions difficult to obtain (18, 77). A 
primary motive behind this investigation was the attainment of a smaller set of 
equations that retains the essential mechanisms of the complete set. The 
solutions demonstrate that the system obtained here can readily be solved 
numerically in conjunction with the flow equations. A numerical procedure 
developed by Lomax and Bailey (53) was found to be of great value for this 
purpose. Application of the method to the present problem is explained in 
Appendix E. 
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CHAPTER II 


DERIVATION OF THE VIBRATION-DISSOCIATION MODEL EQUATIONS 

In this chapter the basic equations required for a study of the effects of 
vibration-dissociation coupling will be introduced. It will then be shown in a 
formal manner how this rather large number of basic equations may be reduced 
to obtain a set of only four "model" equations. The formalism requires separat- 
ing the vibrational energy states into two groupings and then assigning a vibra- 
tional temperature to each group. The resulting system of equations still contains 
the effects of multiple transitions. In single-step transitions (vibrational "ladder 
climbing" model) it is also demonstrated how one may include the effects of 
anharmonicity by introducing additional separations in the vibrational energy 
states and still retain the relative simplicity of a Landau-Teller approach (8,49) 
as regards the form of the equations. 

In the first sections of the chapter only vibrational transitions not leading 
to dissociation are considered. Later the effects of dissociation are added and 
at the end of the chapter the equations to be studied further in the subsequent 
chapters are listed. 

II -A . Introductory Comments 

The mechanisms of vibrational energy interchange with translation (abbrev- 
iated T-V), of vibrational-vibrational energy exchange (V-V), and of dissociation 
(D) may be described by the following chemical equations: 

translation-vibration 

[A 2 ] v + M r - [A 2 ] v+i + M r (II-l) 

vibration-vibration 

[^2^v [A 2 ]j ±m [A 2 ]j ( n— 2) 

dissociation 


[A 2 ] v + M r = 2 [A] + M r 


(II-3) 



The symbol [A 2 ] y denotes a molecule (anharmonic oscillator) in the vth vibra- 
tional energy state; [A] symbolizes the atoms obtained by dissociation of the 
molecules [A 2 ]; and M r represents some species, specifically designated by 
the superscript r (r may denote atoms, molecules, or some chemically inert 
species, e.g., argon), that gains or loses kinetic energy and momentum during 
the reaction. Electronic states are excluded from the above description and it 
is assumed in the work that follows that any effect produced by electronically 
excited molecules may be assumed negligibly small. Such an assumption is 
perhaps not always valid, depending, of course, on the gas temperatures en- 
countered. Taylor (87) and Brau (13) have shown that the electronic states have 
some effect on vibration and dissociation, particularly at high temperatures, but 
that the effect is not large. 

To facilitate the description of the model a sketch showing the intermolecular 
potential of an arbitrary molecule is presented in figure II— 1. The various levels, 
numbered on the right-hand side of the potential curve, are the vibrational energy 
states attainable by the molecule. The lowest level is the zeroth state and the 
(N-l )th corresponds to the uppermost bound energy state. The 'Nth level defines 
the dissociation threshold and the dissociation energy. A molecule having vibra- 
tional energy exceeding this limit will have dissociated becoming a pair of atoms. 
In the model to be introduced the assumption is made that the vibrational energy 
states can be separated into two groups, a lower and an upper, designated by the 
letters A and B, respectively, as shown on the sketch. The separation is indi- 
cated by the dashed line between the levels b and b-1. Once this separation is 
made, it is further assumed that the distribution of energy in each individual 
group is Boltzmann and hence expressible in terms of parameters that are char- 
acterized as the vibrational temperatures of the separate groups. Further, only 
the B-molecules will be allowed to dissociate. Additional assumptions are also 
included that are not severe so far as the construction of the model is concerned, 
but are introduced in order to retain a form of the resulting equations that con- 
tains the fewest parameters possible. These assumptions are as follows: 

1. Dissociation occurs preferentially from the B states, that is, dissocia- 
tion from the higher levels is exponentially weighted over those from 



Dissociation 

energy 



the lower levels (Marrone and Treanor (58)). This is indicated by the 
arrows that are broadened as their length decreases. 

2. The effect of anharmonicity is accounted for by introducing another 
separation parameter a. The energy levels are taken to be equally 
spaced above and below this level and the rates used are harmonic- 
oscillator rate relations that vary linearly in terms of the vibrational 
level number. The spacing for the levels above the separation, however, 
is much reduced and the associated rates much greater than the corre- 
sponding quantities below the separation (see Appendix A and B). 

In applying the second assumption the effect of anharmonicity is actually accounted 
for by approximating the "rate of quantum transfer" (see Appendix B), which oc- 
curs as a result of vibration- translation energy interchange, with linear relations. 
In this thesis the approximation is carried out only for the case of "stepwise" 
transitions and with only two linear relations. 

The significance of the separation b, the boundary between groups A and B, 
is exemplified further in figure II-2. If one were to plot an actual population dis- 
tribution, represented by the ratio n v /n, in terms of a parameter 9 V , defined 
by the equation 

n v /n = exp(-E v /k0 v ) (II-4) 

for some arbitrary point in the region behind a normal shock wave (Treanor (94), 
and Keck and Carrier (48)), one would obtain the continuous curve drawn in the 
figure. The ratio n v /n represents the fractional population of molecules having 
energy associated with the \th vibrational energy state. The significance of the 
parameter 9 V is that, in terms of this parameter, a Boltzmann distribution 
would be represented by a horizontal line. The continuous curve in figure II-2 
illustrates the fact that the mechanism of dissociation perturbs the actual dis- 
tribution from being Boltzmann (48, 84, 94). The model approximates the actual 
non-Boltzmann distribution by two Boltzmann distributions that are described by 
the upper and lower vibrational temperatures, T A and T-g, respectively. 

The separations a and b will allow for a generalization of Landau-Teller 
theory to include, respectively, the effects of anharmonicity and the effects of 
non-Boltzmann-like population distributions. The assignment of values to these 




quantities is the subject of the later chapters and it suffices to state now that the 
b-level separation will always have a value that exceeds or at least is equal to 
that of the a-level separation; that is, 0 < a ^ b < N. 

The equations for the model are obtained in a manner analogous to that used 
by Heims (38) in that energy-moment equations of the population distribution are 
derived. The procedure deviates from Heims’ since the form of the population 
distribution among the vibrational energy states is assumed and hence not as 
many moment relations are required in order to obtain a complete set of equa- 
tions. The zeroth- moment equations denote the number of molecules in the 
lower and upper energy groups (i.e., groups A and B), and the first-moment 
equations give the average vibrational energy per diatomic molecule for each 
group. Four nonlinear ordinary differential equations thus result to describe 
the model molecule. 


The basic equations from which the model equations will be derived are 
described in what follows. The general equations that represent the net popu- 
lation gain or loss per unit time and volume for molecules in the vth vibrational 
energy state (commonly referred to as the "master equations") are given by the 
sum of the following partial contributions. The contributions follow according 
to the interaction processes described, respectively, by the chemical equations 
(II- 1) through (II-3). 
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where 


v = 0, 1, N - 1 


and 
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m r 
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E v 
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net number density of collision partners ^m = 

number density of collision partner of kind r 

number of atoms per unit volume 

probability per second per unit concentration of oscillators that 
a collision with species r will result in the vibrational transition 
from state i to state j 

probability per second per unit concentration of oscillators that 
a collision of a molecule in the ith state with another molecule 
in the Zth state results in the transitions i —* * j, Z — ■ k, 
respectively 

probability per second per unit concentration of oscillators that 
a collision of a molecule in the \th vibrational state with a 
species of kind r results in dissociation 

probability per second per unit squared concentration of atoms 
that a pair of atoms as a result of a three-body collision with the 
third body being the species of the kind r results in a molecule 
in the vth state 

vibrational energy of a molecule in the vth state (= E 0 v (1 - x 0 v 
+ yov 2 ), see e.g., Herzberg (40)) 

translational or gas-kinetic temperature. 



One notes that equations (II-5) and (II-6) involve vibrational transitions between 
bound states only. The terms in the summation over Z in equation (II-5), and 
over & and m in equation (II-6) represent the contributions due to one-step 
transitions, two-step transitions, and so on. These equations are thus general 
in that they include all possible transitions between vibrational energy states. 
The solution of the set of equations given by 



<n-8) 



coupled with the fluid-dynamic equations and the atom-conservation equation, 
would constitute, in principle, a solution of the complete problem of vibration- 
dissociation coupling. As will be pointed out later, such solutions are difficult 
to obtain even numerically (18,77) and have questionable value since the required 
transition probabilities have uncertain values. 

The most often practiced procedure in studies of such equations is to: 

(1) Assume that only single-step transitions are important (although 
dissociation may still involve a multistep transition), and thus only 
the first terms in the sums over & in equation (II— 5) are retained; 

(2) Consider only problems in which the molecular concentration is 
sufficiently small in a medium containing a high relative concentration 
of inert diluent that the nonlinear vibration-vibration exchange reactions 
(eq. (II— 6) ) may be neglected. 

These restrictions are generally valid (Rich and Rehm, 74). In the analysis that 
follows, however, the first condition will not be applied immediately in order to 
increase the generality of the equations that are derived, but will be applied later 
for the system of equations that are solved. The second restriction will be 
adopted since only relatively dilute molecular systems are considered in this 
thesis on the grounds that the inclusion of equations (II-6) make the resulting 
equations too complex for this study. 

In the derivations that follow, each term in equation (II— 8), representing 
equations (II-5) and (II-7) (excluding the term subscripted V-V), will be investi- 
gated separately. On the basis of the declared assumptions these terms will be 
simplified and the equations that result, containing few parameters, will then be 
the equations for the model. 

II-B. Development of the Equations 

It will be very helpful for the analysis that follows to make the notation more 
compact by introducing a function g^ + ^ y , defined as follows: 

g v+i,v = mr k v+J?,v[ n v+£ - n v ex p( kT )j (II_9a) 

= m r ,v n v+i - ky,v+i! n vj (II - 9b) 
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This function denotes the net rate of change of the population density of the vth 
state as a result of transitions to and from the (v + £)th state when these tran- 
sitions occur from translation- vibration energy interchange involving a species 
of kind r. Further, to avoid a cumbersome notation that is not immediately 
relevant to the derivations, the effect of the superscript r will be ignored 
initially. The effect of r can be neglected here since: (1) the kind of collision 
partner causing the vibrational excitation affects the equations in a linear manner. 
The final equations may be readjusted to include r simply by superposition of 
the resulting equations (i.e., by reinserting the superscript and reintroducing the 
summation over r); (2) only dilute molecular mixtures are of interest in this 
study; the effects of considering other than the diluent species (e.g., argon) as 
the collision partner will be small. 



The introduction of the g-function greatly simplifies the notation required for 
equations (II-5). Substituting appropriate values for the subscripts in the relation 
given for the g-function above, we obtain the following simple form: 


r dn vi 

L dt J r j’_'y 


N-v-l 

i 


£=i 


Sv+i,v 


V 

Z s v,v-i 
£=1 


( 11 - 10 ) 


(II-B-1). T-V Energy Interchange, Multiple Transitions. 

(II-B-la), Zeroth Moment Equations 

We now seek a relation that represents the net rate of change of the total 
number of molecules in the A-group of energy states. This is easily done by 
summing equation (11-10) for all the vibrational levels, v, such that 0 ^ v ^ b - 1. 
The equation that results, called the zeroth energy- moment equation for A-states, 
is 

v 

Z Sv,v-£ (E" 11 ) 

0=1 

( 11 - 12 ) 


1-dnA] b-1 N-v-l b-i 

L dt L v " L L Sv+£,v - L 
i “ v v=o v=i 


where 


b-i 

n A = E n v 
v=o 
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Certain elements in the second double summation on the right-hand side of 
equation (II- 11) may be combined with terms in the first double sum. The proce- 
dure for combining these terms may be more easily followed after introduction 
of figure II-3. Each mesh point in this figure represents a pair of coordinates in 
a two-dimensional coordinate space in which one coordinate represents the vibra- 
tional quantum number and the other denotes the number of steps involved in a 
transition. It should be pointed out that transitions involving the points in the 
triangular region above the line v - N-i-1 are not allowed in (T-V)-type reactions 
since the transitions associated with such points would involve dissociation. Such 
reactions are contained in equations (II-7). The mesh points in the crosshatched 
region correspond to transitions that result in a crossing of the separation bound- 
ary between the states b-1 and b, shown in figure II— 1. The first term on the 
right-hand side of equation (11-11) is a double summation involving the g-function 
evaluated at each mesh point in the closed region ABCDA. The second sum in- 
cludes values of g at each point in the region EFDE. We now rewrite equation 
(11-11) to obtain an expression having summations with identical arguments so 
that the g-functions can be related. This is done by noting that each element in 
the first sum is to the left of the boundary defined by the line v - N-i-1 and each 
element in the second sum is above the line v = i. Equation (11-11) may then be 
written 


~ dn A~ | 

dt 


IT-V 


N-b b-l N-l N-4+l b-i b-i 

Z Z ^v+i.v + Z Z *=v+i,v Z Z ®v,v-4 


4=1 v=0 


4=N-b+i v=o 


4=1 v=4 


(11-13) 


The subscripts and summing interval on the third term may be altered without 
affecting the sum to give 

b-i b-4-i 

third term = £ £ g v+ ^ v (11-14) 

J0=1 V=0 


We now observe that the various double sums involve the same argument, and 
we may therefore combine the appropriate individual terms. This will be shown 
graphically. We recall that terms one and two are double sums involving g y+ ^ y 
evaluated at points in the region ABCDA. However, the third term involves the 
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Vibrational energy level number, 





A G N B 

12 3 N-b b-l b N-2 N-l 


Number of steps in transition, l 


Figure II-3. Diagram showing mesh points in a two-dimensional coordinate 
space where one coordinate is the vibrational quantum number 
and the other coordinate is the number of steps involved during 
a vibrational transition 
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same g-function, but evaluated at each point in the region AGHA. We may sub- 
tract these common terms and the terms that remain are those values of g v+ ^ y 
for each point between the pairs of lines v = b-& and v = N-4-1, but below 
v = b-1, that is, the points shown in the crosshatched region in figure II-3. 
Equation (11-13) thus becomes 


N-b b-i 




[M. v - T z ? v +e ,v + t T' ^ 

£=l v=b-£ £=N-b+l v=b-£ 


N-l N-i-l 

+ Z Z % v +£ 9 \ 

i=b+i v=o 


(11-15) 


where the first term involves the mesh points in the region DICD, the second in 
JKLMJ, and the third in NBON. Equation (11-15) may also be more compactly 
written as 



b-i N-v-i 

-E E 

v=o i=b-v 


g V+i ,v 


b-l N-l 


= Z Z 

v=0 £=b 


,v 


(11-16) 


The significance of equation (11-15) is seen if a few of the terms involved in the 
summations are displayed: 

M 

[_ dt Jt_v Sb * b " 1 + ( gb+ l, b -l + g b,b-2) 

‘ ** v“ ■ ^ 

1 step 2 step 


+ (Sb+2,b-i + ^b+i,b-2 + &b,b-3) + * * * - (11-1^) 

' 

3 step 

Each term represents the loss (or gain) in the total number of molecules in the 
A-group of levels as a result of the transitions, one-step, two-step, and so on, 
that cross the separation b. Although equations (11-15) and (11-16) are equiva- 
lent, the first expression, equation (11-15), involving the sum over i, is more 
advantageous in practice since the terms are ordered in a monotonic sequence 



of decreasing values, that is, as Z increases, each succeeding term is less 
than the preceding term (see, e.g., 24, 37, 39). The series converges so rapidly 
that in most cases only the first term is important. 


The equivalent sum representing the net rate of change of the total number 
of molecules with energy in the B-group of states is also readily obtained. As 
noted previously, equations (11-10) (or (11-15)) exclude transitions that result in 
dissociation. It then follows that the T-V interactions conserve the number of 
molecules; that is, 


_d_ 

dt 



v=o 


T-V 


0 


(11-18) 


This identity may also be proved by summing equation (11-10) for all v in the 
range 0 ^ v ^ N-l. From this equation it follows that 


r dn Bi r dn A ] 

L dt l dt Jrp_ 


T-V 


N-b b-i 

-E E 8. 

f-i v=b-£ 


v+i ,v 


b N-je-i 

E E * 

i=N-b-i v=b-Z 


v+Z , v 


N -Z N-i-l 

Yj Yj S v+ £ 9 y 

,0=b+i v=o 


(11-19) 


or 


rdn T ,'i b-1 N-l 

hri v - - 2 E ... 

T-V y=o Z=b 


( 11 - 20 ) 


where 


N-i 

n B = Yj n v (II-21) 

v=b 


(II-B-lb). First Energy- Moment Equations 

The differential equations that yield the change in the net vibrational energy 
(first-moment equations) for both the lower and the upper vibrational energy 
states are also required. These are obtained by first multiplying each equation 



of the set of equations (11-10) by the energy associated with the respective level, 
E v ; the resulting equations are then summed over all the levels associated with 
the particular group. 

Thus for the A-states we obtain 


jd_ 

dt 


b-i 

Z Ev Iv 


■ z 


V=0 


N-v-i 

z 


i=l 


E v S v +i,v 


b-i v 

Z Z E v Sv,v-i 

v=o £=l 


( 11 - 22 ) 


The order of the summations is now interchanged and requires that one keep 
track of the regions involved within the summations. We introduce the expres- 
sion equivalent to equation (11-14) but appropriate to the above equation, and 
there results 


d_ 

dt 


b N-i-i N-i N-i-i 

+ Z Z E vSv+i,v + Z Z E v%v+£,v 

i=N-b-i v=b-i £=b+i v=o 

b-i b=i-i 

“ Z Z E v+£&v+£,v (11-23) 

£= 1 V— o 


z 

|v=o 


E v n v 


b-i b-£-i N-b b-i 

= Z Z E v Sv+g,v + Z Z &v+£,v 

T_v -2 =1 v=0 i=i v=b-£ 


The first and last terms can be combined and equation (11-23) becomes 


_d_ 

dt 


b-i 

Z E v n v 

Ly=o __ 


T-V 


b-i b-£~i 

= Z Z ( E v “ E v+i)Sv+i,-\ 

£=i v=o 


N-b b-l b N-i-l 

+ Z Z E v Sv+i,v + Z Z E v Sv+i,v 

£=i v=b-£ i=N-b+i v=b-£ 


N-l N-i-1 

+ Z Z E v 

£=b+i v=o 


(11-24) 


Interchanging the order of summation and combining the last two terms, we 
obtain the simpler-appearing equation 



i 


d_ 

dt 


b-i 

Z E v n \ 

v=o 


b-l b-i-l 

= Z Z ( E V " E v+i)Sv+i,A 

X-V i=i v =o 

b-l N-l 

+ Z E v Z 

v=o i=b 




(II-24a) 


The last term in this equation yields the gain (or loss) in vibrational energy as 
a result of transitions that shift the molecule from the A- to the B-group of 
energy states (this term is analogous to the terms in eq. (11-15) or (11-16) ), 
while the first term represents vibrational-energy relaxation within the A-group 
of states (note the summation region). This point will be made clear in a later 
section. 

The comparable equation that relates the rate of change of net vibrational 
energy per unit volume for the upper energy states can be obtained by first 
deriving the energy-relaxation equation for the complete molecule and then 
subtracting that portion, equation (II-24a), associated with the lower states. 

The equation for the complete molecule, obtained by multiplying equation (11-10) 
by E v and summing over all molecular energy states 0 ^ v ^ N-l, is 


_d_ 

dt 


N-i 

Z E V n V 


v=o 


N-l N-v-i N-l v 

= Z Z E v ^v+i,v “ Z Z E v ^v,v-^ 


(11-25) 


JT-V 


v=o £=i 


v=o £=i 


Interchanging the order of summation and again altering the indices and summa- 
tion scheme associated with the second term to obtain the same g-function as 
contained in the first term results in the equation 


_d_ 

dt 


N-l 

Z E V n v 


v=o 


N-l N-£-i N-l N-4-1 

= Z Z E v Sv+i,v “ Z Z E v+^Sv+i,v (11-26) 


Jt-v 


i~\ v=o 


£~i v=o 


The summation intervals are identical, that is, both terms involve the same 
mesh points; hence each term may be combined individually giving 


jd_ 

dt 


N-l 

Z E v*H 

v=o 


N-l N-jfl-l 

Z Z ( E v ~ E v+i)Sv+i,'\ 

r j'_y £=l v=o 


(n-27) 
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f 

Y Since dissociation is not contained in the translation- vibration energy- interchange 

equations (i.e., the (T-V) -subscripted equations), the above equation character- 
izes vibrational energy relaxation for the molecule (compare with the first term 
of eq. (11-24) ). Subtracting equation (11-24) from the above relation yields the 
equation representing the rate of change of vibrational energy for the upper 
states as follows: 


N-l 

Z E v n v 

ii 

N-i 

Z E v n v 


_d_ 

dt 

b-i 

Z E v n v 

v=b _ 

T-V 

V=0 _ 

T-V 


v= 0 


N-i N-i-i 

= Z Z ( E v “ E v+^) ®v+i,v 

£=0 V=0 


b-i b-£~i 

- Z Z (Ev E v+i ,v) Sv+i ,v 

v=o 


N-b b-i b N-i-i 

- Z Z Sv+i,v “ Z Z E v§v+£,v 

£=i v=b-i £=N-b+i v=b-i 


N-i N-i-1 

- Z Z EvSv + <,v V-W 

i=b+l v=o 


The first and second terms may be rearranged and we obtain 


_d_ 

dt 


N-i 

Z E V 


v=b 


n v 


T-V 


b-i N-i-l 

~ Z Z ( E v - E v+i)Sv+i,v 

£=1 v=b-£ 


N-l N-i-l 

+ Z Z ( E v - E v+i) Sv+^.v 

i=b v=o 


N-b b-i b N-i-i 

- Z Z E vSv+i,v Z Z E vSv+i,v 

l - 1 v=b-£ j!=N-b+i v=b-£ 

N-i N-^-i 

Z Z E v&v+i,v (11-29) 

&=h+i v=o 
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or, after combining terms. 


in 


j 


d_ 

dt 


N-i 


Z E v 

v=b 


n-i 


JT-V 


b-i N-.0-1 

Z Z ( E v ~ E v+ ^) g v +i,v 

1 = l v=b-£ 

N-l N-i-i 

+ Z Z (Ev - E v+ ^)g v+ £ ;V 

Z=b v=o 


b-i N-l 

- Z EvZ g^,v (n- 29 a) 

v=o b 


The first two terms, associated with the B-group vibrational relaxation, involve 
evaluating the function (E v - E v+ ^)g v+ ^ jV at the points between the lines 
v = b-£ and v = N-,0-1 in figure II-3. The change in vibrational energy due to 
transitions that involve a crossing of the separation boundary between the A- 
and B-groups are given by the last three terms in equation (11-29) or the last 
term in equation (II-29a). 



(II-B-lc). Vibrational Temperatures 

The essential equations required to specify the vibrational- energy relax- 
ation processes for a molecule still excluding the effect of dissociation, are 
given by equations (11-15), (11-19), (11-24), and (11-29) (or equivalently eqs. 
(11-16), (11-20), (II-24a), and (II-29a) ). These equations are a complete set, 
provided the population distribution within each group of energy states is known. 
As pointed out earlier, we assume that this distribution can be reasonably rep- 
resented by a Boltzmann distribution, that is, by 

(ny+Oi = (rMj exp V+ (ZZ Ey ^] (H- 30 ) 

where i denotes the group, either A or B, and T^ is the temperature appro- 
priate to the group. The validity of such an assumption has been shown by 
Rubin, Montroll, and Shuler (62, 75, 76, 81, 82, 83, 84) for single-step transi- 
tions only for which the transition probabilities are those of harmonic oscillator 
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molecules. For such conditions we have T A = T-g. Little is known about the 
effects associated with multiple transitions as a perturbing influence on the 
distribution. Bazley et al. (6), while solving a problem comparable to Montroll 
and Shuler’s (62) but including the effects of anharmonicity and of some multiple 
transitions, has shown that the relaxation processes are comparable with those 
obtained for harmonic oscillators. Thus, existing theories appear to validate 
the use of the Boltzmann distribution except when dissociation becomes impor- 
tant (12, 48, 84, 94). The perturbing influence of dissociation is, of course, 
being considered in the present model by the separation of the energy levels 
into two groups with separate temperatures T^ and T B . 

Assuming a population distribution, in our case a Boltzmann distribution, 
allows one to reduce manyfold the number of dependent variables required for 
a description of the vibrational mode. Ultimately, the problem of determining 
the vibrational state of the molecular system will be reduced to one of obtaining 
values of the four variables n^, ng, T^, and Tg, where n^ and n 3 are defined 
by equations (11-12) and (11-21) and T^ and Tg are obtained on the basis of 
equation (11-30). The population density in any vibrational level can be expressed 
in terms of the total group population by the equations 



n A 

Qa( t a) 

n B 

Qb( t b) 

n B 

%( T B) 



(II-31a) 
(II-31b) 
(II- 31c) 


where the energy of the vibrational level is subscripted to designate its group, 
that is, 


4* 

II 

< 

0 < v < b-1 

: B V = E v 

N-l 
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The A and B subscripts may be dropped if no ambiguity results. The functions 
Q^(T) and Qb(T) (or Q'g(T) ) are the partition functions associated with the A- 
and B-molecules, that is, the A and B grouping of molecular states, and are 
defined by 


Qa(T) = exp (-§) 


N-i j £ > 

Qb< T ) = E exp (”kT 

v=b \ 


(II-32a) 


(II-32b) 


/ Eu ^ 
- QbCD exp(-^ 


(II-32c) 


The partition function differs from Qg by the reference energy. t 


Simplification of the equations is now possible. Substitution of the g- 
function, equation (II-9), into equation (11-16) yields 





b-i N-i ^ b-i N-i ^ 

m Z I ki, v n j ~ m Z Z kg,V n V ex P 

v=o i=b v=o £=n 

(11-33) 



We observe that the left term involves only n^, where b ^ i ^ N-l, and the 
right-hand term only n v , where 0 ^ v ^ b-1, that is, the left-hand term con- 
tains only B-state population densities and the right-hand term only A-state 
densities. In view of equations (11-31), equation (11-33) may also be written 


~ dn A~ 

. dt J T-V 




0I-34a) 


t The energy levels in and Qg are referenced with respect to the molecular 
ground state; in Qg they are referenced with respect to the lowest energy level 
of the B-group of energy states (see Appendix A for discussion). 
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and the rate equation for n B , equation (11-20), is similarly given by 


-m[n B kR } (T B ,T) - n A k ( F } (T A ,T)J (0-34b) 


In both of these equations we have defined the following rate constants: 


b-i N-i 




(E< - E b ) 


(11-35 a) 


, . b-l N-i (” /E-d Ea \ Ea 

* (l) 1 V V ? V a z v) 

kp (Ta ’ T) " «a( t a> l t ’ v<T> " P l" “ “a. 

(II— 35b) 

Similar relations can be obtained for the energy rate equations. We first 
introduce the energy averages defined by the barred quantities given in the 
following identities: 


b-i N-i 


(Tb.T) k^(T B ,T) = q^ T b ) E e v E k i)V (T) exp|- 


v=o i=b 


( E B, - E b) 


(n-36a) 


b-i N-i 


' (T A >T) ^ F (T A >T) = oTTtT) E Ev E k JfV (T) 


Qa(t a ) 


exp - 


( E Bi " E A v ) E^ 


ffl-36b) 


Substituting these relations and equations (11-12) and (11-21) into equations 
(II-24a) and (II-29a), we obtain the energy equations 


I;Ma(Ta)1 


b-i b-i-i 

E E (E v - E v+J g)g v+J g A 


i=l v=0 

— ^ (l) — A (l) 

+ mn B E R (T B ,T)k R (T R ,T) - mn A E F (T A ,T)k F (T A ,T) 


(H-37a) 
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E v+i)Sv+i ,v 


d b-i N-i-i 

*[^(T B jL = Z Z (E v - 

J1 ~ v i=i v=b-^ 

N-l N-J-l 

+ Z Z ( E v “ E v+i)gv+i,v 

i=b v=o 

- m ng Ep (T b ,T) (T B ,T) 

- * a) 

+ mn A E F (T A ,T)k F (T A ,T) (n-37b) 

where 


b-i 

n A^A(^A) “ Zj ^v n v (II-38a) 

v=o 

N-i 

n B^B(^B) “ Zj ^v n v (II-38b) 

v=b 

= n B [^B ( T B) + E b] (II-38c) 

A brief discussion of these equations is worthwhile. The population equations 
(eqs. (II-34a) and (II-34b) ) are precisely analogous to the master equations from 
which they are derived. Each equation contains only two terms, one represent- 
ing the rate of loss and the other the rate of gain. The rate coefficients in these 
equations contain the net effect of all possible transitions whether they be one- 
step, two-step, etc., or any combination thereof. The rate coefficients also de- 
pend on the vibrational temperatures of the respective grouping of levels. These 
temperatures may, in principle, be obtained from equations (11-37). The func- 
tions q^ and qg are known single-valued functions of their arguments (see,e.g., 
Appendix A) and once the value of these functions has been found, one may read- 
ily compute the temperatures T A and T B . It is difficult to discern the signifi- 
cance of the individual terms contained in equation (11-37) at this time. It will 
be shown in the next section that the first term containing the double summation, 
is proportional to the difference q^ (T) - QaC^A) i* 1 equation (II-37a) and 



jg (T) - qjj (Tg) in equation (II-37b). For this reason they will be called the 
vibrational-energy relaxation terms (they are analogous to the Landau-Teller 
expression (8, 49) that describes vibrational relaxation). However, one can de- 
rive such explicit relations only for single-step transitions and only if the tran- 
sition rates are linearly related to the vibrational level v. 

The present problem requiring solutions of equations (11-34) and (11-37), 
although less complicated than the original problem in which one seeks a solu- 
tion of the complete set of master equations, still contains far too many param- 
eters ky^j for easy evaluation (i.e., if one is to compare solutions of the 
equations with experimental data as a means of evaluating these parameters). 
These equations are also too complex to be useful as a model for flow field 
calculations. It will now be shown that further simplification is possible after 
the introduction of additional approximating assumptions. 

(II-B-2), T-V Energy Interchange, Stepwise Transitions 

The formalism used and the equations obtained in the previous sections are 
general. The only limitation in the final form of the equations concerns the 
assumption of the population distributions. Although general expressions are 
much sought after, in practice such expressions often have only limited utility. 
For example, a solution of the problem outlined in the last paragraph of the 
previous section would probably have only limited value, since it would be diffi- 
cult to relate the parameters j to experimental results. Rather than solve 
this problem, additional approximations will be introduced here that will result 
in a final set of vibrational -relaxation equations that are only slightly more com- 
plex than those obtained by Landau and Teller (see, e.g., eq. (1-1) in Chapter I). 

The significant assumptions to be introduced in this section are as follows: 

(1) Vibrational relaxation occurs in a stepwise manner, that is, the effect 
of multiple transitions may be assumed negligible (k^j = 0 for 

|i - jl = £ > 1). 

(2) The effect of anharmonicity regarding both the thermodynamic functions 
and the transition rates, may be approximated by use of equal energy- 
level spacing and harmonic oscillator transition-rate relations above 
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and below some separation a, where large energy -level spacings and 
small transition rates are used for the lower levels and much smaller 
spacings and much larger rates are used for the upper levels. 

It is worthwhile to remark briefly on these assumptions. The first is generally 
valid, except perhaps for the uppermost vibrational energy levels of a molecule 
(24, 37, 39, 67). For the uppermost levels, however, the probabilities for mul- 
tiple transitions from any given level are still less than the probability for a 
single-quantum transition from that level (24, 37, 39). As a matter of interest, 
it will be shown in the next chapter that if the separation b is taken near the 
midpoint of the vibrational potential well, the observed pre-exponential temper- 
ature dependence is best accounted for by a dissociation model that has disso- 
ciation occurring primarily from the very uppermost level. A nonpreferential 
model for dissociation that presumes equal probability of dissociation from any 
of the B-group of vibrational energy levels does not yield results that agree 
with the experimentally observed pre -exponential temperature dependence. 
Hence, there is an apparent consistency between the use of assumption (1) and 
an equivalent assumption for the dissociation model that will be introduced in 
a later section. 

Assumption (2) is discussed in Appendix B. It will only be remarked here 
that the approximation procedure can be further expanded. The approximation, 
as used, depends on the use of two linear functions for approximating the ’’rate 
of quantum transfer” that results from single-step transitions (see Appendix B). 
Including more linear functions and thus additional separation levels would, in 
principle, improve the model equations. It will be shown, however, that the 
additional complication may not be warranted as a result of the insufficient 
accuracy of the available experimental data (considerable "scatter" in the 
data). As more refined experimental measurements become available, addi- 
tional refinement of the model equations may also be required. The formalism 
allows for such improvements. 

As a matter of emphasis, it is worthwhile to point out again that the reason 
for introducing the a-level separation is quite different from that for the b-level 
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separation. The b-level was introduced to take account of the fact that the 
distribution of energy in the various vibrational states is not Boltzmann. The 
a-level allows an approximation for the molecular anharmonicity effects, both 
as regards the energy-level spacing and the transition probabilities. The 
assumptions given earlier will be introduced in turn in the analysis that follows. 


We find on imposing (1) that the resulting equations are greatly simplified. 
Equations (11-35) and (11-37) become 


(l) k b,b-i * T) 

k R( T B. T ) = ^7 


(II-39a) 


k F ( T A>T) 
dt [ n A$A( T A)] T _ v 

<5t [ n B ( t b)] t _ v 


k b-l,b l E b-i \ 

Qa(Ta) eXP \ kT A/ 

b-2 

Z ( E v - E v+i)gv+l,v + 
v=o 


N-2 

Z ( E v - E V+i)g V +l,V 
v=b-i 



(II-39b) 

(II-40a) 


ai-40b) 


In the latter two equations use has been made of the resulting simplification 

that occurs for the averages in equations (11-36), which are now given by 

E F = Ep> = E| 3 _ 1 . Equations (11-34) are unchanged after the application of con- 

-(l) -(l) 

dition (1), except that the rate constants k F and k F are given by equations 
(11-39) rather than by equations (11-35). 

As a result of imposing condition (1), the rate equations (11-34) may be 
expressed simply as the first term of equation (11-17), that is, by 


~ dn A~ | 

_ dt T-V ^b’b-i 



(II-41a) 

(II-41b) 
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With these relations it is a simple matter to remove the explicit dependence of 
equation (II-40b) on the rate of change of the total number of molecules in the 
A-states. This we do by shifting the B-state reference energy by substituting 




■A 


the appropriate expression for q g in terms of q-^, equation (II-38c), and then 
combine terms after using equation (II-41a) for (dn^/dt)«p_y and equation (II-9) 
for the g-function. The alternate form for equation (II~40b) that results is 

N-2 

[ n B^B( T B)] = Tj ( E v " E v+t)Sv+i,v (11-42) 

L JT “ V v~b 



The summation contains no A-state quantities. The reason for making such an 
alteration is that all terms in the series in this equation involve only B-state 
quantities and thus are to be evaluated at the vibrational temperature T-g. This 
simplifies bookkeeping in the derivations that follow. 

Further simplification of equations (II-40a) and (11-42) now requires the 
introduction of assumption (2) given earlier. In addition, a procedure is re- 
quired that is similar in many respects to that used by Bethe and Teller (8), 
described also by Vincenti and Kruger (97) and Clarke and McChesney (22), 
except that rather than consider a molecule that is a harmonic oscillator with 
an infinite number of energy levels we have one that is truncated and thus has 
only a finite number of levels. The procedure depends on the fact that, pro- 
viding the energy levels are equally spaced, substitution of transition proba- 
bilities varying in a linear manner with the level number v, yields series 
expressions in the energy rate equations that are summable and may be identi- 
fied with previously defined functions. The procedure depends on the following 
description for the level energies E v (a diagram of these levels is given in 
fig. A-l): 


E Ai V 


ri 

VI 

> 

VI 

o 

1 

(II-43a) 

E a + Ea 2 (v 

- a) 

i 

rO 

VI 

> 

VI 

1 

(II-43b) 

E b + e A 2 < v 

- b) 

b ^ v ^ N - 

1 

( II-43C) 
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where 


E a = E A* (a - 1) + E a? (H-44a)t 

E b = E a + (b - a) (II-44b)t 

In addition we require the transition probabilities (see Appendix B) given by 


*w,v-i 


= vki / 0 

= (v - a)k$ + ka,a-i 


— exp 
= exp 



k V,V-l 

k v,v-i 


'(2) 

^a.a-i 


- *k (l) 

' ea. 


K a 


i,a 


= exp 


(-£)> 


r (2) 

-a,a-i 


1 < v ^ a - 1 

a < v ^ N - 1 

1 < v < a _ 1 

a < v < N - 1 


(II-45a) 

(II-45b) 

(II-46a) 

(II-46b) 

(11-47) 

(11-48) 


The energy parameters E a , E^, and N and the energy-level spaeings E A 

and E a are formulated and discussed in Appendix A. The rate parameters, 
wi) 2 w 2 ) 

ko/i and kg ?1 are discussed in detail in Appendix B. We will presume for the 
present that the above relations are valid representations of their actual equiv- 
alents and proceed to reduce equations (II-40a) and (11-42) to the simplest form 
consistent with the basic assumptions (1) and (2). 

In what follows the key to successfully reducing the equations to a form 
that involves known functions and very few parameters depends on an effica- 
cious grouping of terms. To this end, rather than work with equations that 
are extremely bulky and contain many terms that can overly complicate the 


TThe subject of the following chapters involves in part the assignment of appro- 
priate values for these parameters. It suffices to state now that 0 < a ^ b < N, 
that is, the b-level, has a value that is above or, in the least, is identical with 
the value of the a-level. 





derivations, symbols will be assigned to certain individual terms. These terms 
will be separately investigated and simplified. One will find that in addition to 
greatly reducing the amount of algebra much insight is also acquired as to the 
expected form of a more general set of equations that would, in principle, even 
better approximate the actual relaxation processes. Such generality, however, 
involves the introduction of more parameters into a problem already plagued 
with too many difficult to evaluate parameters. Such generalization will not be 
explicitly exhibited here but the generalizing method will be made clear. 




The first term will be denoted as T t and the second as T 2 . There is little 
difference between equation (11-42), T t or T 2 , except for the summation limits. 
The difference is significant for Tj only because there are no vibrational 
states below v = 0. The algebra required to reduce these relations to their 
simplest possible form is least for T 1? and therefore attention will be focused 
first on this term. The complexity involved in reducing equations (11-42) and 
T 2 is nearly equivalent. However, it will be worthwhile to investigate equation 
(11-42) rather than T 2 , because more insight is gained as regards the general- 
ity of the method. The resulting equation immediately provides the proper 
expression for T 2 by simply reassigning different values for the summation 
indices b and N and by appropriately changing the subscripts. These param- 
eters, of course, will be implicitly imbedded in the resulting equations. As a 
matter of emphasis it is fitting to point out again that the key to the successful 
application of the following procedure is in the recognition of sets of terms as 
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being previously defined functions. These functions are developed and discussed 
in detail in the last section of Appendix A. Being familiar with these functions 
is important, particularly those relating to the vibrational mode, that is, Q^, 
Qa 2 > Qb> q Al » ^A 2 > ^B* ^ A 2 * ^B anc * tbeir primed equivalents. These quan- 

tities will be reintroduced as required in the analysis that follows to facilitate 
the recognition of specific terms. 

Concentrating now on T 1# we substitute equation (II-9) giving the g-function, 
equation (II-31a) for the number density n v , and, after the uppermost term in 
the series is adjusted, there results 


T l = n Aj 


m 


a-i 

Z 


( E v-i E v) 


, / E y \ a-i (Ev _ Ev+i)r / E v 

kv,v-,exp(- kTA j-^ QAj(Ta) k v,v+l«*p(- kTA 


«A,(T A ) 
(^a-i “ E a) 


J a-i 


Q^0~ K a-.,a=xp^ 


(11-50) 


As already pointed out, the partition functions Qaj(Ta) as well as the other 
thermodynamic quantities to be used are defined in Appendix A. The appropri- 

A 

ate expressions for E v and kjj , equations (II-43a), (II-45a), and (II-45b), 
respectively, are now substituted into the above equation and (k[|g - kg||) is 
removed as a factor. One then obtains 


T i = “ n Aj 


m (k!jo - ko!i) 


a-i 


e A,v / E Aj v 

eXP \ kT A 


J'.Qa.CTa) 


/ e a , v\ 

„ 1 y exp \ ~etx~ ) 

H*!?./! / £ 8\v=. Q A,(Ta) 

‘W 

(a-i) E Al ' 

exp - 
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+ E a a 


*o,i 


kT 


A J 


Al k (1) 

K i,o/ 


'k 0 W A q A!( t a) 


k (1) 

K l,0- 


(n-5i) 
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This expression is readily simplified. We observe in Appendix A that the first 
term in brackets is identified with the quantity 


E A t v ( e Aj v 

qA i (TA) ~ v 5 0 QA^TA) 6XP V kT 


(n-52) 


After substitution of the appropriate exponentials for the rate ratios, according 
to equation (II-46a), and it is noted again from Appendix A that 


V / E Al v\ 

Qa, ( Ta) - l ^{-VrZ) 


k ( ^ 

^ - E A t rg 

k 


E Ai 


(D /r(i)\ 
1,0 / K o,i \ 

~Wd) / 

Xk i,o/ 
(-£) 
(-^) 


1 - exp 


(n- 53 ) 


(II-54a) 


(II- 54b) 


the second and third terms may then be readily identified. There results 




r r Ca-DEAni 

(k f w _ kjwY 

^(Ta) - ^(T) 

, " kT A J 

\ k m K o,iy 

Qa^Ta) J 


( (n-55) 


the factor m(kj^ - k^j), which may also be written mkjjfl - exp(- - )] , is 

defined as the inverse of the vibrational relaxation time and is commonly given 
the symbol 1 /t (see, e.g., 22,97). To differentiate this parameter from a 
similar parameter that will be introduced later, we adopt the notation 


t A a < T ’P> 



(11-56) 


Furthermore, referring again to Appendix A, we observe that the factor given in 
the brackets of equation (11-55) above is the truncation factor defined by 


(n-57) 


i a 1 ( t a) 


a exp 

(a-i) E Al 


L" kT A J 


QA t (T A ) 1 


The characteristics of this function are described in detail in Appendix A. Sub- 
stituting these latter two quantities into equation (11-55) results in the following 
simple expression for the first term of equation (11-49): 


Ti = ^[^(T) X Al (T A ) - q Al (T A )J (H-58a) 

This expression is worthy of brief comment since it will turn out to be 
equivalent to the form used by previous workers studying vibrational relaxation 
in flow fields having relatively low temperature T. Noting equation (A30) in 
Appendix A, we rewrite equation (II-58a) as follows: 


Ai r^oo / \ n 

T ‘ ' ^K (T) " Oa.N W <n-58b) 

The truncation factor is approximately unity provided T A < (a - 1) E ^j/k. With 
this approximation (equivalent to the requirement that the population in the upper 
levels be relatively small or that the effect of truncation be negligible), the 
equation is then identical in form to that of Bethe and Teller (see, e.g., (8) or 
Chapter I) and is the equation used by many researchers studying the effects of 
uncoupled vibrational relaxation (10, 22, 29, 30, 32, 97). We will note later that 
the other terms in equation (11-49) and also the term in equation (11-42) are 
negligible in this case. The above equation shows the correct use of relations 
for truncated harmonic oscillators. To illustrate the differences between this 
relation and relations used by other researchers (57, 58, 90, 91, 93) the equation 
can be cast in a different form. After the introduction of equation (A30), the 
above equation may also be written 


Ti 


— 

r A 



•^(Ta) 

Wt) 


%( t a) 


(II-58c) 
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and this can be identified as the vibrational relaxation term used by the other 
researchers (see, e.g., Marrone and Treanor (58)), if the ratio (T^) /z^ (T ) 

is set equal to unity (for this comparison we consider a to be equal to N). 
The quantity q A is the truncated harmonic oscillator function that was used 
by the other researchers. This point will become clear as the other terms in 
equation (11-49) as well as the term in equation (11-42) are simplified. It should 
be evident, however, that the ratio Z A ^(T p^/ z^ (T) is not unity in the case 
studied by Marrone and Treanor; furthermore, the use of more accurate quanti- 
ties for the energy functions q^ (e.g., see 57, 58) is of questionable value. 
Further discussion of this point will be deferred until later (see last page of 
section B-4 of this chapter). 

Before proceeding with the simplification of equation (11-42) it will be help- 
ful to abbreviate the notation by temporarily assigning the symbol T 3 to the 

summation for — [n^q^ (T^)] , that is, 

dt T -V 

N— 2 

T 3 = ^ (E v - E v+ i)g VjV+ i (11-59) 

v=b 

The procedure to reduce T 3 to the simplest functional form parallels that just 
employed to yield equations (11-58). We obtain 


fN- 


T 3 = mn-g 


S 


1 ( E V-1 - E v) £ 

Qr-, k v. 


( E v\ N^1 (Ev _e v+i) 
exp - — : I- / , — - — — — : — 



( E b-i - E b) f ( E b\ (%-! - e n) r / %-iV 

Q b (T b ) b » b - 1 eXP \ kT B ) + Q b (T b ) exp \ kT B ) 

(n-60) 

as the expression analogous to equation (11-50). We note that the third term 
occurs here as a result of the non- zero lower summation limit. The fourth 
term is analogous to the third term in equation (11-58). The last two terms are 
required in order to obtain the appropriate summation limits for the set of 
energy states being investigated. Substitution of the appropriate form of the 

A 

energies E v and transition probabilities k^j, given by equations (II-43c), 
(II-45b) and (II-46b), yields 
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T 3 = 


fN-i t? 

n B m I Z 
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B 


(11-61) 


The denominator in each of the terms is Qg (T B ). It is helpful now to list the 
additional thermodynamic quantities that will be required here. They are also 
given in Appendix A and are as follows: 


Qb ( t b) 
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Z exp 
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(v-b) E Aj 
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(n- 62 a) 


(II-62b) 

(II-63a) 

(II-63b) 


(H-63c) 
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(11-64) 
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N- 1 E a (v-b) T (v-b) E a 

,b(Tb) - «mo - E b ^pgr 


(II-65a,b) 


We now combine the terms in equation (11-61) in such a manner that those con- 


^(2) ^(2) 

taining the coefficients k 01 and kJ j0 are grouped together. The set of terms 


designated T 3 may now be written 


rn _ f /gfe) I'M) V EA 2^ V ~^ 

* * ' B ( < M ' M [4 o5?rij- 


exp - 


(v-b) E a , 


~ + E^Cb-a) 


A 2\ r( 2 ) 


1 - —j 4 — exp - 

Qb( t b) 


(N-b-i) Ea 2 


q b ( t b) 


exp - 


(N-b-l) E A , 


+ “Ea (kj/a.t - ka-i,a) 


lv o — 1 fk 

+ E A m 4 . exp - 
* ®b( t b) L 


(N-b-l) E A n E Ag m kb,b-i 
kT B _ ^B ( t b) 


(11-66) 


To obtain this expression the partition-function representation Qjj (Tg) was 
substituted in place of its series as given by equation (II-62a). Certain of the 


terms above may be immediately identified. The first term is the quantity 


represented by equation (II-65b). The third and fourth terms together are the 
truncation factor £b(Tjj), equation (11-64), and their coefficient is given by 
equation (II-63c). We thus obtain 


T 3 n i$j m ( k l,o - k o,i) ^b( T b) + E A 2 ( b " a ) ~ 9 b( t ) 4 < T b) 
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(11-67) 
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^ (o) A (2)\ ^ (2) 

The factor m^k}^ - k^}), which may also be written ink} J 


1 - exp 


(-^)] 

\ kT /J 


is similar to the factor that appears in equation (11-55). This factor has the same 
significance and will be similarly represented by 


(T,p) 
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/ e aA1 


mk M 

1 - exp t "Et ) 
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( 11 - 68 ) 


which is the characteristic vibrational relaxation time associated with the upper 

A (2) ^ (2) 

states. With this newly defined quantity, the terms mE^ a _j - a) ma y 

be written — 1 ’ V . 

r B t(2) 

1,0 


Further reduction of T 3 is not possible without first introducing a special 
identity, the significance of which will become apparent only when T ls T 2 , and 
T 3 are resubstituted into equations (11-49) and (11-42). This identity, which will 
presently be introduced, evolves from an investigation of the quantity m 

x ^b-l,b ex P( E A 2 /^ t b) / Qb( t b)> which may be identified as the rate-of-energy 
transfer from the vibrational state b-1 to the state b when these states have 
Boltzmann populations with a vibrational temperature Tg. Substituting the 

A 

appropriate values for the transition probability k^_ t as determined from 
equation (II-45b), into this quantity, we obtain the identity 


n B E A 2 m Vi,b / E A 2 \ n B E A 2 [( b - a ) + ki-i,a] / E A 2 \ 

Qb ( t b) 6XP fe/ = Qb( t b) eXP \ 

(II- 6 9a) 

-(2) -(2) -(2) -(2) 

If we express k 0jl and k a _ lja in terms of k lj0 and k aja _j by means of equations 
(II-46b) and (11-48), respectively, introduce the upper-state vibrational relax- 
ation time given by equation (11-68), and transpose the left-hand term, the above 
identity becomes 
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We may add this quantity to T 3j as given by equation (11-67), without affecting 
the value of T 3 . We then observe that the terms containing exponentials may 
readily be identified with a previously defined thermodynamic quantity, as will 
be demonstrated by what follows. When the above identity, equation (II-69b), is 
introduced into equation (11-67), there results 
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With equations (II-62b) and (II-63b) we identify the coefficient 1 1 - exp 
(N-b)E A2 Tl /E A \ 

eXP (kT^ )/ Q B( T B) as E A 2 /^b( T b)- Furthermore, if we 


kT 


B 

express the quantity contained in the first bracket of the above equation in a 
form similar to that exhibited by equation (II-58b) by substitution of equation 
(II-65a), we can remove the quantity as a factor and obtain 
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This equation has nearly the desired form and would be satisfactory except 
that it would be desirable if the second factor of the first term could be related 
to a truncation factor, say £b( t b)> in order that the entire first term be 
similar to that obtained for T ls equation (11-58). Such a representation is 
possible if we note, using equation (II-44a), that (b-a) = Efr - E a and intro- 
duce the following equation (A41) from Appendix A: 


*b ( t b) = £ B ( t b) + cn-71) 

The first two terms in the second factor may be identified immediately; if we 
then multiply both the numerator and denominator of the terms that remain by 
X B (T B ) (this is a valid operation since ^b( t b) is never zero) and introduce 
equation (II-65a) to represent the denominator, we obtain the following desired 
form: 
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(11-72) 


The parameter A(T), which has been introduced here, will be one of the 
essential parameters in the model equations and is defined by 


A = E a - E A! a(5 


where 
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or 
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and 


K i,0 


(II- 73b) 



(II- 73c) 


The second representations for both A and 6 as given above follow directly 
after equations (11-44 a) and (11-47) are introduced. The quantity x(T) may be 
either equation (B14), or X F , equation (B23). The parameters A(T) and 
6( T) will be called anharmonicity parameters (1) and (2), respectively; the 
values of these parameters are discussed in Appendix B. It is worthwhile to 
point out, however, that when E Aj = and 6 = l/ the anharmonicity param- 
eter (1), A(T), is zero. Furthermore, both parameters are temperature de- 
pendent, but the dependence of (1) on temperature is only through (2). It will 

-A. 

turn out that when the effect of anharmonicity can be ignored, <5 is unity and A 
is zero, and when the effects of anharmonicity are large 0 ^ 6 < l and A(T)«E a . 


It is helpful now to recall the definition of T 3 , as given by equation (11-59), 
and rewrite equation (11-72) as follows: 


N— 2 


Tj (E v - E v+ i)g v ,v + l “ t-o j^B^ T ^B( T B) “ <5 b( T b) 


v=b 


A(T) 

^B ( T B) - 


- I ^2 \ 

” Qb( t b) b-1 > b exp ykT fi y 
n B e A 2 m r 

+ Qb( t b) kb,b_1 


(Il-74a) 


The last two terms in this equation are not readily reduced to a more simple 
representation involving known functions. One notes, however, that when the 
upper-state temperature Tg is equal to the translational or fluid temperature 


t It will be shown in Appendix B that 6 = 1 implies E A ^ = E A ^ and vice versa. 
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T, their sum is zero. As a result of vibrational nonequilibrium these terms 
are related to the net transfer of energy across the lowest state in the respec- 
tive group of energy levels under consideration (in this case the group desig- 
nated B). Comparable terms do not exist for a similar energy transfer across 
some uppermost state. In fact, the upper summation limit is implicitly embed- 
ded only in the truncation factor ^g(T-g) (recall that 3 b( T b) - ^b^b^b^b) 
and that Qg(Tg) is independent of any summation limits). These points are 
best exemplified by equation (II-74a) rewritten in such a manner that the depen- 
dence of the various quantities on the parameters are exhibited. This dependence 
has not been included in the notation heretofore in order to avoid excessive com- 
plication of the derivations and of the notation; in order to know the parameter 
dependence, one must refer to Appendix A. Rewriting equation (II-74a), using 
the appropriate functional representations, we obtain 


N-2 

E ( E v - E V4l)£v,V+l 
v=b 


-IT-V 


B 


’■( T;E a 2 4m) 


- 5b( T r 4 A ,)] Ip,[T B ;(N-b»,E A! ] 


/ 




1 


./ e a 2 k?i\ 


vv *&) 

n B e A 2 m ^b-i,b 

<3b( T B-’ E A 2 ) £B[ T B;(N-b),E A J 

■J 

Q' B [T B ;(N-b),E A2 ] 


x exp 



n B E A 2 m kb,b-l 

Q' B [ T B;(N-b),E A2 ] 


(II- 74b) 


We note that changes in (N - b) , the difference of lower and upper sum 
parameters, affects the truncation factor and the last two terms are af- 
fected also by changes in b. It should be emphasized that A(T) is independent 
of either summation limit. If we desire to approximate the effects of anharmo- 
nicity better by introducing a new level separation, say a t , above which the energy 
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levels have a different energy spacing, and also have different linear transition 

A 

probability functions, we will note that all the functions, including A, will be 
different, but only to the extent of their dependence on the parameters. The 
general functional character of equation (11-74) will be unchanged. No further 
remarks will be made concerning possible generalization of the method. 



From these remarks the correct expression to use for T 2 (see eq. (11-49) ) 
should now be evident. It is written as follows: 


b-2 


t 2 - X ( E v - E v+i) Sv+i,v 


v-a 


- - «1< T a)] • c A,( T A> 

J3 


a(t) 

( T A) £ A 2 ( T A) 


n A 2 E A 2 m . / E A 2 \ n A 2 e A 2 m . 

" Qa 2 ( t a) ka_1,a eXP \ kT A/ + Q'a 2 ( t a) ka,a_1 


(n-75) 


We recall that the vibrational temperature of these states is T A , since the 
states lie below b. Furthermore, in order to retain consistency in the notation 
we define q^ (T = q^ (T A ) and write 

%< T > - < t) 


« n - 76 > 

All the required relations are now available such that the rate equations for 
the relaxation of vibrational energy, equations (II-40a) and (11-42), may be appro- 
priately found. Substitution of equations (II-58a) and (11-75) into equation (11-49) 
yields for equation (II-40a) 
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d 11 A | 

dt[“AU(TA)] T _ v " ^Pa. (T) f A.< T A) - 3 Ai Pa)] 


“A, r.. 

+ ^b 


~°° (T) £. 


,( T A) 


M 


i - 


A(T) 

^1( T A). 


, v n A 2 E A 2 m ~ / E A 2 

+ V^a-l “ E a) Sa,a-i n ' /T . k a-i,a ex Pl up 
V QA 2 ( T A) \ kT A 


+ 


A 

n A 2 e A 2 mk a,a-i 
Q Aj ( T a) 


+ 


E b-i 


dn A~ 

dt 


T-V 


(11-77) 


Investigating the separation term (E a-1 - E a )g a a _ t , we obtain the following 
expression after substitution of the g-function (eq. (II-9)), and the number- 
density ratios (see, e.g., eq. (A17) in Appendix A): 


( E a-i “ E a) Sa,a-i - m ( E a-i ~ E a) 


k a,a-i n A 2 




I_Q'a 2 ( t a) Q Aj( T A) 


exp 


J a-i 


kT 


A > 
(11-78) 


Since the vibrational temperature of the states A t is the same as that of the 
states A 2 , namely Ta, we may relate n^ and n^ as follows: 


n 


a-i 


% 

«a,( t a) 


exp 



n A 2 

Q'a 2 ( t a) 


exp 



(11-79) 


There is a degree of arbitrariness as to whether the reduced spacing first 

occurs between the levels a-1 and a or between levels a and a+1 without 
affecting the evaluation of the series while the closed-form expressions are 
being derived for the partition functions (see Appendix A). However, the choice 
taken (eq. (11-43) ) is such that E A = E a - E a _ 1# It follows that the separation 
term may be written 

W 1 , "A^A, 1 ”- (V\ 

(E >- - Ei>) g ^-‘ ■ - «vfT) + qW a ' 1>a exp ^ 

(11-80) 
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Furthermore, one observes that this term is just the negative of the two adjacent 
terms in equation (11-77). The lower-state energy rate equation thus becomes 



a( t a)] t _ v = t^Pa^) WTa) - ^(Ta)] 


n A, 


+ r-g [ C) A 2 ^ T ^ ^A 2 ( T A) ^A 2 ( T A)] 


A(T) 

%( T A)_ 


+ E b-i 



T-V 


( 11 - 81 ) 


The reasons for introducing the term leading to the identity given by equation 
(II-69b) should now be partially evident. It leads to a relatively simple repre- 
sentation for T 3 and, furthermore, to the cancellation of the separation term 
in equation (11-49) (the term in the series involving A-state levels that separates 
the Aj- and A 2 -groupings of energy levels). It will be shown next that the ex- 
pression introduced will also lead to a term in the upper-state vibrational- 
energy relaxation equation that is rate limiting in its effect on dissociation. 

Before making this identification we substitute the summation represented 
by T 3 , equation (11-72), into equation (11-42) to obtain the upper-state energy 
relaxation equation 


Jt[ n B^B( T B)] T _ v ~ T B pB( T ) £ B( T B) " 3 b ( T B)] 


1 


A(T) ' 
^B ( T B) 


n B Ea 2 m A /E Az \ n B m 

Qb( T b) b_1 > b 6Xp \kT B / Q b (T b ) 


k b,b-i 


(II- 8 2 a) 


The significance of the last two terms in this equation becomes apparent if we 
substitute the appropriate expression relating the quantities and Q;b(T b ), 

using equation (II-38c), and introduce equations (11-41) into the resulting 
equation: 
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&[»b«b<Tb)] t _ v - ^[q|(T)£ B (T B ) - QB(Tb)] [l - 

+ ( E b-i - E b) gb,b-i - E b-i L"drj T V 


n B e A 2 m 
( T b) 


*b-*.b e*p(j^) 


+ 


n B e A 2 m 

q b ( t b) 


k b,b-l 


(II-82b) 


In tins case, the quantity in the above equation involving the g- function is given 
by 


- E 


A 2 


k b,b-i n B 

Q B ( t b) 


n A k b-i,b 

Qa(t a ) 



These terms do not lead to a complete cancellation of their equivalents in the 
above rate equation in the manner that the analogous terms lead to cancellations 
in equations (11-77), but lead to the expression displayed as the last term in the 
equation that follows : 


^[ n BqB( T B) ^ Tr 
ut L J T-V 


[^b (t) £© (t b ) - q B (T B ) 

B L 


1 


±(T) ' 

^B ( T B) 


p [ dn A~ 
b_1 - dt _ T-V 


(l) 


E A 2 n A mk F 


( t a> t ) 


X 


1 


^B Qa ( T A) 
n A Qb( t b) 




(II-82c) 


-( 1 ) 

where we have introduced the rate constant k F (T^,T), defined by equation 

(II-39b), to simplify the notation. The dependent variables n^ and n^ 

are not necessarily related in terms of a ratio of partition functions in the 

manner that we found for n A and n A . Their values depend on the complete 

A i a 2 

solution of the model equations. Only when the upper and lower temperatures, 
T^ and Tg, are equal will the ratio of upper and lower population densities be 
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related to a ratio of partition functions. The last term in the above equation will 
then be zero. It will be shown in Chapter IV, when the parameter b is appro- 
priately evaluated, that the last term has the effect of limiting the rate of disso- 
ciation; for this reason the level b will be denoted the "rate-limiting" level. 

At this time the term represents the net rate of transfer of energy across a 
level b that separates those levels with populations with a temperature 
from those levels with a temperature Tg. 



The equations considered thus far are applicable only for translation- 
vibration energy interchange in a gas containing a relatively dilute mixture of 
molecules, that is, for vibrational nonequilibrium processes that may be prop- 
erly described by the first chemical equation, equation (II— 1) . The equations 
that finally result (eqs. (II-34a,b), (11-81), (II-82c), and eqs. (11-39) ) are suffi- 
cient to describe the internal properties of a gas, provided the effects of disso- 
ciation are negligible (and also the effects of resonance vibration- vibration 
energy interchange, which will not be considered in this paper). When there is 
appreciable dissociation and/or recombination, additional complication will 
ensue. In the next section these complications will be taken into account. 


(II-B-3). Dissociation 

This section deals with obtaining the set of rate equations that are valid for 
describing the translation- vibration energy interchange involving dissociation. 
Four rate equations, subscripted D, will result from the analysis, and will be 
combined with the equations found in the previous sections to yield the relations 
needed for the study of vibration-dissociation coupling. 

The procedure for obtaining the model equations for the dissociation pro- 
cess is similar to that already used. In the previous sections the effect of 
vibration-translation energy interchange, resulting in vibrational relaxation, 
was investigated by obtaining zeroth and first moment equations of the popula- 
tion distribution using equations (II-5) as the basic equations. These moment 
relations were then simplified by the introduction of various assumptions yield- 
ing the final equations (11-34), (11-39), (11-81), and (II-82c), which contain a 
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relatively small number of parameters. The procedure involved a method of 
accounting for molecular anharmonicity effects both in the energy level spacing 
and in the transition probabilities. 

This section will be based on equation (II-7). The procedure for simpli- 
fying this equation will be similar to that used in the previous sections in that 
the zeroth and first moment equations of the population distribution will again 
be obtained. The procedure will deviate to the extent that the preferential dis- 
sociation model of Marrone and Treanor (58) will be introduced to obtain simple 
relations describing the dissociation. Anharmonicity effects will be accounted 
for by using the energy-level spacings introduced in the previous sections. The 
preferential dissociation model of Marrone and Treanor is described in detail 
in Appendix C. Suffice it to say here that the Marrone and Treanor model pro- 
vides a method whereby one may reasonably approximate the level transition 
probabilities kv,N anc * k N,v- Their scheme is simple and its employment intro- 
duces few additional unknowns into the problem. The essence of the model is 
that the probability of dissociation is exponentially weighted in such a manner 
that dissociation from the uppermost levels is most likely. In the derivations 
that follow the Mar rone -Treanor model will be slightly modified. In particular, 
it will be assumed that dissociation and recombination involving molecules in 
the A-group of energy states are so improbable that they may be neglected. 

The zeroth moment expression of the population distribution for the inter- 
actions involving dissociation is obtained by summing equation (II-7) over the 
subscript v, including only those states in the sum appropriate to the group of 
energy states under consideration. The result for the A-states is 

dn Al r dn vl V 1 r ^r -r 1 

= L hjH ” 2j 2j [~ k v,N n r m r + k N,v n a m r (11-83) 

J u v=o L J u v=o r=i J 

As already mentioned, however, dissociation and recombination of molecules in 
these states is assumed to be extremely improbable. Therefore 




0 


v < b - 1 


(n-84) 
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and we have 


Ldt J D 




Similarly, the first moment relation of the A-state population distribution, 
obtained by multiplying equation (II-7) by E v and then summing, is also zero; 
that is. 


d [ n A^A( T B)] D 


( 11 - 86 ) 


The direct effect of dissociation will then only be felt by the upper group of 
molecular energy states. 

Before proceeding with the derivations it will be helpful to remove from the 
equations the effect of the dependence on r, that is, on the kind of species in- 
volved in the collision. This may be done by either temporarily ignoring the 
effect, as was done earlier, or by removing it by an appropriate averaging 
procedure. We do the latter by writing 


Ldt J D 


= - n v mk V)N + n|mk N v 


(11-87) 


where the rates ^ and v are defined by the averages 


^v,N = E 


nr t* y. 

In" k v,N 


(II-88a) 


“ Li m k N,v 


(II-88b) 


As previously remarked, only dilute molecular mixtures are being considered; 
that is, we assume that the mixture contains relatively small concentrations of 
molecules and of atoms obtained by molecular dissociation, plus some inert 
species that is the major constituent (e.g., argon). If the population fractions 
m r / m are sufficiently small for the molecules and reacting atoms compared 
with that for the inert species , equation (11-87) may be accurately represented 
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if the quantities k are those of the inert species only. To simplify notation we 

* 

will henceforth ignore the bar on k . 


Using the equations given above we now write the zeroth and first moment 
expressions for the B-group of molecular energy states as follows: 

N-i ^ N-l ^ 

= - E ^v,N n v m + E k N,v n l m (n-89) 

v=b ’ v=b 



d [ n B q B( T B)L N-l a N-l^ 

dt = - L k v,N E v n v m + k N)V E v n|m (11-90) 

v=b v=b 

where n^ and are the quantities defined previously by equations (11-21) 

A .X 

and (II-38b). The transition rates and k N ^ y peculiar to the preferential- 

dissociation model, are derived in Appendix C (see eq. (C14)). The expres- 
sions given in the appendix may be substituted into the above equations. In 

N-l N-l 

addition, if one notes the identities Z P v = Z E v p v = ^b( T f)> ancJ 

v-o v=b 

N-l 

Z E vPveq = ^B (“ T|j), we °^ a i n the following equations: 
v=b 

|"dn B 1 * k-ri (T) 

[— J D = - kFeq.( T ) V ( T B.T)»Bn> + K *( T ) "l"> d 1 -^) 


dn B q B ( T b) 


dt 


J D 


~ “ q B( T F) n B m ^Feq/ T ^ V ( T B’ T ) + q B^ T u) 


k Feq.( T ) n l m 

V T > 

ai-92) 


where the quantity ^Yeq ( T ) equilibrium dissociation rate (i.e., the rate 

applicable when the vibrational mode is assumed to be in equilibrium). This 
quantity is discussed further in Chapter III. The function Kg( T) is the partial 
equilibrium constant (see Appendix C, eq. (C13)); the quantities p v are norma- 
lized probabilities that are also discussed in Appendix C; and the factor V(T B ,T), 
called the vibrational coupling factor after Marrone and Treanor (58), is a known 



function containing ratios of upper-state vibrational partition functions as 
follows (see Appendix C for discussion): 


v( B > " Qb( t b) ( Sb(- t u> 


(n-93) 


The variable Tp is the characteristic temperature at which energy is removed 
by dissociation (distinct from the dissociation temperature, 0 jq; see Appendix A) 
and is defined as follows (from eq. (C7)) 


1 





(11-94) 


The negative quantity T^j (a constant) in the above relations may be considered 
as the characteristic "vibrational temperature" (pointed out by Marrone and 
Treanor (58)) at which energy is returned to the vibrational mode after recom- 
bination. The value of T-q also describes how rapidly the dissociation prob- 
ability drops off for the lower vibrational energy levels. The assigned value of 
T-y is discussed in the next chapter. 


(II-B-4). Combined Effect of Vibrational Excitation and Dissociation 


The model equations obtained in the earlier sections, equations (11-34), 
(11-81), and (II-82c), may now be generalized to be valid when the effects of 
dissociation and recombination are important. The generalization for the num- 
ber densities requires the superposing of equations (11-34) , (11-85), and (11-91) 


(n-95a) 
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(II- 8) and 
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dt 

= m |n B : 
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dt 
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- (T)nlm 

“IT " k Feq.(T)V(T B ,T)n B m + 


(II- 95 b) 


The energy rate equations are obtained by superposing equations (11-81), 
(II-82c), and (11-86), (11-92), respectively. We obtain 
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(II- 95c) 
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b- 1_ dT 


+ E A 2 n A m ^F\ T A ,T ) 
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^B( T F) n B m ^Feq. V ( T B ,T ) 


^B ( T U) n a m ** Feq. ( T ) 


(II-95d) 


In equation (II-95c) the densities n Ai and n A ^ have been expressed in terms of 
the dependent variable n A and appropriate ratios of partition functions (see, 
e.g., eq. (A19) in Appendix A). 

The above equations are cumbersome in their present form, and a more 
concise formulation is desirable. With this in mind the following notation will 
be adopted. 


^(T,T A ) = q^(T) - q Aj (T A ) 

(II- 96 a) 

6I 2 ( T > T A) = %( T ) - %( T A) 

(II- 96b) 

^b( T > T b) = q£( T ) ‘ «<*b) 

(II- 96c) 
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J 


%( T * T a) ~ %( T ’ T a) ^Aj ( t a) 
%( t ’ t a) = ^(t,t a) *A 2 0 a) 

%( t > t b) = ^( t > t b) £ b( t b) 


(II- 97a) 
(II- 97b) 
(II- 97c) 




There is some redundancy in the above notation, equations (II- 96), since 
q~ = qg and hence fij^ = fig when T A = Tg. However, this redundancy is 
removed in equations (II- 97). The quantities fi^ and fig are not necessarily 
equal even when the temperatures T A and Tg are equal. The quantities fi 
defined above may be considered as measures of the departure of the respective 
grouping of vibrational energy levels from local thermodynamic equilibrium. 
When the argument temperatures are equal, these quantities are zero, and as 
the differences of the argument temperatures increase in a uniform manner the 
quantities themselves also increase in a uniform manner. In addition to the 
foregoing we also write 



= k Feq> (T) V(T B ,T)n B m 



4 m 

k b ( t ) 


(H-98a) 
(II- 98b) 


This notation has been used previously by Marrone and Treanor (58,93). The 
first expression denotes the rate at which molecules dissociate and the second 
is the rate at which molecules are formed as a result of recombination. That 
these expressions apply appropriately to the molecules and not just to the upper 
grouping of energy states is evident if we sum equations (II-95a) and (II-95b) 
and note that n = n A + ng is the total molecular density including all vibrational 
states. 

We now rewrite equations (11-95) using the notation of equations (11-96), 
(11-97), and (II- 98). There results 

TT * m ["B i R > ( T B’ T ) - "A i £ ) ( T A' T )] 
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(II- 99b) 
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(II- 9 9c) 
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(II-99d) 


The above equations together with the atom conservation equation 

n a n a 

n + ~ =n A + n B +- 2 - = constant 


(II-99e) 


constitute the complete set of equations required to determine the internal 
properties (vibrational mode) of the gas and thereby the effects of vibration- 
dissociation coupling. Since these equations are a complete set it is perhaps 
fitting to summarize the physical significance of the various terms and then 
compare these equations with those used by the other researchers. 

The individual terms in the first two equations are easily identified. The 
terms in equation (II-99a) represent the gain and loss of the A-state molecules 
due to the transitions yielding B-states. The first term in equation (II- 99b) 
represents the converse processes, that is, the net effects on the B-states due 
to the transitions that yield A-states. The second and third terms represent 
the loss and gain of B- state molecules as a result, respectively, of molecular 
dissociation and recombination. 



The energy equations (II- 99c) are slightly more complex than the popula- 
tion equations just described. The terms within the braces of the A-state 
vibrational- energy equation (II-99c) represent the relaxation of vibrational 
energy for the subgroups A t and A 2 of the A- states where each sub-group has 
different associated characteristic relaxation times. The first term in the 
braces is for the subgroup of states denoted that has an associated charac- 
teristic relaxation time and the second term is for the subgroup denoted 
A 2 that has a relaxation time It is worthwhile to remark that although 

t A 2 < TAj* t ^ lis e ^ ect is counteracted by the partition- function fractions 
Q Ai /Q a an d Qa/ Q A- a resu ^» the lower state vibrational relaxation de- 
pends most strongly on the first term, in particular, when the temperature T^ 
is not large (e.g., see Chapter IV). The third term in equation (II-99c) repre- 
sents the rate of gain (or loss) of vibrational energy from the A-states due to 
transitions from B to A states (or the reverse). The first term in equation 
(II-99d) describes the vibrational-energy relaxation occurring in the B-states. 
The second term is the converse of the similar term in equation (II- 99c), while 
the third and fourth terms represent the loss and gain of the vibrational energy 
of the B-states due to dissociation and recombination. The last term represents 
an additional gain (or loss) of net vibrational energy due to the interaction with 
the translational mode. Such an effect arises when molecules make transitions 
from a system at one temperature T^ to a system at a different temperature 
T-g. This term is zero when the two temperatures are equal. It is difficult to 
ascribe a direct physical significance to this term, because of the manner in 
which it was introduced. It will be shown in Chapter IV that a similar term 
yields a rate-limiting effect on the rate equations. In this respect the term may 
produce the effect analogous to the "bottle-neck” discussed by Bray (14) and by 
Pritchard (71). 



It is worthwhile to compare these equations (11-99) with those used by 
Marrone and Treanor (58) in their study of vibration-dissociation coupling. In 
the Marrone-Treanor model the assumption was made that the distribution of 
energy in the vibrational energy states can be described by a single vibrational 
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temperature Ty. Here we have introduced an additional temperature into the 
description. Furthermore, in the Marrone-Treanor model the effects of 
anharmonicity were accounted for only to the extent that more accurate quan- 
tities for the thermodynamic energy function were used (it will be shown that 
this is not essential). It was not possible with their formulation to account for 
effective decreases in characteristic vibrational relaxation time for the upper 
levels that occur as a result of molecular anharmonicity. 

Although the model presented here contains additional detail as regards the 
description of the vibrational degrees of freedom, it also includes additional 
complication. There may be problems where the additional complication is not 
warranted and the Marrone-Treanor model will be of value. With this in mind 
it is worthwhile to investigate the functional form of the rate equations derived 
in this paper as the complexity is reduced by further assumptions. 

We first assume (58,93) that a single temperature is appropriate for a 
description of the vibrational mode, that is, Ty = T^ = T-g. We observe that 
the last term in equation (II-99d) then vanishes. The subscripts A and B are no 
longer significant and we, therefore, can combine equations (II- 99a) and (II- 99b) 
to obtain 

[*]-[£!,♦[£]» 

Likewise we combine equations (II- 99c) and (II-99d) to obtain 
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(II- 101) 



where we have used 


nq(T v ) = n A^A^ T v) + u B^b( T v) (11-102) 


A more useful formulation of the above equation is obtained by first relating the 
densities n^ and n B in terms of the net molecular density n and appropriate 
ratios of partition functions. There results 
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We then combine the second pair of terms in braces to obtain the desired 
formulation 


dn q( T v) n ( Q Ai(Ty) a Qa 2 -b( T v) a 

dt “ Q v ( T v ) | t Ai fi A 1 ( T * T v) + t b ^A 2 -B ( t > t v) 
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3 B (r v)[S] p + 5 b<- T U>[f] E (II-103b) 


where the following relations have been introduced, consistent with the previ- 
ously defined similar quantities: 

Q A 2 -b( T v) = Qa 2 ( T v) + Qb( T v) (H-104a) 

Qa 2 -b( T v) £ A2-b( t v) = Qa 2 ( t v) £ A 2 ( t V) + Qb( t v) ^B^v) (II- 104b) 

%-b( T v) = q A 2 -B( T v) £ A 2 -b( T v) (H-104c) 

The subscript A 2 -B is used to denote the union of the vibrational states 
contained in the groups A 2 and B. 
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Equations (11-100) and (II-103b) are analogous to those used by Marrone and 
Treanor (58). We observe, however, that to account for the effect of anharmo- 
nicity it is not sufficient merely to use more accurate values for the vibrational 
energy function q in the same manner as Marrone and Treanor, but an addi- 
tional term must be introduced. This term is the second term given in the 
braces of equation (II-103b). If we are able to neglect the effects of anharmo- 
nicity (that is, if A(T)/q^_ B (Ty) ~ 0 and E^/E^ ~ l) the terms in braces 
maybe combined (see, e.g., eq. (A19) in Appendix A). Equation (E-103b) may 
then be written as follows: 

^ S(T,T V ) - q B (T F ,[f] F + q B (- Tu )[^] r 01-105, 

where 

fi(T,T v ) = [q°°(T) - q°°(T v )] *(Ty) (II- 106a) 

£(T V ) 

= q(T)-y^y - q(T y ) (II-106b) 

The exclusion of subscripts on the relevant quantities denotes that the vibra- 
tional states are all considered together as a group. The representation given 
above, with the exception of the dissociation and recombination terms, is simi- 
lar to the form of equation (II-58c), which prompted the previous remarks 
(given after eq. (II-58c) ). In this equation the effects of anharmonicity are 
assumed negligible (or ineffectual), and it is only the existence of the truncation 

factor ratio X(Ty)/ ICE) that causes the equation to differ from that used by 
Marrone and Treanor (58). The differences in the equations are not important 
provided that the temperatures T and T y are both much less than the charac- 
teristic dissociation temperature 0j) = E n -i/k. The truncation-factor ratio is 
then of order unity. For this case, however, there is negligible dissociation 
and the problem is N not of great interest. At the higher temperatures, say 
T = 16,500° and Ty = 300° K, the factor £(Ty)/ £(T) for the example case of 
molecular oxygen is about 3, and hence the differences become important and 
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the Marrone-Treanor formulation questionable. This is discussed further in 
Appendix D. 




n~C. The Rate Equations Applicable to Fluid-Flow Problems 

In the previous sections the model rate equations have been derived that 
are applicable for studies of vibration-dissociation coupling. In these equations 
the time derivatives have been introduced without any discussion regarding 
their significance. To apply such equations to problems involving a moving 
fluid requires that some consideration be given to these derivatives. Formally, 
the equations that have been derived apply to a constant-mass, constant- volume 
system such as provided by a closed rigid container. To be applicable to a 
moving fluid this constraint must be removed. The procedure for doing this is 
described clearly by Vincenti and Kruger (p. 248 (97) ). Briefly, we let a rep- 
resent any one of the variables (per unit volume) n A , n B , (n A q A ), or (n B qg). 
The equations applicable to fluid flows, if we neglect diffusion effects, are given 


D2 1 da 
Dt “ P dt 


(11-107) 


where 2 is the intensive (per unit mass) counterpart of the variable cr (<r = p2), 
and ^ is the substantial or material derivative The variable 

? is the mass velocity of the fluid. Since the equations for the derivative da/dt 
have been derived for a closed fixed-volume system, the density p may be 
taken inside or outside this derivative as we wish. 

The equations derived previously, equations (11-99), are readily modified to 
a form that is applicable for flow field calculations by simply changing the de- 
pendent variables to the 2 quantities. The equations will be listed below in 
terms of the appropriate variables, and this list will serve as a summary of the 
results heretofore obtained. The intensive variables and corresponding 
to the 2 quantities are defined by 

n m 

r m P N 0 

m = A, B, or a (n- 108 a) 
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r 

e i = TiQi = YiMi 


i = A, A u A 2 , or B 

(II- 108b) 

where N 0 is Avogadro ! s number. The rate parameters in the equations are 
defined according to the relations 

4 eq .(T) = N 0 k| eq> (T) 

(II- 10 9a) 

k R (l) ( T B’ T ) = N okR’ (l) (T B ,T) 

0I-109b) 

4> (1) (T A> T ) = N 0 k^W(T A ,T) 

(II- 109c) 


With these quantities the rate equations become 1 ’ 
DV A 


Dt p ^m 


r B Xf4V)-r A ^ 4 ' ( Vb.t) 

_ r T m r 

dy b dy a _ /py\ + /dt\ 
yDt J-p \Dt 


Dt 


IT ‘ ^ | Q A l ( T A)»A 1 (T.T A ) * Q A! (T A )(^j V T,T A ) 

A r (T) 


Dt 


(n-110a) 


(II- 11 Ob) 


1 - 


%( t a). 


dy a 

+ R o e b-i dt 


(II-llOc) 


tin equations (11-110) the effects of the type of collision partner have been re- 
introduced into the equations. This requires introducing the superscripts and 
sums for r. These equations are the model equations for the reactions (II— 1) 
and (II- 3) when the effects of vibration-vibration energy interchange (II— 2) may 
be considered negligible. 
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Dt y B n B( T ’ T B) E 


1 - 


A r (T) 

q b( t b) 


- R 0 ^b-l 


Dt 




- OB( T F)(^) r + M- T U)(^) R + ^A/WmE^ 

- - — r 

'b t a ; 


v (i) , I. r B QA( T A) 


- 0 b-n Tl 




(II-110d) 


(II-110e) 


The last equation is obtained by differentiating the atom-conservation equation 
(II-99e). The molecular rate terms subscripted F and R are defined by 


(tj) f ' PV<T B .T)r B T m E £ kFeq.(T) 

(T) 


Dy\ _ P*y a y m 
Dt / R 


2 -v 2 -v V 

y — k 

/c b (t) y r m 


r 

Feq. 


(II- 11 la) 
(II-lllb) 


The constants 0^,-1 and © A are characteristic temperatures defined, respec- 
tively, by 0 b _j = Ej^/k and = E^/k. The quantity R 0 = kN 0 is the 
universal gas constant. The symbols containing a caret (-) differ from the like 
quantities without the caret by the factor N 0 , Avogadro's number. 

An alternate form of these equations that turns out to be more useful for 
numerical work (see, e.g., Appendix F) is obtained by re-expressing them in 
terms of the dependent variables J A , T B , T A , Tg, and 7 a . The bar here de- 
notes that these quantities are natural logarithms of their unbarred equivalents, 
that is , 


y A = iny A 

(II-112a) 

y B = 

(II-112b) 

y a = *ny a 

(II-112c) 
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The rate equations expressed in terms of these quantities are 


Dt 


DVr 


= py m [exp ( r B - y A ) Z^ ^^(t b »t) - Z ^ kj’ (l) (t a ,t)] 


dT = " exp ^ A - y B ) 


Dt 


-(-B)[(f) F - <g) J 


(II-113a) 

(n-113b) 


DT, 


m C v a ( T a) 


5 v^I Q a ‘ (Ta)Sa ‘ <t ' Ta) 


Ai 


+ Q A ( t a) Tr °A ( t ’ t a) 


B 


1 - 


A r (T) 


% ( t a) 


+ [Ro9b-, - 1 a( t a>] 


£16 

Dt 


(II- 11 3c) 


DT 


B 


01 Cy B (T B ) 


Z 77 ^B ( t » t b) 

r T B 


' A r ( T ) 

. "MV 


by\ 
Dt ) F 


- DT A i 

q B( T B) eXp ( V A " y B)^T “ ex P(- y B)j[qB( T F) " q B( T B)] 
f- q B^ _ T U) " q b( T B^ (dt) r | + H ° G A 2 p y m I 1 " exp [j 0 A 2 ^T " T^) 


* Z TT k R (l) ( T B.T) 


r y m 


D^ 

Dt 


= - 2 exp (- y a ) 


[-(©,*©,] 


(II-H3d) 


(H-113e) 


t 


: Thi s equation may be obtained from equation (II-82b) or by appropriate 
arithmetic combination of equations (11-110). 
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where the specific heats C v (T^) and C v (Tg), defined by the derivatives 

A Jj 

dq^(T A )/dT A and dq^T^/dT-g, respectively, are discussed in Appendix A. 



The following table lists the quantities contained either explicitly or 
implicitly in the rate equations (11-110) or (11-113); alongside are the equation 
numbers of the formulas that may be used to evaluate the indicated quantities. 


Table II- 1 

Table of Formula Numbers of Quantities Contained in the Rate Equations 


Quantity 

Description 

Equation No. 
(Quantity) 

Equation No. 
(Quantity Primed) 

k^(T B ,T) 

Reverse transition rate, 
lower levels 

(II- 39 a) 

- 

k?(T A ,T) 

Forward transition rate, 
lower levels 

(II- 39b) 

- 

*Feq.< T > 

Forward dissociation rate, 
upper levels 

(III- 13) 

- 

*b< t > 

Partial equilibrium constant 

(C13) 

- 

K( T) 

Equilibrium constant 

(A8) 

- 


Partial vibrational partition 
function 

(A14a) 


q a 2 


(A14b) 

(A15) 

Qa 


(A14a)+(A14b) 

- 

Qb 


(A14c) 

(A16) 

Q v 


(A18) 

- 

A°° 

fi Aj 

Measure of departure from 
equilibrium of nontruncated 
harmonic oscillator 

(11-96 a) 

- 
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Table II- 1 (Continued) 


Quantity 

Description 

Equation No. 
(Quantity) 

Equation No. 
(Quantity Primed) 

C *A, 

Vibrational specific heat 

(A44a) 


S 


(A44b) 

- 

°v A 


(A47a) 

- 

C v B 


(A44c) 

- 

c v v 


(A47b) 

- 

V(T b ,T) 

Coupling factor 

(n-93),(C-10) 

- 

t f 

Characteristic temperature 

(11-94), (C7) 

- 

A(T) 

Anharmonicity parameter 1 

(II-73a),(B26) 

- 

6(T) 

Anharmonicity parameter 2 

(II- 73b) 

- 

X R (T) 

Slope function 

(B28a) 

- 

X f (T) 

Slope function 

(B28b) 

- 

T Ai 

Vibrational relaxation time 
lower states 

(B25a) or (in-lla) 

- 

t b 

Vibrational relaxation time 
upper states 

(B25b) 

- 

( t a,/ t b) 

Ratio 

(B25c) 

- 


The rate equations (11-110) or (11-113) are still sufficiently general to be 
applicable to flow fields that are spatially multidimensional and unsteady. Such 
generality, however, will not be considered in the chapters that follow where 
only steady one-dimensional flows are investigated. A more convenient notation 
for the rate equations is given by 


(n-114a) 


_ d2 (n-114b) 

dx 

where 2 denotes either y A , y B , y a , e A , £ b» T A’ or T B’ The Q 2 quantities 

represent the differential change of the associated quantities along streamlines 
and are obtained directly by dividing the right hand side of either equations 
(11-110) or (n-113) by the mass velocity v. 

The rate equations also contain a number of basic parameters that are for 
the most part dependent on the molecular species chosen to be investigated, in 
particular, the quantities k^p , k^p , A(T), E^, E A , N, a, b, k Feq (T), and Tjj. 
Once appropriate values are assigned to these quantities, one can investigate 
the effects of vibration-dissociation coupling on fluid flows. The next chapters 
will be directed primarily toward a study of these quantities. 
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CHAPTER ffl 


SPECIAL SOLUTIONS OF THE MODEL EQUATIONS 



A set of rate equations applicable for studies of vibration-dissociation cou- 
pling was derived in the previous chapter. There remains, however, the problem 
of assigning values for the embedded parameters. As an additional complication, 
many of the parameters are not constants, but are functions of temperature. 
Unfortunately, there is no simple procedure available for this problem. One 
must specify values for all of the quantities before a solution can be obtained, 
and only after one has compared the solutions with observed (experimental) 
data can the parameters be confidently evaluated. It turns out that the values 
of the parameters have to be tentatively assigned and the results studied. These 
values are altered and the effect of the change is then investigated. In this man- 
ner one may hope to obtain a final set that effects good agreement when the 
solutions are compared with experimental data. The problem is further com- 
plicated by the fact that the complete solutions of the rate equations are rela- 
tively costly both from the point of view of expense and of computer time, since 
the solutions require an integration of the model equations and in addition the 
equations of flow. These equations are nonlinear and require electronic machine 
integration. Fortunately, there are a set of special solutions that do not require 
numerical integration and, therefore, are relatively easy to obtain. In particular , 
in the quasi-steady region behind a normal shock wave, to be described later, 
certain combinations of the rate equations may be set equal to zero (it will be 
shown that in this region the quantities t a> t b> and r B / are constant) and 
the resulting equations solved. An investigation of these solutions yields con- 
siderable insight into the effects that occur as a result of varying the values of 
embedded parameters and thus are of great value. This chapter concerns an 
investigation of the effect of changes of the embedded parameters on the quasi- 
steady solutions. 
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HI -A. Introductory Comments 


It is well known (see, e.g., 11, 12, 13, 48, 85, 90, 94) that the region behind 
a normal shock wave, subject to vibration-dissociation coupling, can be sepa- 
rated into three rather distinct zones, one of which will be discussed in this 
chapter. The justification for this separation will be made clear (Chapter IV) 
when the complete numerical solutions are exhibited showing the full relaxa- 
tion to equilibrium behind a normal shock wave. For the present, the zones 
will merely be described with an eye to lending purpose to the solutions to be 
obtained. In the order of their occurrence behind the shock wave the zones are: 
(1) the transient zone, immediately behind the shock, in which vibrational 
relaxation occurs and where chemical effects are negligible; (2) the quasi- 
steady zone where both the changes in the fractional molecular population, 
n v /n, of the vibrational energy states, and the effect of atom recombination 
may be considered negligible, and where dissociation is observed; and (3) the 
final relaxation zone in which the nonequilibrium vibrational and chemical 
effects diminish rapidly until the gas reaches equilibrium. 

The transient zone (1) and the final relaxation zone (3) persist for rela- 
tively short periods of time, so that experimental study of these regions is 
difficult. Furthermore, an analytical study of these zones requires an integra- 
tion of the complete set of master equations (or of the equivalent set of equa- 
tions introduced in this paper), and such investigation is not simple. As a 
result the rate constants and molecular kinetics are least understood in these 
regions. The quasi-steady zone, however, persists for a relatively long period 
and has been the subject of considerable experimental study (19, 22, 56, 86, 99). 
An analytical investigation of the master equations for the flow region charac- 
terized as quasi-steady requires solving the system of equations obtained by 
setting the derivatives d(n v /n)/dt equal to zero. This, compared with the 
integration of the complete set of rate equations, is a much simpler problem 
and thus also a subject of considerable study (see, e.g., 13, 48, 73, 85, 94). 

For these reasons comparisons between theory and experiment are most 
easily made in the quasi-steady zone, the temperature dependence of the 
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dissociation rate being that function used as a basis for comparison. Treanor 
and Marrone assessed the validity of their nonpreferential (90, 91) and prefer- 
ential (58) vibration-dissociation models by making experimental comparisons of 
the pre-exponential temperature dependence of the effective dissociation rate 
constant. Similar studies were also conducted by Keck and Carrier (48) and by 
Snider (85). Before pursuing the subject further it is worthwhile to digress 
briefly to obtain somewhat of an understanding of the temperature dependent 
qualities of the dissociation rate constant. 



An intuitive feeling for the significance of the temperature dependence of 
the dissociation rate can be obtained from "Collision Theory" (see, e.g., Vincenti 
and Kruger (97), p. 221), in which an expression is obtained for the dissociation 
rate. Basically, it is assumed that only those collisions are effective in pro- 
moting dissociation for which the net energy of collision (i.e., the relative line- 
of-center energy plus a contribution of internal energy (vibrational and/or 
rotational) ) is greater than the molecular binding energy. This is a quasi- 
equilibrium theory in that the vibrational mode is assumed to be equilibrated 
along with the rotational and translational modes. It is shown that the tempera- 
ture dependence of the dissociation rate constant is given by T 3 / 2 “ s exp(-e Q /kT) 
where the constant e Q is the dissociation (or activation) energy and the quantity s is 
one-half the number of square terms that contribute effectively to the energy of col- 
lision (see, e.g., (97), p. 227). The pre -exponential temperature dependence is 
given by the factor t 3 / 2_s . For example, if only the relative transitional 
motion is effective in the dissociation process, then two quadratic terms are 
important (i.e., the two coordinates required to specify the line of centers) and 
s = 1. The dissociation rate is then proportional to T 1 / 2 exp(-e 0 /kT); the pre- 
exponential factor is, of course, T 1 / 2 . If additional internal energy modes 
must be considered then s > 1. It may occur that only a fraction of the energy 
of some modes contribute to dissociation in which case fractional values of s 
may be observed. It should be emphasized that Collision Theory is a quasi- 
equilibrium theory and hence of value only in that it provides some insight to 
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T 

the effect of the internal modes on the temperature dependence of the dissocia- 
tion rate constant. If, as an upper bound, we presume that the translational, 
rotational, and vibrational modes are all effective during an atom-molecule 
collision, we might expect on the basis of this theory a rate constant with a 
temperature dependence of T“ 3 / 2 exp(-e 0 /kT) (see, e.g., (97), pp. 223 and 227). 

In the case of argon -oxygen collisions Camac and Vaughn (19) obtained 
experimentally a dissociation rate constant with a pre -exponential factor of 
T -1± * 2 over the temperature range 3,400° < T ^ 7,500° K (in this case 
s = 5/2±.2). Wray (99), in an independent experiment over a wider temperature 
range, 5,000° < T ^ 18,000° K, suggested a dissociation rate without any pre- 
exponential temperature dependence (inferring s = 3/2). It will be shown later 
that there is considerable scatter in the experimental data and the problem of 
accurately determining a pre -exponential temperature dependence is difficult. 

It will be shown, however, that a T“ 1//2 factor is not unreasonable for either 
Camac and Vaughn r s or Wray r s data. 

The non-preferential dissociation model of Treanor and Marrone (90) is 
consistent with a T" 1 factor (91), but fails to predict a satisfactory induc- 
tion time (58), that is, delay time before the onset of molecular dissociation 
(this will be discussed further in Chapter IV). The preferential dissociation 
model of Marrone and Treanor (58), although providing qualitatively correct 
transient characteristics, yields a T~ 2 pre -exponential variation (58, 90). 

The exponent is too large on the basis of collision theory or experiment. It 
was partly the inadequacy of either of the two models to describe quantitatively 
both the transient and quasi -steady behavior that this thesis study was suggested 
by Doctor Charles E. Treanor of the Cornell Aeronautical Laboratory. 

m-B. Quasi -Steady State in a Dissociation Experiment 

In the work that follows the effects associated with the model will be 
investigated in a manner similar to that of Treanor and Marrone (58, 90, 91, 

94) and of other researchers (13, 48, 73, 85). First the conditions that char- 
acterize the quasi-steady zone will be pointed out, and from these conditions 
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the required equations will become evident. Since only systems that contain a 
dilute mixture of molecules are being considered, the translational temperature 
may be assumed constant behind the shock. A brief explanation will be included 
on the procedure used to solve the system of equations for the quasi-steady 
zone. A discussion will then be given of the values of the parameters that will 
be used, and finally the numerical solutions will be presented and discussed. 



In the introductory comments it was pointed out that the quasi-steady zone 
is that region of the flow field where the ratio (n v /n)* (the asterisk denotes 
values in the quasi-steady zone) is not a function of the independent variables 
of the problem (time or distance measured from the shock). Furthermore 
the effects of molecular recombination may be neglected, that is, 



The first constraint yields the relations 



(- 




I^A (Ta) + h?h B <Tb> 



(N 0 E V ) 


(m-la) 

(in-ib) 


where 


= constant 


(III-lc) 


" 7 A + y B 


(IH-ld) 


These relations are satisfied by the conditions^ 

T^* = constant (III -2 a) 

T b * = constant (III-2b) 

tin the next chapter it also will be shown that equations (HI-2) are indeed 
valid conditions for the quasi-steady zone. 
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The equations required for the investigation of the quasi-steady zone are 
readily obtained by substituting equations (in -2a) and (in -2b) into equations 
(II-113c) and (II-113d), respectively, and substituting equations (m-2c) into 
the equation obtained after subtracting (H-113a) from (II-113b). Neglecting 
the atom recombination terms , we thus obtain the following set of transcen- 
dental equations to be solved for T A , Tg, and (Vg/Y A )*: 


< t a> % <t,t£) + Qa 2 (T^) S2a 2 (T,Ta) 1 - 

+ Qa (T^) [ R 0 ©b-1 - <*A (Ta)] T A t (vQ rA )*= 0 

(III-3a) 

(^) + «B ( T B> ex pry- A -7 B )*T Ai (vQ ,A ) 

-[iB * T F> "4 < T B>] V < T B’ T > ('a,' 5 r m k Feq. < T >) 

+ K o%[ T A 1 <»'m k B > < T B' T >] * j - 0 (in-3b) 


[l + exp (7 a - 7 B )*] ( v Q 7A ) * + V ( t b.T) p? m k Feq. ( T > = 0 ( ni ~ 3c ) 

where 

(vQ 7A ) = py m [exp (y B - 7 A )*4 1} (Tg,T) - k£? (T^.T)] (IH-3d) 



(III-3e) 
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The above equations are sufficiently complex that it is worthwhile to alter 
their form before considering a procedure for their solution. To simplify nota- 
tion we first introduce the definitions 



T ° ( pr m k Feq.) 1 
T b s (^mVb-i)" 1 

and the ratio of the rates given by equations (11-39) 


(III-4a) 


(III-4b) 


(IH-4c) 


< T A.T) _ IVii <4 <tb) exp / fbn 

4* (T|,T) ■ Q a <t£ V T A 


Q'b ( t b> 

■3 A < t A> 


exp 




(in-5) 


, / e A 2 \ 

This last equation results after the substitution k^_ 1 = exp I - • 

With the above quantities we rewrite the equations as follows: 

"j — 

< *A 1 (TA>0 Al (T,T J J) + QA 1 (TA^« S >^)jn Al (T.T^l-^<2_J 


'3A< T A)[ E o e b-,-lA( T A)J [ t A 1 ('' , 3 ,A ) . * 


(III-6a) 



( % CT.T® [l - ^j>] + <6 (T B *) ^ <-?|> 



“ [Qb < T |) " q B ( t b>] V ( T B’ T ) 



Ro©A 2 {TAi\ 
+ \%) 


1 - exp 




= 0 


(in-6b) 


[1 + exp (-7^)] (vQ 7a ) 


M* v ( t b> t > 

vQ a / + = 0 


(in-6c) 


where 



exp(y^) / t aA 
T AiQB(TB) V T b ) 


( k ( F } (T^,T) 

( k^T^.T) 


exp (-y 



(III-6d) 


exp(y^) / t aA 
t A^B ( t B^ \ T b / 


QB(Tg) r /_!_ A 
QA(^) eXP [ A 2 W“t] 


fb _* 

~ y R 


(m-6e) 


This system of equations, although still appearing quite complicated, can readily 
be solved numerically by the following procedure: First obtain an explicit rep- 
resentation for y-p as a function of the temperature T and the unknown quan- 
tities and Tg . Such a function is found by solving the quadratic equation 
in 7 g that results from the substitution of equation (III-6d) into equation (III-6c). 
The result is 



' 2AC 
B + VB 2 + 4AC 


(III-7a) 


where the identities 

AC = k^ (T^,T)/4 Vb,T) 


(III-7b) 


B = 




Q'b( t b) V(Tb,T) - AC 


(in-7c) 



are introduced to simplify notation. To obtain the above formulation, the 
root of the quadratic equation containing the negative radical was discarded 
since > 0 . The remaining root with the positive radical, of the form 
( -B + VB 2 + 4AC )/2 , was modified to obtain the equivalent form given by the 
equation above so as to avoid the numerical difficulties that occur due to sub- 
traction of the quantities B and VB 2 + 4AC when their values are very 
nearly equal. 

Having an explicit representation for y-^ simplifies the problem of solving 
equation (ITt-6) , since the problem may now be considered as one of solving a 
pair of coupled transcendental equations for the two unknowns Tj£ and Tg . 
Equation (III-6a), indicated symbolically by 

ft [ T. Tj£, Tg, 7 r (T, T I, T|)] = 0 (III-8a) 

is one of the equations and equation (III-6b), written 


f 2 [ T, T^, Tg, (T, t £, Tg)] = 0 (in-8b) 

is the other. The second equation, when solved numerically on the assumption 
that values of T and are given, may also be written 

Tg = T B (TA,T) (III-9a) 

When substituted into equation (III-6a) this yields the equation 

fi {T, T^, Tg (T|, T) , 7g [T, T^, Tg (T^, T)]} = 0 (HI-9b) 

This equation may then be solved numerically for T^ once T is specified. 
Although this outlined procedure contains embedded iteration loops that are 
very inefficient from the computational time standpoint, it was found satis- 
factory since only a limited number of cases were run. About 8 seconds were 



required with an IBM 7090/7094 DCS computer system to obtain a solution for 
one value of the temperature T . Both equations (lit- 9b) and (m-9a) were 
solved by means of a nested pair of iterations, which were based on the con- 
ventional Newton-Baphson method. It should also be pointed out that the order 
in which the equations are solved (as implied above) is important since the 
procedure of solving equation (HI -8a), assuming values for T and Tb , may 
not yield a solution unless the iteration procedure is started with relatively 
accurate bounds and initial guesses for T A and Tg . 


It is worthwhile to discuss equations (HI-6) briefly. One observes that 
these equations are independent of density, and thus their solutions will also 
be independent of density. The equations depend parametrically on the func- 
tions A(T), t a /r b , T Ai /rg , and t a / t d • Since r is inversely proportional 
to pressure, the latter three ratios are independent of pressure and are func- 
tions only of the temperature T . Also embedded in equations (IH-6) are the 
separation constants a and b and the constant parameter Tg . As a result 
of the procedure used to evaluate the transition probabilities (see Appendix B), 
the temperature -dependent quantities listed above, except , can be related 
to X(T) (see Appendix B) as follows: 


r A 



t a 



(b-a)X(T)- a 


Ea, 

e A a 


1 - exp 



1 - exp 




(Hl-lOa) 


(III- 10b) 


and equation (11-73) relates the dependence of A(T) on X(T) . Also given in 
Appendix B are two suggested equations for evaluating X(T) , that is, Xr(T) 
and Xf(T) . Solutions that depend on each function were obtained in order 
to assess the sensitivity of the resulting solutions to this parameter. 


82 



Solutions will also be exhibited for the hypothetical case X(T) = 1 that would oc- 
cur if the transition probabilities were strictly those for a harmonic oscillator. 

The characteristic vibrational relaxation time associated with the 

lower vibrational states (for discussion see Appendix B) is taken to be that 
quantity measured experimentally by Camac (20); that is. 


= Maniac 


= m C l T 


i/e 


’ / 2228V1 I" /C 2 V / 3 

Ll-exp V"— ;J ex P[-^Tj 


(Hl-lla) 


where 


Ci = N 0 Ci 


= 7.22 x 10 16 cc/mole-sec-(°K) 1 / 6 

(III- lib) 

C t = 1.20 x 10 _7 cc/particle-sec-(°K) 1 / 6 

(HI-llc) 

C 2 = 1.04 x 10 7 ±30% °K 

(III- lid) 

To assess the sensitivity of the solutions to , two values of C 2 

are used: 

cl = 1.04 x 10 7 "K 

(IH-12a) 

C\ = 1.40 x 10 7 °K 

(m-i2b) 


Here c\ is the value suggested by Camac, and c\ is slightly larger than the 
value associated with the upper error bound of 30 percent reported by Camac 
(see eq. (II-lld)). 

The equilibrium dissociation rate k-p e( ^ contained in the characteristic dis- 
sociation time t d (see eq. (III-4b)) is evaluated in two ways for the results that 
will be shown. First, it is simply presumed that we have a T 1 / 2 pre -exponential 
factor (recall discussion of Collision Theory given earlier in this chapter). It 
then turns out that the vibrational coupling factor V* obtained from the quasi- 
steady solution effectively adds an additional temperature dependence. 




[ k Feq. <T>] MT = 9 * 00 x loW ex P (" V T ) (in-13b) 

In both of the above equations the constant coefficients were evaluated so as to 
agree with Camac and Vaughn* s experimentally determined rates at 
T = 4,000° K, as was done by Marrone and Treanor (58). 


The parameters a and b are chosen somewhat arbitrarily. The value of 
a (discussed in Appendix B) is such that the deviation of the approximate "rate- 
of-quantum -transfer tf relation is not greater than about a factor of 2 (when 
compared with similar relations obtained from Morse oscillator quantities) at 
the vibrational energy level corresponding to the energy E a . A value of a = 9 
approximately satisfies this criterion. The choice of the value for b will ulti- 
mately depend on how accurately the solutions can be made to agree with 
experimental data. The choices b = 16, 24, and 30 result in the fractional 
energies E^/Ej) = 0.544, 0.786, and 0.968, respectively. The latter choice 
yields unsatisfactory pre-exponential behavior for the effective dissociation 
rate, as will be shown, and is only briefly considered. It is to be emphasized 
that these values of b are arbitrary, but appear reasonable on the basis of 
the results given in references 48, 84, 94. 

A series of plots are given in figures III-l through IQ-12 showing the solu- 
tions obtained from equations (HI-6) . Separate curves are given in each plot, 
each associated with a different value for the separation constant b, for each 
of the quantities Tj£/T, Tg/T, V*, and k Feff _ (T)/exp k Feff> (4,000). 
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For the reader! s convenience Table HI~1 lists the values and identifies the par- 
ticular temperature-dependent quantities used for the solutions. For example, 
for figure ni-1 the parameters have the values Tg = 00 (actually, Tg = 10 23 ), 

X = X R , k Feq = 9 x 10 14 exp C 2 = 1.4 x 10 7 , and a = 9 . Associated with 

these values, figure Ill-la shows six curves, of which three give T^/T for 
b = 16, 24, and 30, while the other three similarly give Tg/T . Figure in-lb 
contains three curves for the coupling factor V*, with each curve again for a 
different value of the separation constant b . Figure UI-lc shows the pre- 
exponential temperature dependence of the dissociation rate constant. The 
circles and squares are experimentally measured values, the circles from 
Wray* s data (99) and the squares from Camac and Vaughn* s data (19), The 
solid-line curve, labeled C & V, is the dissociation rate computed from the 
formula given by Camac and Vaughn (19) (representing a least squares fit of 
his data), but appropriately normalized. Since Wray (99) gave a formula with- 
out any pre -exponential temperature dependence, the horizontal line associated 
with 10° would be appropriate for his formula. The dotted-line curves are 
guides from which one may estimate the functional behavior of the various 
curves regarding temperature dependence. The lower dotted line represents 
a T” 2 pre-exponential temperature dependence, and the upper dotted curve 
corresponds to T~^ 2 . The two solid lines labeled b = 16 and b = 24 repre- 
sent the pre -exponential temperature dependence as obtained from the present 
quasi-steady solutions. 

An important point evident from figure III-lc is that the experimental 
values have less than the T~ 2 dependence obtained by Marroneand Treanor (58). 
Although Camac and Vaughn fit their data with a function that varies as T” 1 , it 
would appear from Wray* s higher temperature data that a variation of T” 1 / 2 
would be more reasonable. It is to be emphasized that a different normaliza- 
tion would cause a shift of the dotted T~ 1//2 curve and that a different place- 
ment of this curve can be made. Because of the scatter in the experimental 
data the most appropriate factor is still uncertain. The published data, how- 
ever, do not show the same degree of scatter because the exponential 
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Table m-1 

Parameters associated with the quasi-steady solutions, 
figures m-1 through III- 12. 


Figure no. 

T u 

X 

k Feq. 

c 2 

b 

a 

m-i 

oo 

R 

MT 

(2) 

16,24,30 

9 

m -2 

QO 

R 

WR 

(2) 

16,24,30 

9 

m-3 

©d/3 

F 

MT 

(2) 

16,24,30 

9 

m-4 

©d/3 

F 

WR 

(2) 

16,24,30 

9 

m-5 

©d/io 

F 

MT 

(2) 

16,24,30 

9 

m-6 

©d/io 

F 

WR 

(2) 

16,24,30 

9 

m-7 

© d /io 

F 

MT 

(1) 

16,24,30 

9 

m-8 

©d/20 

F 

MT 

(1) 

16,24,30 

9 

m-9 

0j)/5O 

F 

MT 

(1) 

16,24,30 

9 

m-io 

© d /100 

F 

MT 

(1) 

16,24,30 

9 

m-n 

— 

F 

MT 

(1).(2) 

16,24 

9 

m-i2 

e D /io 

X (T) = 1 

MT 

(1) 

16,24,30 

9 


X: R see eq. (B14c) 

F see eq. (B23b) 

^Feq. : see ec l- (ni-13a) 

MT see eq. (in-13b) 

C 2 : 1 see eq. (IH-12a) 

2 see eq. (HI-12b) 

factor exp (-©pj/kT), which actually dominates the measurements, is included. 
The pre -exponential factor compared with the exponential factor represents a 
relatively small variation. 

The problem of determining the appropriate values for the parameters em- 
bedded in the equations is difficult because of the relative uncertainty of the ex- 
perimental data, the limited amount of such data, and the rather large number of 
parameters to be evaluated. The essential features exhibited by the solutions are 

1. The relative values of T^/T and Tg/T, 

2. The temperature dependence of V*, 

3. A comparison of the pre -exponential temperature dependence of the 
effective dissociation rate constant. 



It will turn out that these features together with certain facts obtained by other 
investigators will be sufficient for estimating appropriate values of the embed- 
ded parameters. 

It is known that one effect of dissociation is to perturb the population in the 
upper states to values less than those appropriate for a Boltzmann population 
distribution (48, 84, 94). We would therefore expect that Tjj/T < TjJ/T, espe- 
cially at the higher temperatures. Furthermore, we seek values for the 
parameters that will yield a pre-exponential temperature dependence that 
varies roughly as T" 1 / 2 . Searching through the figures we note that the com- 
bination of parameters associated with figures Itl-l may be discarded since, in 
this case, Tg/T > T^/T (see fig. Ill-la) and, in addition, the pre -exponential 
temperature factor varies too weakly with temperature (see fig. III-lc). For 
the case in figure in -2 kp e ^ (T)j^rp is replaced by ^-peq (^WR anc * 
results show that the perturbing effect of dissociation is now accounted for 
(see fig. III-2a). The pre -exponential temperature dependence, however, does 
not agree with the experimental result (see fig. III-2c). One observes that the 
T 1 / 2 factor contained in k Fe( ^ (T)^^ produces effectively a positive pre- 
exponential temperature exponent rather than a negative exponent as is implied 
by the experimental data. Reducing T-q (see, e.g., figs. HI-2, m-4, HI-6) 
reduces the positive exponent, but not to such an extent as to effect satisfactory 
agreement with the experimental data. However, we also observe that reduc- 
ing Tjj reduces T-g/T relative to T^/T and increases the pre -exponential 
variation with temperature (cf. figs. HI-1 through HI-6). 

As explained earlier, the parameter Tjj is a measure of how rapidly the 
dissociation probability drops off for the lower vibrational levels (see 
Appendix C for discussion). The larger values of T u correspond to a rela- 
tively equal probability of dissociation from all the B-state levels and the 
smaller values to a relatively greater probability of dissociation from the 
uppermost levels. The quasi-steady solutions seemingly imply that a large 
pre -exponential temperature dependence requires small values of Tu* This 
is consistent with Treanor and Marrone* s result (91, 93) in that T^j = °o 


(nonpreferential dissociation model) resulted in a T _1 pre -exponential tempera- 
ture dependence while T|j = ©d /6 (58) (preferential dissociation model) 
resulted in roughly a T" 2 temperature dependence. Aside from including the 
vibrational temperatures T A and T B and the effect of anharmonicity, the 
model differs from Marrone and Treanorjs preferential dissociation model (58) 
in that here only the vibrationally excited molecules, characterized as being in 
the B-states, are allowed to dissociate. The effect of decreasing b would 
result (although such calculations are not shown) in a p re -exponential tempera- 
ture factor as large as theirs. 

The discussion thus far has centered on figures HI-1 through IE-6, which 
show the solutions associated with the value of C 2 given by equation (III-12b) 
(see also table in-1), that is, which correspond to a relatively large character- 
istic vibrational relaxation time r A 1 - Decreasing this parameter by use of the 
expression for C 2 suggested by Camac and Vaughn (i.e., using eq. (III-12a)) 
increases the values of T3/T relative to T^/T and also decreases the pre- 
exponential temperature factor. It thus has an effect opposite to that desired 
on the basis of the experimental data. This may be countered by further reduc- 
tion in the value of Tjj and the appropriate effects are illustrated in fig- 
ures HI-8, III-9, and IE-10. However, comparing figures IE-9 and III-10, we 
observe that there appears to be a small value of T u below which smaller 
values have little effect on the pre-exponential temperature dependence. That 
this is true is illustrated in figures IE-11, wherein the coupling factor V* is 
plotted as a function of T|j for T = 10,000° and 20,000° K,and for b = 16 and 24. 
Comparing figures El-lla with IE-llb we observe that a change in the value of 
C 2 of only about 36 percent (a factor of between 2 and 3 in the characteristic 
vibrational relaxation time over the range of temperatures of interest) 

has a rather large effect on the coupling factor or the pre -exponential tempera- 
ture dependence. These figures also show that any incremental change in T-q , 
when this parameter is greater than 10 5 or less than 2 x 10 3 , has little effect on 
the coupling factor V* (also see discussion in Appendix C). In the former 
limiting case we have nonpreferential dissociation from any B-state level and 
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in the latter case highly preferential dissociation — in essence, dissociation only 
from the very uppermost level. We note further that the latter limiting case 
yields the smallest values for the vibrational coupling factor V* for any given 
temperature T and hence the greatest pre -exponential temperature dependence. 



It then follows that the largest pre -exponential temperature variation 
obtainable with the assigned values of b is that shown in either figure IH-9 
or ni-10. A greater pre -exponential temperature dependence will require 
smaller values of b . The results shown, however, are within the experimental 
data scatter; hence, 16 is the largest value of b that will yield favorable agree- 
ment with experiment. 


No comparisons are shown from which one can estimate the effect of using 
the two specified values for X(T) (see eqs. (B28a) and (B28b)). It turns out 
that the solutions are little affected by whether Xp (T) or Xg(T) is used for 
X(T) . The greatest effect is a difference in the values of T^/T and Tg/T , 
amounting to about 2 percent, and occurs at the lower temperatures where the 
differences between the values of Xp and X R are greatest (see Appendix B) . 
At the higher temperatures the differences were observed to be less than 
1/2 percent. One additional effect of X is worth pointing out. The reader 
will observe that Tg/T crosses T^/T for a few of the solutions (see, e.g., 
fig. IH-7a). This crossover occurs because the approximations introduced for 
Xg(T) or X F (T) in Appendix B are inaccurate at the lower temperatures, and 
hence it is probable that Tg is also in error at those temperatures. It is to 
be expected that an additional segment (that is, an additional separation com- 
parable to the a-level separation) would eliminate the crossover and improve 
the value computed for t b- It will be shown in the next chapter that such an 
additional segment may be advisable for other reasons also. On the other hand, 
the effect of such improvements is not expected to be large since the relative 
molecular population of the states that are affected is small (see discussion in 
Appendix B and also the discussion of the component derivative terms in the next 
chapter). This is further attested by the fact that the use of either Xg or Xp 
has only a small effect on the solutions. This is not to say that the effects of 
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anharmonicity may be ignored. One observes by comparing figures m-7 and 
m-12 that differences do occur depending on whether one includes or ignores 
the effect of anharmonicity on the transition rate, although the differences for 
the case of a normal-shock wave may be sufficiently small to be unobservable 
experimentally. It will be shown in the next chapter that for nozzle flow, the 
effects of including anharmonicity may be large and should not be ignored. 

The investigation in this chapter has provided insight on the effect of the 
parameters and the effect of incremental changes in the values of the param- 
eters on the quasi-steady solution. It should be apparent that the values of b 
corresponding to energies that are greatly above the midpoint of the 
vibrational-energy potential well do not provide a sufficiently large pre- 
exponential temperature variation as compared to that found experimentally. 
For each b below the midpoint value, the parameter Tjj can be adjusted to 
provide the proper pre -exponential temperature variation. Consequently, a 
unique choice of values for the embedded parameters cannot be found simply 
on the basis of the comparisons with the quasi-steady solutions. Additional 
comparisons must be made and this, in part, is the purpose of the investiga- 
tions discussed in the next chapter. 

Solutions will be obtained that show the details of the relaxation to equilib- 
rium behind a normal shock wave. The effect of changes in C 2 , b , and Ty on 
the relaxation behavior will be investigated further, particularly, the effect of 
these variables on the induction time. Finally, an example will be given for 
nozzle flow together with a discussion of the distinguishing features associated 
with such flow. 



Temperature, T °K/IOOO 


(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 

Figure III— 1 . Variables obtained from quasi-steady solution versus temperature 
(associated parameters are listed along first row table III— 1) 







Temperature, T °K/IOOO 


(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 

Figure III— 2 . Variables obtained from quasi-steady solution versus temperature 
(associated parameters are listed along second row of table III-l) 




Pre - exponential temperature factor 




(c) Pre-exponential factor 
Figure III— 2 Continued 
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Coupling factor, V* Vibrational temperatures, dimensionless 



Temperature, T, °K/IOOO 

(c) Pre-exponential factor 
Figure III— 3 Continued 
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2 - 


b 



Temperature, T °K/1000 

(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 


Figure III-4. Variables obtained from quasi-steady solution versus temperature 
(associated parameters are listed along fourth row of table III-l) 






(c) Pre-exponential factor 
Figure III— 4 Continued 



Temperature, T °K/IOOO 


(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 

Figure IE-5. Variables obtained from quasi-steady solution versus temperature 
(associated parameters are listed along fifth row of table III— 1) 




Temperature, T, °K/IOOO 


(c) Pre-exponential factor 
Figure III-5 Continued 




Coupling factor, V* Vibrational temperatures, dimensionless 




Temperature, T, °K/IOOO 


(c) Pre-exponential factor 
Figure III— 6 Continued 









Pre - exponential temperature factor 



2 6 10 14 18 22 26 

Temperature, T, °K/IOOO 


(c) Pre- exponential factor 
Figure III— 7 Continued 
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Coupling factor, V* Vibrational temperatures, dimensionless 



Temperature, T °K/IOOO 

(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 



Figure III— 8 . Variables obtained from quasi-steady solution versus temperature 
(associated parameters are listed along eighth row of table III-l) 
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Pre -exponential temperature factor 




(c) Pre-exponential factor 
Figure III— 8 Continued 
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(c) Pre-exponential factor 
Figure III-9 Continued 
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Coupling factor, V Vibrational temperatures, dimensionless 




Pre- exponential temperature factor 



2 6 10 14 18 22 26 


Temperature, T, “K/IOOO 

(c) Pre-exponential factor 
Figure III— 10 Continued 
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Coupling factor, 
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(a) Relatively rapid vibrational relaxation time 

(b) "Correct" characteristic vibrational relaxation time 

Figure III-ll. Vibrational coupling factor obtained from quasi-steady solution 
versus temperature (associated parameters are listed along 
eleventh row of table III— 1 ) 
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Coupling factor, V Vibrational temperatures, dimensionless 




(a) Dimensionless vibrational temperatures 

(b) Vibrational coupling factor 

Figure III— 12. Variables obtained from quasi-steady solution versus tempera- 
ture (associated parameters are listed along twelfth row of 
table m-1) 


Temperature, T °K/IOOO 


(c) Pre-exponential factor 
Figure III-12 Continued 




CHAPTER IV 


COMPLETE SOLUTIONS OF THE MODEL EQUATIONS 
IV- A. Introductory Comments 

The preceding chapter was based on a study of the effects produced by 
changing the values of various parameters contained within the model equations. 
The investigation was greatly facilitated by the availability of easily obtained 
special solutions. Just as important, these special solutions had a particular 
significance in that they represent the conditions relevant to a comparison with 
many of the experimental rate observations. The solutions are particular to 
the quasi-steady zone behind a normal shock wave. The purpose of the investi- 
gation was primarily directed toward an assignment of values for the param- 
eters. It is not possible, however, to evaluate the parameters in the model 
equations with certainty purely on the basis of comparisons with the quasi- 
steady solutions. As was observed, one may obtain equivalent results using 
different combinations of values for the parameters. In addition, because of 
the scatter in the experimental data, there is a broad range of values that one 
may assign to any one of the parameters and still effect satisfactory results. 
Additional comparisons are therefore advisable, and this is one purpose of the 
discussion in this chapter. In particular, we take advantage of the fact that 
measurements have been made on the delay time (incubation or induction time) 
that occurs before the onset of dissociation (99). These experimental measure- 
ments provide another and completely different basis for comparison. In 
addition, several measurements have been made of the characteristic vibrational 
relaxation times associated with expanding (nozzle) flows (98). However, 
because there is a considerable variation of the reported values of this quantity 
(44, 45, 61, 78, 79, 80, 98) (in some cases the values differ by several orders 
of magnitude), quantitative comparisons will not be made. It will be shown 
that the model yields qualitative agreement with these studies and, for the 
more recent references (98), the agreement appears quite good. 



The essential difference between this and the preceding chapter is that 
here the discussion will involve the simultaneous solution of the model rate 
and fluid-flow equations . A rather complicated numerical integration procedure 
was required to obtain these solutions and the development of this procedure 
was a major task. A discussion of the numerical aspects of the problem is 
given in Appendix E, and the explicit form of the equations that are solved is 
given in Appendix F. 

After the final comparisons are made, a single set of values can be defined 
for the embedded parameters that effect reasonable agreement of the model 
with available experimental data. The model, in this respect, is complete. 

As a result of the formalism reported in Chapter II and used to derive the 
equations, the various terms constituting the rate equations have physical sig- 
nificance, and their investigation therefore is worthwhile. Such an investigation 
yields insight as to why differences should be expected when comparing, for 
example, characteristic vibrational relaxation times obtained from normal- 
shock wave and nozzle data. 

In the section that follows, the discussion will first concern the shock 
structure that results from the complete solution of the model equations. 
Initially no emphasis will be made on the effects that result from any particular 
choice of parameter values. Qualitative features of the shock structure appear 
to be due to the structure of the equations regardless of specific values of the 
parameters. The relationships of the various terms to these features will first 
be discussed to promote an understanding of the vibration-dissociation coupling 
process. After the discussion on the shock structure, attention will be given to 
the differences caused by changing the values for the parameters. The effects 
of the parameters on the induction time will also be considered. Later sections 
provide a description of the events occurring in nozzle flow. An explanation is 
also given of the significant differences observed between shock-wave relaxa- 
tion processes and the relaxation processes associated with expanding flows. 



IV-B. Normal- Shock Wave Solutions 



The results of this section are obtained by setting the derivatives dA/dx in 
equations (F16) equal to zero and then integrating the resulting equations ac- 
cording to the numerical integration procedure described in Appendix E. (The 
explicit method is used for about the first 125 computational points, and the 
implicit method is used thereafter.) The conditions downstream of the shock, 
which provide the initial values for the integration procedure, are obtained by 
an iteration procedure similar to that given by Vincenti and Kruger (p. 179 in 

97). Here the nonequilibrium variables 7 , 7,7 , T , and T are assumed to 

A. d a a Jj 

be frozen across the viscous shock. The gas mixture used for these solutions 
was 4 percent molecular oxygen and 96 percent argon. The effects of consider- 
ing molecules other than argon as a colliding partner may therefore be taken as 
negligible (that is, the fractions y^/y^ and y Q /y m in eqs. (H-113) may be con- 
sidered small and only the species argon considered for the summation index r). 

A sampling of the results is given in figures IV-1 through IV-7. For the 
reader’s convenience, table IV-1 lists the values for the temperature-dependent 
parameters in the cases plotted. The cases illustrate the effects occurring for 
several values of the parameters and for shock waves having values for the 
downstream temperature T 2 (kinetic temperature directly behind the shock 
wave) of 6,677° and 16,570° K. One sample case is plotted for the lower tem- 
perature purely to illustrate the differences that occur for weaker shocks. The 
abscissa scale for these figures corresponds to the time after passage of the 
shock (laboratory time) appropriately adjusted^ so as to remove the dependence 


f The abscissa scale is computed according to the formula 


(^lab) adjusted “ ^lab x ( n 2A* 28xl0 17 ) (IV 1) 

where n 2 = p 2 /kT 2 , k is Boltzmann’s constant, and the subscript 2 denotes 
values immediately downstream of the shock wave (see, e.g., p. 56 in 10 for 
discussion of t^). This adjustment is sufficiently accurate since the density 
behind the shock is nearly constant. 
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Table IV- 1 


Parameters Associated With Complete Solutions, Figures IV- 1 Through IV- 7 


Figure no. 

T 2 

T u 

X 

IV- 1 

16,570 

©jj/lO 

F 

IV-2 

16,570 

©d/10 

F 

IV- 3 

16,570 

©d/10 

F 

IV-4 

16,570 

©d/10 

F 

IV- 5 

16,570 

© d /50 

F 

rv-6 

6,677 

e D /io 

F 

IV- 7 

16,570 

©d/10 

X = 1 


k Feq. 

©2 

b 

a 

Quasi- steady 
Fig. no. 

MT 

(1) 

16 

9 

m-7 

MT 

(2) 

16 

9 

m-5 

MT 

(1) 

24 

9 

m-7 

MT 

(2) 

24 

9 

m-5 

MT 

(1) 

16 

9 

m-9 

MT 

(1) 

16 

9 

m-7 

MT 

(1) 

16 

9 

in -12 



X: 

k Feq. : 

CV 


F 

MT 


( 1 ) 

( 2 ) 


see eq. (B23b) 
see eq. (III-13b) 
see eq. (Ill- 12a) 
see eq. (Ill- 12b) 


of this scale on density. In this manner one can compare the effects occurring 
for different shock- wave solutions regardless of the density behind the shock. 
This procedure parallels that used by Wray (99) and in effect relates the effect 
of relaxation to a common density behind the shock of 1.28 x 10 17 particles/cc. 

Before beginning the discussion of the solutions, a few comments of a 
general nature should be made concerning these figures. First, the reader is 
cautioned to observe that the time scales may be different for the different plots. 
The figures were generated automatically with a General Dynamics S-C 4020 
Electronic Plotter (for a description of the system used, see 50). The scaling 
parameters were also obtained automatically and since the integrations were 
not always terminated at the same relative reference time, some plots encom- 
pass longer times after passage of the shock wave than others (e.g., compare 
figs. IV-1 and IV-2). The figures show the complete shock profiles; for the 
longer times, once equilibrium is reached, no further effects are observed and 
the values of the dependent variables remain constant. The plots show an 
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erratic behavior in the first transient zone that in some cases is of greater 
magnitude than for others (e.g., compare figs. IV- 1 and IV-6). This behavior 
depends purely on the numerical integration procedure and represents the 
effects of numerical instability. It has no effect on the final values of the inte- 
grated variables. One has full control of the amplitude of the instability, and 
the effects are shown only to lend additional weight to the discussion on the 
numerical integration procedure given in Appendix E. For the present dis- 
cussion one may ignore the instability phenomenon, since they do not affect the 
final results and since they are of no consequence to a physical description of 
the results. The correct values of the variables are made evident by a smooth 
curve fitted through the lower locus of computed points. 


Given in each figure are plots of the variables T ,T A , T B,r A ,y B , r a , and the 
ratio (this latter quantity becomes constant in the quasi-steady zone). In 

figures IV-1,3,5,6 the vibrational coupling factors are also plotted and in fig- 
ures IV-1 and IV-7 are shown, in addition, the derivatives and their component 
terms. 


(IV-B-1). Shock- Wave Structure 

We first consider in detail the shock structure resulting from a solution of 
the model equations, ignoring any consideration of effects peculiar to a particu- 
lar choice for the embedded parameters. Figures IV-1, for example, give a 
complete picture of the relaxation processes (these processes are associated 
with a single set of parameters having the values listed along the first row in 
Table IV-1). The three characteristic zones mentioned previously are labeled 
in this figure by the numerals 1,2, and 3 corresponding, respectively, to the 
transient, the quasi-steady, and the final relaxation zones. In the first region 
the vibrational temperatures, T^ and T B’ and the upper-state and atom con- 
centrations, and 7 a 5 rapidly increase in value, attaining their quasi-steady 
values in about a quarter of a microsecond; that is, in the transient region, the 
upper vibrational states are rapidly populated, thermally "heated" and, at the 
same time, partially dissociated as evidenced by the increasing atom concen- 
tration (recall that only the B molecules are allowed to dissociate). In the 
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Adjusted shock passage time, sec 


(a) Temperatures 

Figure IV-1. Variables obtained from shock solution versus adjusted 
laboratory time (associated parameters are listed along 
first row of table IV-1) 


Concentration, moles per unit total mass 



(b) Concentrations 


Figure IV-1 Continued 
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Concentration, moles per unit total mass 



(c) Concentration, Lower States and Molecular 


Figure IV-1 Continued 
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Vibration coupling factor 




Adjusted shock passage time, sec 

(d) Vibrational coupling factor 
Figure IV- 1 Continued 
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transient zone, the lower-state molecular concentration remains relatively 
constant, and the vibrational coupling factor rapidly increases from its initially 
small value. 

The quasi-steady zone is indicated in these figures, and the values obtained 
from the quasi-steady solutions of the preceding chapter are given!" These 
values are labeled by the variables superscripted with an asterisk. The quantity 
V* compares accurately with the quasi- steady value obtained here, as can be 
seen in figure IV-ld. There is a slight discrepancy of about 2 percent, however, 
between the results for both T^ and T^ as computed in the previous chapter 
and as found here. This difference is attributed to the fact that the temperature 
T is not exactly constant as was assumed in the previous chapter, but changes 
about 10 percent as a result of the effects of dissociation (note that the tempera- 
ture T varies only in the quasi-steady zone). The dissociation energy has a 
relatively large effect on the kinetic temperature even though the molecular 
concentration is small (initially 4 percent of the total mixture). The density and 
velocity, however, are not so greatly affected by dissociation. These quantities 
are nearly constant and are therefore not plotted. (They varied about 4 percent 
as a result of dissociation.) One notes that in the quasi-steady zone the ratio 
7 b / and the temperatures T^ and T-^ are nearly constant. The basic 
quasi-steady assumptions are thus in fact valid in these cases. It is also 
interesting that, although the concentrations y ^ and are varying, their 
ratio is indeed constant. 

An additional feature characterizing the quasi-steady zone is that in this 
zone the effects of dissociation become appreciable, as is evidenced by the 
rapid decrease in the molecular concentration and by the effect on T. The 
entire molecular concentration is approximately contained in the states 
denoted "A" as a result of the relatively large difference in the concentration 


t The temperatures T 2 and the applicable quasi-steady solutions are listed in 
Table IV-1. With this information the respective quasi-steady values for the 
variables are readily obtained from the figures of the preceding chapter. 
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variables 7 ^ and y^ . This is exemplified further in figure IV-lc, where 
the lower- state concentration y A is plotted along with the total molecular con- 
centration y = y A + y B . It is also worthwhile to point out that the tempera- 
ture Tg initially increases in the transient region at a more rapid rate than 
> but near the end of this region T A overtakes Tg . In the quasi-steady 
region the effect of dissociation is to cause t b to be lower than T A . 

As regards the passage time of the latter two zones, for the case exempli- 
fied by figure IV-1 the quasi-steady zone persists about 8 yusec or about 30 
times longer than the transient region. The final relaxation region persists 
about 7 fji sec, during which time final relaxation to equilibrium occurs. From 
these results one can appreciate why many of the experimental observations 
are conducted during the rather long-time passage of zone 2. Although the final 
relaxation zone also requires a relatively long passage time, experimental 
measurements in this region are often hampered by other effects such as, for 
example, the arrival of a driver gas (shock-tube studies, see, e.g., 10, 70, 86). 


We have so far been discussing the structure of a shock wave according to 
the behavior of the basic dependent variables. This discussion was based on a 
single !T complete n solution. In figures IV-le through IV-li are plotted the com- 
plete derivatives and their separate terms associated with this solution. To 
facilitate a description of these figures the following notation is introduced for 
the terms that appear in equations (11-113): 


f A = p y mkffi (T A’T) 

^ i v 


(IV-2a) 


q'a . £> * q (Tb , T) 

Ql = 2 la V(T B ,T) k Feq m/v 

• Q 



-2P 2 y 


m 


k FeqF> 

k b( t ) 


/v 


(IV-2b) 

(IV-3a) 


(IV-3b) 
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This notation is consistent with that used by Lomax and Bailey (53) where, for 
example, Q A denotes the production of lower-state molecules, (T A = in7 A ), 
per unit distance along a stream line and has c.g.s. units of cm" 1 . The reader 
should not confuse the superscripted Q-quantities with the earlier Q-quantities 
without superscripts, the notation used for the partition functions (cf., e.g., 
eqs. (IV- 5 ) ). The Q-quantities with superscripts only represent the sum of 
the Q-quantities having like superscripts and integer subscripts (note eq. (IV-7)). 
Except for the factor l/v and also the factor /3 in equation (IV-6b), the reader 
can identify the above quantities as being the individual terms in equations 
(3T-113) of Chapter II. The symbols not defined here have the same meaning as 
in Chapter II. The reason for introducing the factor ft will be explained 
shortly. 

In figure IV- le are plotted the separate terms of the derivative for T A as 
well as the derivative itself. The values for these terms were obtained at the 
same time that the variables shown in figures IV- la through IV- Id were found. 
The sign associated with the separate terms, as well as their sum, is indicated 
parenthetically after the labeling symbol (the logarithm of the absolute values 
are plotted, hence the sign must be indicated separately). The three zones are 
also labeled here for reference. In the transient region we see that the deriv- 
ative and the component terms have very small values (resulting in no effect 

'Y 

on r A ) and that the loss term (negative term) Q x A is dominant. This rela- 
tive behavior continues throughout the transient zone. Near the end of this 

region we observe that, although the gam term (positive term) overtakes the 

y. 

loss term in value, their difference Q , or algebraic sum, approaches a 
maximum value that remains relatively constant throughout the quasi-steady 
region. At the end of this region, or at the beginning of the final relaxation 
zone, the gain term approaches even closer the value of the loss term, as is 
exemplified by the decrease in the value of the derivative Finally, at 



Derivatives Q A , cm 




(e) Derivative, q7a. ( and component terms, qTA 


Figure IV-1 Continued 
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equilibrium the derivative is constant, but with a value sufficiently small so as 
to have only a negligible effect on the concentration variable y^ 

y 

The derivative Q A thus behaves in the manner expected. Since the deriv- 
ative contains only single gain and loss terms such behavior could have been 
inferred from an examination of the behavior of the variable y ^ (fig. IV-lc). 
Likewise an equivalent statement may be made about the behavior of y a (see 
figs. IV- lb and IV-lf), where the processes are different only to the extent that 

the signs of the dominant terms are different. There is, however, an additional 

y T T 

complexity associated with the derivatives Q Q A, and Q B that requires 
comment. 


y y 

The separate terms of Q ^ are plotted in figure IV-lg. The term Q 4 B , 
associated with the recombination effects, is at first negligibly small (initially 
nearly 40 orders of magnitude smaller than the next larger term) and remains 
so throughout both the transient and quasi-steady zones. (This same effect can 

Jy y 

be noted for all the terms Q/ a , and Q 5 B that represent the effects of 

atom recombination.) The dominant term in the transient region is the gain 

y_ j. 

term Q X B that represents the gain in y B due to losses from y A > 


y 

^The derivative Q A ( as we ll as the other derivatives) should become zero 
when the final relaxation to equilibrium has occurred. The fact that the deriva- 
tive has a non-zero small value, that is, small relative to the values of its sep- 
arate terms (four orders of magnitude less), indicates that the numerical 
integration is not exactly precise. In Appendix E it is pointed out that these 
solutions are only accurate to about 0.01 percent. This is further exemplified 
here by the fact that at equilibrium the difference between the separate terms 
is of this order. 

y y 

tit is interesting to observe that i s proportional to A (see eq. 

(IV-4a) ) and that in the transient zone each term is the dominant term of the 
derivative in which it appears. The term Q^B however, causes a rapid 
variation in the variable y ^ while the term Q^A has no effect on y A until 
T^ is sufficiently large. 
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Figure IV- 1 Continued 


129 













The depletion terms, in the rate equation for the upper-state population, are 

y y 

Q 2 b and Q 3 b . The latter term, associated_with molecular dissociation, has 

y y y 

only a very small effect on the derivative Q B (Q3 B ~10~ 5 Q B ), and continues 
to have only a small effect until the gain and loss due to the transitions ^ y ^ 

are nearly balanced. This occurs near the end of the transient region. In the 

y y y 

quasi-steady region the relative values of Q2 > and Q3B remain nearly 

y 

constant and all have an equally important net effect on the derivative Q B # 

y 

Note also that the derivative Q B is smaller than any of its component terms 

y 

(except the recombination term, Q 4 B , of course). The behavior in the quasi- 
steady region remains unchanged until the effects of atom recombination become 
important, after which the final relaxation takes place. In the final relaxation 

y 

zone the derivative Q B has a slightly lower value resulting in negligible 
changes in the value of This variable may thus be considered to have 

attained its final equilibrium value. 


There are a few essential features that will be pointed out for emphasis. 

We note in particular (1) that the effects of atom recombination are negligible 
in both the transient and quasi-steady zone; and (2) that in the quasi-steady 
region the transitions causing molecular excitation are balanced by the effects 
of dissociation. These points are not original to this work, but have been dis- 
cussed by previous researchers (see, e.g., 85, 94). Because of the relative 
simplicity of the present model, however, the effects are more readily exempli- 
fied here. 


The rate equations discussed thus far are those relevant to the changes of 
the population variables y A , y B , and y . Additional insight may be gained 
from a detailed examination of the energy relaxation processes from which, for 
example, the influence of anharmonicity effects may be explained. 

In figure IV-lh are plotted the derivative Q A ^ e associated com- 

ponent terms Q^ A - In this case the dominant term is the gain term Qj Ta 
that represents vibrational energy relaxation for the lowest energy states, 
denoted . The only loss term is Q 3 which has only a negligible effect 
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Adjusted shock passage time, sec 


T T 

(h) Derivative, Q A, an d component terms, Q 

Figure IV-1 Continued 
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on the derivative Q^A ^ the transient region. This term represents, in 
essence, the energy losses due to the transitions ^ . Perhaps the 

most interesting observation that can be made is that the effects of the vibra- 
tional relaxation processes Q 2 T A associated with the states denoted A 2 are 
also negligible in the transient region. This apparently means that in the 
transient region behind a normal- shock wave the effects of anharmonicity, 
although potentially large, because the characteristic relaxation time may 

be small, are not in fact important. A close examination of the term Q^A sug- 
gests the reasons for this assertion. First we note that the factor 
[l-4(T)/q (T a )] is small since the ratio A(T)/q^(T A ) is initially of order 

unity. The factor increases monotonically as this ratio decreases (T is constant 
and as T^ increases, ^ 2 ( T ^) also increases), reaching a maximum value as 
q A2 (T A ) reaches its maximum. Even when the effects of anharmonicity can be 
considered negligibly small, however, (that is, when A(T)~0 and the factor 
[ l-A(T)/q^ (T^)] is unity; see, e.g. , discussion after eq. (11-73) of Chapter II), 
we have in addition the Boltzmann population factor Qa 2 ^ T A^^a( T A^ * This 
factor is also small and has the effect of lessening the influence of the charac- 
teristic relaxation time , hence, the effect of Q 2 T A on the derivative 

Q T A . The noteworthy conclusion from this discussion is that the term Q a T A # 
which is analogous to the Landau-Teller expression (see discussion after eq. 
(II-58a) in Chapter II), governs the vibrational relaxation in the transient region 
for the A-states. 



The foregoing behavior of the component terms continues throughout the 
transient zone. In the quasi- steady zone, where T^ becomes constant (see 
fig. IV-la), the value of Q T A decreases abruptly and all the terms become 
important. (Note also that in the quasi- steady zone the sign of Q T A alternates 
between positive and negative values, as indicated by the oscillatory behavior 
of this quantity.) Actually, even in the quasi- steady region where Q 2 ^A is 
slightly greater than Q^A f it turns out that the effects of anharmonicity are 
also not appreciable (compare, e.g., solutions of the previous chapter, figs. 

HI-7 and m-12). 
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An explanation of the interplay of the various terms that constitute the 
rate equation for Q T B (fig. IV-li) is considerably more complex since 
more terms are involved and the processes are more complicated. One 
observes that there are only two gain terms, and Q^B • The latter 

term, however, has a negligible effect on Q^B i n the transient and quasi- 
steady regions behind a normal shock. Thus, only one gain term has any effect 
on the upper-state energy relaxation processes, and it turns out that this term 
when combined appropriately with all the other terms yields (when p = 1) 
values for the derivative Q T B that are slightly negative. This condition 
causes the system of equations that are solved for the region behind a normal 
shock wave to be inherently unstable. For this reason the factor p was 
included, and was so adjusted that the sum of all the separate terms yield a 
value of Q T B that is very slightly positive. In particular, p was arbitrarily 
adjusted so as to satisfy the condition 


|qi T b I /(Iq 2 Tb I + Us B \ + Iq? 5 ! + Iq 5 Tb I) - i = i x io - 5 (iv-8) 

This condition introduces a problem to the extent that we are, in essence, arbi- 
trarily defining the starting values for the derivative Q^B # jn practice, how- 
ever, this problem was not serious for the following reasons: (1) the factor p 
was always of order unity (never smaller than 0.96 in value for the cases inves- 
tigated), and imposing an uncertainty of this magnitude is certainly not great as 

(i) 

compared with other uncertainties (e.g., the value of k-^(Tg,T) contained in 
Q 2 T B). (2) The effect of this correction factor is quickly lost as the other terms 
become important. The introduction of the factor p was only required for the 
normal- shock problem. In obtaining the expanding flow, that is, the nozzle solu- 
tions that are to be discussed later, the factor was always set equal to unity. 

The addition of the factor p therefore has the effect of making the system 
of equations that are solved stable. It may be seen, however, that the rate 
equation Q^B still shows an unstable behavior (its value oscillates several 
orders of magnitude in value and, in addition, there are changes in its sign). 
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(i) Derivative, Q and component terms, Qj® 


Figure IV- 1 Continued 


135 







Also, in contrast to the preceding derivative Q T A ? the beginning of the quasi- 
steady zone is not marked by any distinctive change in the value of \\r e 

recall that in the preceding case the quasi-steady zone was indicated by a 
rather abrupt decrease in the value of (Note also that the derivative for 

T a is smooth throughout the entire flow except for a relatively small variation 
in the quasi-steady region.) 

The interpretation of the oscillatory behavior of Q^B is complex. It re- 
quires an understanding of the numerical integration procedure (see Appendix E), 
in particular, as regards the capability of the procedure to resolve the effects of 
"parasitic eigenvalues" (see also 53 and 53a). In addition, an understanding of 
the concept of "local equilibrium state" (see, e.g., p. 236 of 97) is also helpful. 

It turns out that the oscillatory behavior of the complete derivative is largely a 
characteristic of implicit methods and occurs when the effects of large- 
magnitude negative eigenvalues (parasitic eigenvalues) are important. Such 
behavior is more often exhibited in the case of nozzle flow (see, e.g., 29, 53a) 
that is discussed in a later section. 

To interpret this effect we reformulate Q^B i n the following manner (see 
p. 236 of 97): 

dT s = Tb _ Xo( T ’P, T A, T B,y A ’ y B , r a ) 

dx _Q VT (T, P, T a , T b , y A , y B , 7 a ) (IV-9) 

where r is some relaxation time that characterizes the relaxation of Q Tb , 
and X 0 is some function that is zero when r is zero. (The condition r = 0 
is a formal characterization of local equilibrium.) In the formal case of local 
equilibrium Q^B is equal to 0/0 and therefore is indeterminant. In such 
cases t b is customarily found by solving X G for Tg (see p. 237 in 97). 

For the case at hand such a procedure complicates the numerical computation 
and is not necessary since r is always non- zero, although of small relative 
value (i.e., vr corresponds to a relaxation distance that is small compared to 
other distances in the shock). One important aspect of the numerical integration 



fi 
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procedure is that it can be used in such cases. The oscillatory behavior of the 
derivative shown in figure IV-li has a negligible effect on the solution and the 
resulting solution corresponds to that which would be obtained by setting X Q 
equal to zero. An essential feature of these solutions is that the variables are 
smooth (see, e.g., fig. IV-la) although their derivatives show an oscillatory 
behavior. 



The state of local equilibrium for the temperature T B prevails throughout 
the entire region behind the normal shock, irrespective of the relaxation zone 
(this relative behavior was noted for all of the cases studied). The dissociation 
terms assume an increasing importance in the quasi-steady zone and, finally, 
as the effect of recombination becomes significant, we have final and complete 
equilibrium. 

It is of interest to observe the effect of in these solutions. In 

Chapter II this term was designated as "rate limiting" and it may be seen from 
figure IV-li that its effect is one of limiting the rate of "heating" of the upper 
vibrational energy levels. The term has the opposite sign from and 

serves to cancel the effect of that term. In the quasi-steady zone it also has 
the effect of lessening the influence of the dissociation term. Because of the 
inter coupling of all the terms and their complexity, it is not yet clear whether 
q 2 Tb is rate-limiting in the sense used by Pritchard or Bray (14, 71). 

The discussion until now has centered on a description of the shock-wave 
structure without concern for the differences that occur as a result of using 
different values for the embedded parameters. In figures IV-2 through IV-7 
the plotted solutions illustrate the results for other values of these parameters. 
The values used for these solutions are listed in Table IV-1. It is to be recalled 
that the abscissa scale represents the time after arrival of the shock wave 
normalized in such a manner that the effect of differing densities on the relaxa- 
tion processes is removed. 

We observe first, comparing figures IV-1 and IV-2, that the relative effects 
that occur as a result of increasing the value of the constant C 2 (the constant 
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Concentration, moles per unit total mass 




(c) Concentration, Lower States 
Figure IV-2 Continued 
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in TA t ; note eq. (HI-12) ) by about 40 percent are important to and Tg and 
of less importance to the concentrations. We recall also from Chapter HI that 
the effect of changing C 2 in the manner indicated was effectively to increase 
the characteristic vibrational relaxation time (actually, by a factor of about 2 
for a temperature T 2 of 16,500° K). This effect is observed here as roughly a 
factor of 2 increase in the passage time of the transient zone. Also, and 
T-g are different from the values found for the previous solutions in the quasi- 
steady zone. The effect on the transient time is illustrated even more explic- 
itly in figure IV-8, which will be discussed later. The increased transient 

passage time occurs because r A , in effect, defines the time scale for energy 

A i 

transfer to the upper levels; hence, increasing r A increases the transient- 
zone time. The rather large effect on the quasi-steady values for the tempera- 
tures T^ and T^ was already discussed in the previous chapter. The reader 
will also observe that there is also an effect on the time when final relaxation to 
equilibrium occurs, but the effect is not large (less than a factor of 2). No 
attempt will be made to find the correct values of T^ and T B other than to 
show what happens when different values for the parameters are used. However, 
the importance of having experimental measurements of the relative amounts of 
energy contained in the vibrationally excited energy states is dramatically 
apparent here (in particular, regarding fluid flows). Such measurements would 
be sensitive indicators for the problem of inferring more accurate values of the 
transition rates (or vibrational relaxation times). The author is not aware that 
any such measurements have been made. 

In figures IV-3 and IV-4 are plotted similar solutions that result after the 
value of the constant b is increased. One recalls that this constant denotes 
the separation between the vibrational energy levels having different tempera- 
tures T^ and Tg. The previous pair of solutions, figures IV-1 and IV-2, was 
for b = 16 = 0.544) and the present pair is for b = 24 (6^/0^ = 

0.786). A comparison of figures IV-1 with figures IV-3 or figures IV-2 with 
IV-4 will show that increasing b reduces the difference T^ - T^, as was 
noted in the previous chapter, and has an opposite effect to that observed from 
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(a) Temperatures 

Figure IV-3. Variables obtained from shock solution versus adjusted 
laboratory time (associated parameters are listed along 
third row of table IV- 1) 
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(c) Concentration, Lower States and Molecular 
Figure IV-3 Continued 
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Temperature, °K/IOOO 


18 




Adjusted shock passage time, sec 


(a) Temperatures 

Figure IV-4. Variables obtained from shock solution versus adjusted 
laboratory time (associated parameters are listed along 
fourth row of table IV- 1) 
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Concentration, moles per unit total mass 




(b) Concentrations 
Figure IV-4 Continued 
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Concentration, moles per unit total mass 
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lO -2 



Adjusted shock passage time, sec 

(c) Concentration, Lower States 
Figure IV-4 Continued 
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increasing Cg (and hence 7^). On the other hand, increasing b produces the 
same effect as increasing C 2 on the passage time of the various zones, as well 
as on the time required to reach equilibrium. In particular, these times become 
longer. 

A comparison of figures IV- 1 and IV-5 shows the result of decreasing the 
value of from © D /l0 to © D /50. In the latter case the embedded param- 
eters were considered best from the standpoint of the quasi-steady solutions. 

It is to be recalled that the effect of reducing T^j to e^/50 is in the direction 
of making dissociation "highly preferential" so that only the most highly vibra- 
tionally excited molecules dissociate. One notes that although the differences 
between figures IV-1 and IV-5 are difficult to discern, this change does increase 
the time required for relaxation to equilibrium (also see fig. IV-8). 

The effects occurring for weaker shock waves are illustrated in figure 
IV- 6. Here the parameters have values identical to those used for the solution 
illustrated in figure IV-1 except that the temperature T 2 is lower. We note 
that in this case a relatively long time is required for the passage of the 
transient zone, and in the quasi-steady zone the temperatures T, T^, and 
are all equal. Once the transient zone has passed the entire vibrational popula- 
tion may thus be considered Boltzmann at the temperature T. Also, the 
internal processes relax in a sufficiently rapid manner that as soon as dis- 
sociation begins (note fig. IV-6c) the perturbing effect of dissociation has no 
effect on the vibrational population distribution. 

One particular feature to be noted from the figures IV-1 through IV-6 is 
the comparison of the quasi- steady values of the variables T^, T^, and V 
with the values obtained from the previous chapter. The reasonable agreement 
that is observed in all cases lends support to the basic quasi-steady assump- 
tions that were introduced in the previous chapter in order that those solutions 
could be obtained. 

A comparison of figures IV-1 and IV-7 shows the relative effects that 
result when X = Xp and X = 1, respectively. This parameter accounts for 
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(b) Concentrations 
Figure IV-5 Continued 
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(a) Temperatures 

Figure IV-6. Variables obtained from shock solutions versus adjusted 
laboratory time (associated parameters are listed along 
sixth row of table IV- 1) 
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Concentration, moles per unit total mass 




(b) Concentrations 
Figure IV- 6 Continued 
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Adjusted shock passage time, sec 
(d) Vibrational coupling factor 
Figure IV-6 Continued 
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Temperature, °K/IOOO 




Adjusted shock passage time, sec 


(a) Temperatures 

figure IV-7. Variables obtained from shock solutions versus adjusted 
laboratory time (associated parameters are listed along 
seventh row of table IV-1) 
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Concentration moles per unit total mass 



Adjusted shock passage time, sec 



(b) Concentrations 
Figure IV- 7 Continued 
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(d) Vibrational coupling factor 
Figure IV- 7 Continued 
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Derivatives Q cm 
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(f) Derivative, Q a , and component terms, Q j a 
Figure IV- 7 Continued 


163 



r 



y y 

(g) Derivative, Q B , and component terms, Q^B 
Figure IV-7 Continued 
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(h) Derivative, Q A, and component terms, Qj A 

Figure IV-7 Continued 
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Derivatives Q B t °K/cm 




T T 

(i) Derivative, Q and component terms, Qj® 

Figure IV-7 Continued 
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the effects of anharmonicity in the transition rates; in the case that such effects 
are negligible this parameter has a value that is unity. In figures IV-7 are also 
given the derivatives and, in addition, the component terms for comparison with 
the similar quantities of figure IV- 1. The solutions for this latter case are 
analogous in some respects to those of the Marrone and Treanor model (58) in 
that no account is made for the more rapid transitions occurring in the upper 
vibrational levels. The closer level spacing for the B-states is still retained, 
however. One sees that the temperature T^ is perturbed to values much lower 
than T^. This was also apparent in the quasi-steady solutions. The more 
interesting feature in the comparison of figures IV-1 and IV-7 is the very close 
similarity of the behavior of the separate terms contained in the derivatives. 

In particular, one may infer the importance of the Boltzmann factor in reducing 
the term Q^A in figures IV-lh and IV-4h. This term is relatively much 
smaller in the transient region for this case than in figure IV-1. One concludes 
that the effect of the Boltzmann factor is appreciable for the case of figure IV-1 
and hence greatly reduces the effect of the small relaxation time r A . The 
relative dominance of qJ'A and Q^A are interchanged in the quasi-steady 
region, but Qj'A i s still the major term. Thus, although anharmonicity has 
some effect on the solutions, it is not a large effect. 

(IV-B-2). The Induction Time 

In the preceding section the shock-wave structure was described in some 
detail, and the effects caused by changing the values of certain of the parameters 
were discussed. The section thus provided insight and familiarity with the 
model as regards normal- shock- wave flow. Here an additional comparison will 
be made with experimental data by using the experimental measurements of 
Wray (99) to check independently the evaluation of the parameters in the model 
equations. Wray measured the lag time before he observed molecular dis- 
sociation; this lag time is related to the elapsed time, after passage of a shock, 
during which the molecular oxygen concentration remains essentially constant 
(see, e.g., figs. IV-lc, 2c, etc.). 


To obtain his values Wray (99), using ultraviolet absorption techniques, 
measured the intensity of the 1470 A ultraviolet line during the passage of a 
shock wave. The absorption of this radiation, known as Schumann-Runge 
absorption (12, 99), is dependent on the number density of oxygen molecules 
in low-lying vibrational energy states. The radiation absorbed, measured after 
the shock wave passes some observational slit, is thus a measure of the 
molecular-oxygen concentration. The results reported by Wray were found 
by extrapolating back to the initial value the exponentially decreasing absorbed 
intensity, appropriately scaled. The shock-passage time associated with the 
intersection of the extrapolation and the initial value were defined by Wray as 
the "incubation” time. These results are denoted by At w and shown in figure 
IV-8. The smooth curves were obtained from the solutions of the model equa- 
tions. To obtain these theoretical results two types of induction times are 
defined. In one case the time, denoted At t , is the time for the concentration 
variable y ^ to decrease 0.1 percent. These results are plotted as curves (1) 
through (5). In the second case the straight-line portion of the oxygen- 
concentration curves is extrapolated back to the initial value, and the associated 
abscissa coordinate, denoted by Atg , is taken as the second induction time. 

The curves (6) through (8) correspond to this procedure. Each of the pairs of 
curves (1) and (6), (2) and (7), and (5) and (8) are obtained from the same solu- 
tion and are thus for identical values of the embedded parameters. 

Curves (1) through (5) illustrate the effect on At* of using different values 
for the embedded parameters. From curves 1 and 2 we observe first that 
changing the value of T u> the measure of how rapidly dissociation drops off 
for the lower vibrational levels in the preferential dissociation model (note 
definition in Appendix C), has no noticeable effect on Atj. The reason is 
evident if one examines the separate derivative terms associated with the dis- 
sociation processes as given in the previous section. The effects of dissocia- 
tion that depend on the parameter T^ are negligible within the time scale Atj . 
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Figure IV-8. Induction time computed using two different methods 
and compared with experimental data of Wray 
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A comparison of curves (1) and (4) shows the effect of changing the param- 
eter (see eq. 331-12). The change in C 2 , as was pointed out previously, 
corresponds to an increase by a factor of 2 in the vibrational relaxation time 
at the higher temperatures. The effect is observed here as a similar increase 
in the induction time at the higher temperatures and little difference at the 
lower temperatures. A comparison of curves (1) and (3) shows that an increase 
of 50 percent in the level separation parameter b results in a 50 percent in- 
crease in the induction time at the higher temperatures and an order of magni- 
tude increase at the lower temperatures. The result of using X = 1 rather 
than x - Xy is shown in curve (5). In this case there is only a 20-percent 
effect. 

The most interesting observation to make in figure IV- 8 is that with the 
specified values for the parameters, the group of curves Atj span values that 
nearly enclose Wray T s data. More importantly, the temperature dependence of 
the computed values also agrees favorably with Wray's experimental data. The 
choice of the induction time At x , associated with a 0.1 percent differential 
decrease in the concentration variable y ^ , was made for several reasons: 
first, it was desirable to have a method for determining an induction time At 
that would not require considerable plotting and that would serve as a natural 
point at which the numerical calculations could be terminated without the neces- 
sity of computing the entire shock profile. Second, it was desirable to have a 
scheme that is readily applicable even when appreciable dissociation is not 
evidenced. The agreement, however, is accidental since the values assigned to 
the embedded parameters were chosen somewhat arbitrarily. It was later 
realized that the induction time At A is not directly relatable to the experimen- 
tal time At w obtained by Wray? 


T These results may, however, be of value in an experiment that does not involve 
extrapolation, but where small percentage changes in the initial molecular 
concentration of the lower states are measured. 



The method for determining Atg is similar to that used by Marrone and 
Treanor (58) but is difficult to apply unless there is considerable dissociation. 
As was mentioned, the Shumann-Runge absorption involves primarily the low 
lying vibrational levels of molecular oxygen. By taking this appropriately into 
account a measurement of induction time At can be obtained that would be 
directly relatable to At^ and would lie between At t and At^ 11 ’ 12,13 ^. Be- 
cause the low lying levels are most significantly populated, however, it is to be 
expected that an appropriate evaluation of an induction time, for comparison 
with Wray T s data (99), would yield results that are close to the values Atg . 

A comparison of At 2 and At^ is therefore the correct comparison for an 
assessment of the model equations. 

Curves (6) through (8) show At 2 as calculated with different values for 
only a few of the embedded parameters. A comparison of curves (6) and (7) 
illustrates the effect on At 2 that results from changing the value of T^j. 
Although changing T-jj had little effect on At* , here the effect is appreciable 
(causing the curves to differ by nearly a factor of 2 at the higher temperatures). 
The principal reason for the different behavior here is that the procedure for 
determining At 2 , described earlier, also depends on the slope of the curve for 
y ^ and hence on the rate of molecular dissociation. Reducing the value of Tjj 
is equivalent to reducing the number of levels from which dissociation effec- 
tively occurs and hence reducing the rate of dissociation. The parameter At 2 
therefore not only measures the transient time but is also dependent on the rate 
of dissociation. A comparison of curves (6) and (8) illustrates the effect of 
anharmonicity on the induction time At 2 ; in this case the effect is appreciable. 
Curve (8) is included only to illustrate the effect of using x = 1 rather than 
X = Xjp- 

The comparisons that have been made illustrate the differences due to 
changing the values of certain of the embedded parameters. A comparison of 
curves (6) and (7) with the induction time At w of Wray, leads us to conclude 
that the model yields induction times that are about an order of magnitude too 
long. This is not too disturbing when one recalls that the initial values assigned 


to b were arbitrary and the values used for T^j were found so as to depend on 
b. In this connection it is worthwhile to review briefly the procedure that led to 
the values assigned to the embedded parameters (see preceding chapter). The 
value of the separation b was chosen to yield a separation energy that is 
roughly midway up the vibrational potential-energy well. Specifically, we have 
used E^/Ej) = 0.544 and 0.786. These values were chosen arbitrarily, although 
they were thought to be reasonable on the basis of the results reported by 
Treanor (94), Shuler (84), and Keck and Carrier (48). The values for the re- 
maining parameters, in particular T^ , were chosen (recall discussion in 
Chapter IQ) to provide a pre-exponential temperature dependence for the effec- 
tive dissociation rate that would be in reasonable agreement with experiment. 

It was also noted that the small values of the parameter T^ that were required 
for the proper pre- exponential temperature behavior, led to the conclusion that 
dissociation occurs only from the very uppermost vibrational energy level (see 
also discussion in Appendix C). In this respect the preferential dissociation 
model is consistent with other theories (7, 48, 73, 84, 85) since it is strictly 
"ladder climbing." In the other theories vibrational excitation occurs stepwise 
and dissociation (or recombination) involves the very uppermost level only. In 
summary, a choice of b at the midpoint of the vibrational-energy well requires 
highly preferential dissociation in order to produce an appropriate pre- 
exponential temperature dependence, and the induction time that then results 
is about an order of magnitude too large. An opposite result was obtained (see 
(58)) from a non-preferential dissociation model in that the resulting induction 
times were more than an order of magnitude too small. 

The parameters can be readjusted to yield more representative values for 
the oxygen molecule induction behind a normal-shock wave. Reducing the value 
of the separation parameter b will reduce the induction time. Such a change 
will yield a larger pre-exponential temperature factor for the effective dis- 
sociation rate constant (specifically, the vibrational coupling factor V) unless 
Tjj is also increased. Increasing the value of T u as b is reduced will main- 
tain the present agreement of the temperature dependence of the vibrational 



coupling factor and at the same time increase effectively the number of the 
uppermost vibrational energy levels from which dissociation occurs. One notes 
from figure IV- 8 that larger values of (compare curves (6) and (7) ) also 
reduce the induction time and would thus effect more favorable agreement with 
experiment. After the appropriate value is found for b and T^. , dissociation 
will likely involve principally just a few of the very uppermost vibrational 
energy levels, and the ladder-climbing concept will probably still apply. That 
the Marrone and Treanor model (58) yields induction times (when T^j = °°) 
that are too small while this model yields results that are too large provides 
confidence that the parameter values can be adjusted so as to obtain closer 
agreement with experimental data both as regards the temperature dependence 
of the effective dissociation rate and the induction time. 

Before concluding this section there is one additional comment that should 
be made. In Appendix B the procedure used to account for the anharmonic 
effects associated with the transition rates is discussed. The procedure in- 
volves the use of two linear segments to approximate the rate-of-quantum- 
transfer relation. This two-segment approximation, although it may be 
considered superior to a single-segment Landau-Teller approach, is perhaps 
too coarse. Introducing another segment into the approximation (and thus 
another separation parameter of the same type as "a") will also have the effect 
of reducing the induction time. Such an improvement may be of value although 
it is not expected to have a large effect on the induction time. 

IV-C. Nozzle Flow 

(TV-C-l) Example Solution 

In this section the flow equations, given by equations (F16) of Appendix F, 
are solved by means of the implicit numerical-integration method (eqs. (E22) ). 
The dimensionless streamtube area is obtained from A(x) = 1 + (x/3.20) 2 , 
which represents an axisymmetric hyperbolic nozzle. At large distances 
(x » 3.20) the geometry is approximately that of a cone with a half angle of 10°. 
The gas is taken to be a mixture of 20 percent molecular oxygen with a diluent 




of 80 percent argon. For this calculation the effects of considering other than 
argon as a colliding partner are again considered negligible. (The solution is 
thus obtained from the same equations as in the previous section except that 
the fraction of the mixture that is molecular oxygen is greater and the stream- 
tube area is not constant.) Only one solution will be given; this solution applies 
to nozzle flow when the reservoir has a temperature of 8,000° K and a pressure 
of 3,000 atmospheres. The starting conditions for the integration procedure are 
obtained by assuming that the flow is in equilibrium (isentropic) up to and 
slightly beyond the throat. The solution includes some of the effects of recom- 
bination that occur before the atom-recombination processes freeze, as well as 
effects of the ,f cooling" of both the upper and lower vibrational energy states as 
a result of vibration-translation energy interchange. The solution also demon- 
strates that the vibrational temperatures and the kinetic temperature are all 
equal for a relatively large distance downstream of the throat. It will be shown 
that this characteristic allows one to estimate an effective vibrational relaxa- 
tion time that is smaller than the relaxation time of Landau-Teller theory. Its 
value is quantitatively consistent with recent experimental data (98). It also 
turns out, however, that the smaller relaxation time (i.e., smaller as compared 
with the time obtained from Landau-Teller theory) is realized only when the 
vibrational temperatures are large. The nozzle solutions will be discussed in 
detail first, and later the procedure for estimating the characteristic vibra- 
tional relaxation time will be described. 

The results of the calculations are plotted in figure IV-9. Figures IV-9a 
through IV-9g show nearly all the dependent integration variables (T, t a- t b, 

7a > Tg, ^a’ P’ anC * to & e ther with the pressure p and the vibrational coupling 
factor V, which a**e found from the thermal equation of state (E3) and equation 
(CIO), respectively. Also plotted in these figures are values computed for the 
equivalent case of equilibrium nozzle flow, where the rate constants are all 
assumed to be infinite. These values are subscripted with the letter e. The 
derivatives and their component terms are plotted in figures IV-9h through 
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X, distance measured from throat, cm 

(b) Concentrations 
Figure IV-9 Continued 
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X, distance 


(c) Concentration, 
Figure 







IV- 91. The notation used in these figures is identical to that used previously 
(eqs. (IV-2) through (IV-7)). 

Beginning our discussion with figure IV-9a, we observe first that the tem- 
peratures T, T^, and are all equal to a distance of 7 centimeters down- 
stream of the throat. At this point the lower state temperature starts to 
differ from the other temperatures, and at 10 centimeters becomes constant 
(freezes). The upper-state temperature T^ and the kinetic temperature T 
continue to remain equal throughout the flow region shown. Other calculations 
(not shown) indicate that the upper-state temperature will also freeze when 
lower stagnation (reservoir) pressures are used. The freezing of T^ always 
occurs farther downstream than T^. We also note in figure IV-9a that T, T^, 
and T-g start to differ from their equivalent equilibrium-flow temperature at 
about 1 centimeter downstream of the throat. The flow is nonisentropic beyond 
this point. This is also evident from figures IV-9b and IV-9c, which show that 
the concentration variables y^, and y^ start to deviate from their 
equilibrium-flow values at about this point. At a distance of 3 centimeters we 
notice that recombination no longer occurs, as is evidenced by the constant 
atom concentration. The upper-state concentration variable continues to de- 
crease throughout the entire flow region, but since its value is so small there 
is no noticeable increase of the lower- state concentration. Changes of the 
lower- state concentration occur only where atom recombination is observed. 
That the decreasing y ^ affects the lower states and is not caused by dis- 
sociation will be apparent when the derivatives are investigated. 

The pressure, density, and velocity are compared with their equilibrium- 
flow equivalents in figures IV-9d through IV-9f. The pressure appears to be 
greatly affected by the chemical effects and not appreciably by the freezing of 
T^ , since the rapid freezing of the lower-state temperature T^ causes no 
further change of the pressure variable. Density, however, is affected by both 
the freezing of the chemical effects and of T^ . In particular, we have a 
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smaller flow velocity and larger mass density than is obtained on the assump- 
tion of complete equilibrium. (Notice that velocity is plotted with a linear 
scale while density is plotted with a logarithmic scale.) 

Figure IV-9g shows the vibrational coupling factor V computed for the 
nozzle problem. This factor differs at most by only about 2 percent from unity 
and for all practical purposes has little effect on the flow. This is in contrast 
to the situation for normal-shock-wave flow, where this quantity was very 
important. There the coupling factor delayed the onset of dissociation. For 
the quasi-steady zone (see, e.g., Chapter HI) this factor varied as T~^ and 
such behavior is also not exhibited here. 

An understanding of the processes that occur during the expansion in the 
nozzle is, perhaps, more easily obtained by examining the various terms con- 
tained in the model rate equations, as was done for the shock-wave problem. 
The values of these terms are displayed in figures IV-9h through IV- 91). The 
labeling used for the curves is defined by equations (IV-2) through (IV-7); the Q- 
quantities without subscripts are the derivatives while the quantities with sub- 
scripts are the components of the derivatives. One general remark should be 
made about the apparent unstable behavior shown by the complete derivatives. 
As was mentioned earlier, the implicit method was used for these calculations, 
and a characteristic of implicit methods is to show such behavior for the com- 
plete derivatives (notice that the separate terms are smooth). This effect does 
not appear in the solutions where, as may be seen, the curves in figures IV-9a 
through IV-9f are quite smooth. Although there appears to be a slight noise 
effect superimposed on the curve representing the vibrational coupling factor 
in figure IV-9g, the fact is that the scale is greatly expanded and the noise 
amounts to only about 0.02 percent. The error of the calculation is about this 
order of magnitude (see discussion in Appendix E). 

Figure IV-9h illustrates the behavior of Q^A . Near the throat the gain 
and loss terms are nearly equal and result in a derivative that is small, 
although sufficiently large to cause about a 30-percent change in 7^ . As 
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(j) Derivative, , and component terms, Qp* 
Figure IV-9 Continued 
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the temperature and density of the flow decrease the derivative decreases, 
and at distances downstream of the throat greater than about 3 centimeters, 
the derivative is so small that no further effect is observed for the value of 
the variable y ^ . This quantity is then said to be frozen. 


In figure IV-9i are shown the derivative and component terms for the 

y 

derivative Q a . Here we observe that the behavior near the throat is similar 
in that the gain and loss terms are nearly equal, yielding smaller net values 
for the derivative. The derivative Q^ a is negative, indicating that atom re- 
combination is occurring. This behavior continues downstream for about 2 

y 

centimeters, after which the gain term Q( a rapidly decreases in value. Atom 

y 

freezing occurs near the point where the value of Q ' a is approximately the 
same as that of its term Q7 a . The chemical "sudden-freeze" behavior, dis- 
cussed by Bray (14, 17) and resulting in his sudden-freeze approximation (14), 
is clearly evident here. The sudden-freeze, as normally applied, occurs where 
the dissociation term and the derivative have about the same values. 


In figure IV-9j are exhibited the terms for the derivative . Because 

the net effects of the processes near the throat are the result of the small dif- 
ferences between the terms q/ b and q/B ? and and q/B ? it is difficult 

to separate the dominant processes. One can only infer that all of the terms 
are important. At distances greater than about 3 centimeters the dissociation 
and recombination terms qJb an d qJb ? respectively, are small; and the con- 
tinued decrease observed for the value of y B (see fig. IV-9b) is due to the 
transitions y^-*- y^ . 

Figure IV-9k shows the behavior of the component terms for the vibrational- 
temperature derivative Q^A. The derivative oscillates erratically for this 
case, having both positive and negative values. The value at the peaks is about 
10 5 °K/cm, while T^ itself varies between 4,000° and 6,000° K. This general 
behavior also occurred in the normal-shock problem (see fig. IV-li), and for 
that case the comment was made that the upper states were in local thermo- 
dynamic equilibrium. The upper-state temperature derivative Q T B an d its 
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component terms are displayed in figure IV-91 and we observe that oscillatory 
behavior is again evident. In the present nozzle problem both the upper and 
lower vibrational energy states are in local thermodynamic equilibrium, con- 
sequently, Ta and Tg are equal to the kinetic temperature T (fig. IV- 9 a). 
When the variable freezes, at about 10 centimeters downstream of the 
throat, the associated derivative is small. 

Turning to the component terms in the derivative, we observe that the term 
Q^A , which represents energy addition to the lower states, is large. The 
vibrational relaxation terms Q^A and Q^A # however, although appearing to 
be of smaller magnitude, also have a large effect on the relaxation processes. 

It is also interesting to note that when T^ becomes constant, the terms Q^A 
and Q 3 T A are small relatively to the term Q^ A . We may conclude that the 
effect of anharmonicity is small when the dominant vibrational mode freezes, 
and the vibrational model yields results that are nearly the same as would be 
given by Landau-Teller theory. This point will be discussed further when the 
characteristic vibrational relaxation time associated with nozzle flows is 
introduced. 

For the upper states we observe that although the dissociation term 
and the atom-recombination term are large relative to the derivative 

Q T B , their net effect is small. The term , representing vibrational re- 

laxation in the upper states, is the dominant term as a result of the very small 

time scale associated with this equation (the associated equation is considered 

T 

as being in local thermodynamic equilibrium). The TT rate-limiting tr term, Q 2 B , 
appears to have only a small effect near the throat region, but it is very impor- 
tant after the lower-state temperature freezing occurs (about 7 cm downstream 
of the throat). On the basis of the calculation the significance of this term is 
difficult to assess. This term results from having different temperatures for 
the lower and upper groupings of energy states (see Chapter II); it may be re- 
lated to the similar rate-limiting terms discussed by Pritchard and Bray 
(14, 71). 







This example of nozzle flow has been introduced to illustrate the behavior 0* 

of the model equations for other than shock-wave flow. One important aspect of 
these solutions is that the upper vibrational energy states continue to relax even 
though atom recombination no longer occurs. If vibrational relaxation is gov- 
erned by a single term analogous to that of Landau-Teller theory (e.g. 3 as given in 
the Marrone-Treanor model (58) ) a different effect is observed: All the energy 
states freeze simultaneously. Bray (14) commented that a study of the complete 
set of vibration-dissociation population equations (master equations) in which 
harmonic-oscillator transition quantities are used for the rate relations, will 
yield the same result; that is, the vibrational -population distribution can be 
represented as a continuous sequence of Boltzmann distributions described by 
a single vibrational temperature, and the freezing of all the energy-level popula- 
tions occurs simultaneously. In his study (14) in which he allowed for multiple- 
level transitions for dissociations that occur above some upper vibrational 
energy level that he calls "rate-limiting/ 1 Bray obtained different behavior. 

The rate-limiting level or "bottle neck" freezes first, and subsequent vibra- 
tional relaxation is thereby prevented.. For this case he commented that 
appreciable energy may be frozen in the upper vibrational energy levels. This 
effect was not observed for the solution presented here. In this case the upper 
vibrational energy states continue to relax even after the freezing of the lower 
energy states. 

The reader will recall that in the previous section where the shock -wave 
flow was discussed, it was commented that the value used for the separation 
level b for the shock- wave solutions was probably too large. Reducing b 
might yield results analogous to the "bottle-neck" effect discussed by Bray. 

This speculation is based on the fact that the relative value of the rate limiting 
term Q 2 B is quite large for the case studied, particularly in the latter stages 
of the flow. As mentioned in the previous section, however, additional modifi- 
cation of the model equations may also be warranted in that more segments may 
be required for an approximation of the rate-of-quantum-transfer quantity 
introduced in Appendix B. It is not yet clear whether the bottle-neck effect will 
occur when the additional refinements are made. 
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(IV-C-2). Estimate of the Characteristic Vibrational Relaxation Time 
Applicable for Nozzle Flows 

One can estimate a characteristic vibrational relaxation time associated 
with the model equations for flows in which the upstream conditions are near 
local equilibrium. The procedure for obtaining this relaxation parameter is 
analogous to that described by Vincenti and Kruger (p. 236 in (97) ). The study 
depends on the results given in the previous section in that we observed that 
the vibrational temperatures are equal to the kinetic temperature throughout 
the greater portion of the flow and, once freezing begins, it takes place quite 
rapidly. This effect allows us to answer in a very simple manner the question 
as to what effect anharmonicity has on the process of vibrational freezing. 

We shall not include the dissociation and recombination terms in our inves- 
tigation, since the effects being considered are largely independent of these 
processes. We base our investigation on equation (II- 103a) of Chapter II and 
seek the characteristic vibrational relaxation time associated with this equation. 
Rewriting the equation, ignoring the last two terms, we have 
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We recall that the quantities measure the departure of the respective 
grouping of vibrational energy states from local equilibrium and are defined 
according to the relation 


6i(T,T v ) = [q“(T) - q”(T V )] -M T V) 


(IV- 11) 



Since we are concerned with only small departures from equilibrium, q“(T^) 
may be expanded about the local value of the temperature T. We thus obtain 
the relation 

«i(T,T v ) = - C~(T) AT Z^Ty) (IV- 12) 


where 


T y = T- AT 

and C^(T) = dq?°(T)/dT is the vibrational specific heat associated with an 
harmonic oscillator molecule having an infinite number of vibrational energy 
levels. (For a description of the notation the reader may refer to Chapter II or 
Appendix A.) The essential point is that for the flows having kinetic tempera- 
tures that are large compared with the vibrational temperature 0 i associated 
with the level spacings, the above quantity is approximated by k AT Xi(Ty) 
where k is Boltzmann T s constant. In this relation the dependence on i is only 
through the truncation factor Xi(Ty). Equation (IV- 10) may then be written 
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When equation (IV- 10) is displayed in this form, the characteristic vibrational 
relaxation time, denoted by Ty, is readily recognized. We have 
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The effect of anharmonicity insofar as the model equations are concerned may 
be easily inferred from this relation. Landau-Teller theory yields the vibra- 
tional relaxation time r A . A study of the quantity in braces (equivalent to a 
study of the ratio ^ / T y) thus allows one to examine the effect of anharmonicity 
in causing deviations from Landau- Teller theory. Such deviations have been 
experimentally observed (see, e.g. , 43, 44, 45, 78, 79, 98). 

In figure IV-10 are plotted the ratio ^ / T y and the contribution of each of 
the three separate terms to this ratio. The ratio is functionally dependent on 
the separation parameters a and b, and on the temperature-dependent param- 
eter X(T) (see list of equations at the end of Appendix B). The quantity A(T) 
and the ratio ^/r B are directly dependent on this parameter. The values of 
these quantities for the curves plotted are 9, 16, and X = X^ 9 respectively. 

One observes that anharmonicity causes a reduction by a factor of between 1/3 
and 1/4 of the characteristic vibrational time t as compared with the time 
obtained from Landau-Teller. This reduction (maximum curve value) occurs 
at relatively large vibrational temperatures, that is, at Ty ~ 8,000° K 
(or Ty/0^ w 3.6 where = 2234). At the lower temperatures, 

T ~ Oa » the model yields results equivalent to Landau-Teller theory (8,49). 

V Aj 

An investigation of the various terms illustrates the effect of the contribu- 
tions of the various groupings of energy states on the net effect of anharmonicity. 
The first term is associated with the lowest grouping of vibrational energy 
states, A 4 ; the second term is for the states, denoted A 2 , that are midway up 
in energy; and the last term is the contribution associated with the uppermost 
grouping of vibrational states, denoted B. A close investigation of the various 
terms shows that the population factors Q^(T)/Qy(T) have an important 
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Figure IV-10. Anharmonicity effects from model equations 
compared with Landau-Teller theory; charac- 
teristic vibrational relaxation time 
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influence on the magnitude of the contributions. Although this is not explicitly 

shown, it turns out that changing the value of b changes the relative magnitude 

of (T^Ay^ an d ( t a/ t y) h as no effect on their sum. Increasing the 

value of the anharmonicity parameter X (see Appendix B), however, increases 

the net effect of anharmonicity, and, accordingly, reducing x reduces this effect. 

Also, it may be shown that setting x = 1 results in a ratio r fr of unity, 

A* V 

indicating that the model is equivalent to Landau-Teller theory in this case. 

Reducing the value of the parameter a has the effect of making the parameter X 

effective over a greater number of vibrational energy states (see Appendix B) 

and thus increases the net effect of anharmonicity. On the whole, the effects of 

anharmonicity are only crudely approximated by the model equations, a point 

that is best exemplified by figures B-4 and B-5 of Appendix B. Introducing 

additional segments for an approximation of the "rate-of-quantum-transfer" 

relation would yield more accurate relations for the ratio r /t if as many 

A 1 v 

segments are introduced as there are energy levels, in principle a precise value 
for the quantity Ay could be obtained. (Such an approach may not be justi- 
fied, since the values of the transition rates k i that would appear in the 
equation are still not certain.) It is expected, however, that with the values of 
ki used in this paper more accurate computation of Ay will not yield 
results that are greatly different from those reported here. 
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CHAPTER V 



RESUME 

With the master equations as a basis, a simplified mathematical mode has 
been derived that is useful for studying the effects of vibration-dissociation 
coupling in fluid flows. The effects of changing the values of the various param- 
eters contained within the equations has been investigated. Results have been 
given for the quasi-steady zone behind a normal shock wave, for the structure 
of the complete shock-wave profile, and for the flow in a nozzle. It has been 
shown that the equations contain the minimum structure required to obtain 
agreement with existing experimental data. As more refined experimental 
measurements become available, additional complication of the model rate 
equations may be necessary. The formalism is sufficiently general that modi- 
fications can readily be accomplished. 

The most important results of this study are the recognition of the neces- 
sary modifications required to Landau-Teller theory to account for 

(1) The finite number of vibrational energy levels, 

(2) Molecular anharmonicity as regards both the energy-level 
spacing and the collisional transition rates, 

(3) Non-Boltzmann population distributions. 

To take account of (1) the Landau-Teller equation was modified by the addition 
of a "truncation factor." In regard to (2), the "collisional rate of quantum 
transfer" was approximated by linear functions. The results yield more rapid 
vibrational relaxation of the upper energy states and closer energy-level separa- 
tion in these states. To approximate (3) the distribution of energy was repre- 
sented in terms of two vibrational temperatures, one for the lower energy 
states and one for the upper states. 


The thermodynamic quantities that result (partition functions, energy, and 
specific heat) accurately approximate the equivalent quantities computed on the 
basis of spectroscopically determined data. The temperature dependence of the 
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pre- exponential factor in the effective dissociation rate is accurately repre- 
sented. The assignment of values to the parameters in the model equations to 
effect agreement with experimental values of the induction time behind a normal 
shock wave is by trial and error. This process was not carried to completion 
and as a result values of the induction time too large by a factor of 10 were 
computed with the preliminary choice of parameters. The procedure for the 
required readjustment of the parameters has been discussed. The present 
results show qualitative agreement of the characteristic vibrational relaxation 
time associated with nozzle flows (expanding flow) in that this time is about a 
factor of 1/4 of that given by Landau-Teller theory. The equations reduce to 
the form of Landau-Teller theory when the temperatures are not too large. 

Although further adjustments of the parameters are indicated, the present 
study shows that agreement with available experimental data can be achieved 
with a single assignment of values for the embedded parameters. These param- 
eters (a, b(Tjj) and X) depend only on the choice of species in a mixture. This is 
in contrast to simpler models where different values were required for the 
parameters in order to effect agreement with different experimental results. 

The fact that the model equations can readily be solved in conjunction with 
the equations of flow has been demonstrated. The numerical integration pro- 
cedure that was used is also discussed. 
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APPENDIX A 


THERMODYNAMIC QUANTITIES 

The thermodynamic quantities (partition functions, internal energies, 
specific heats, etc.) that are used to obtain the solutions described in the text 
are given in this section. Of these functions, those that pertain to the internal 
energy modes — translation, rotation, and electronic excitation — are listed 
for the reader’s convenience with little comment. These functions are com- 
monly used, and if supplementary information is desired, the reader can refer 
to other references (see, e.g., 9, 59, 97). The functions relating to the vibra- 
tional energy mode, however, are discussed more extensively, because of their 
bearing on the vibrational model that is the subject of this thesis, and because 
they are not commonly used. 

(A-l). General Relations 

In the list of relations that follow, the subscripts T, R, Ei, and V refer, 
respectively, to the translational, rotational, electronic, and vibrational parts 
of the indicated quantities. To abbreviate notation the integers 1, 2, 3, and 4, 
respectively, may be used in place of the above subscripts. The subscript j 
denotes the species for which the quantity is defined. The symbols appearing 
in the various relations are defined in the list of symbols. In the case of the 
relations referring to the vibrational mode, the additional subscripts A lf A 2 , 
and B indicate groupings of vibrational energy states to which the molecule 
belongs. In the list that follows the quantities are separated according to their 
type, the "complete" quantity of each type is listed first and the component 
parts are listed thereafter. 



Partition Functions 



Qj = nQy 


(Ala) 

Q T J = V ^T,j 





(Alb) 

Q - ( T V j_I 

“ y RJ/ 


(Ale) 

Q Eij = E exp y ' 

e E£,A,j \ 

T / 

(Aid) 

Q VJ = Q A t ,j(T) + Q A2)j (T) 

+ q bj (T) 

(Ale) 


(The last equation applies in the case of vibrational equilibrium only; see eqs. 
(A14) for a description of the formulas.) 


Internal Energies 


e j = Z e i,j 

i 

(A2a) 

-b„t* 8t ! 

(A2b) 

e T,J = f R oT 

(A2c) 

e R,j = R o T ( a j " x ) 

(A2d) 

g E.M,j / e E4J,j\ 

eE « =E °^ E<J T ) 

(A2e) 


e V,j -[ Q Ai,j( T ) q A lt j (T) + Q A 2) j + Q B,j (T) 


(A2f) 

(The last equation applies in the case of vibrational equilibrium only; see eqs. 
(A30), (A38), and (A39) for a description of the formulas.) 
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Enthalpy 


h = e j + R o T + h ° 



(A3) 


Specific Heats at Constant Volume 



(The last equation applies in the case of vibrational equilibrium only; see last 
section for description of the formulas.) 


Specific Heats at Constant Pressure 


<°p>j = 


'ahj\ 
i 8tJ 


( c v)j + Ro 


(A5a) 

(A5b) 
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Temperature Derivative of the Specific Heats 

(Note: Asterisk denotes that vibrational contribution is excluded) 



d * 

(c >>j 


dT 



dT ^T.j 

_d_, \ 

dT (c v) R>j 
d 


Rp 

T 


2(Cv) E^,j + ^ 




i Q Ei,j 


T / 6XP \ T 


"El,l,j 

R C T 


3(Cv) E^,j ( e E,j,j 

+ l RnT 


(A6a) 


(A6b) 

(A6c) 


(A6d) 


Entropy Function 


The contribution of each internal energy mode to the entropy can also be 
defined. Such relations, however, are not useful in this paper and will not be 
given here. The entropy associated with a specific species, say j, will be 
used and is 



1 + 


e j 

R 0 T 



(A7) 


where the quantity Qj denotes that the volume V (see eq. (Alb) ) is omitted 
from the translational partition function prior to the formation of the product 
denoted by equation (Ala). The variable Pj denotes the partial pressure 
associated with the species j (Pj = nijkT) . 


Equilibrium Constant 

A general form for a chemical equation (see, e.g., 97) is 

Ht v ij = Z/ ^ij i = 1 

j j 


(A8a) 
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where the and are the stoichiometric coefficients associated with the 
reactants and products, respectively, and Mj (chemical symbol) denotes the 
j th species. In terms of this notation the equilibrium constant /c^(T) , associated 
with the i th reaction, is: 

^(T) = 


where 

and h? is the heat of formation of the j th species. 



(A- 2), Derivation of the Relations for the Vibrational Mode 

Early researchers (see, e.g., 8, 30, 36, 62) assumed that a molecule could 
be described as a harmonic oscillator with an infinite number of vibrational 
energy levels. Such a description is attractive because of the simple thermo- 
dynamic relations that result (see, e.g., 97). These relations are reasonably 
valid provided that the temperature T assigned to the vibrational mode (the 
subscript V on T denoting vibrational temperature is dropped in this section 
since it only complicates the notation) is very small compared with the char- 
acteristic dissociation temperature 0 D (« E N _j/k). At temperatures com- 
parable to or higher than the characteristic dissociation temperature, such a 
description fails since the vibrational partition function and the vibrational 
energy increase with increasing temperature without limit. It follows, in 
particular, that at high temperatures the amount of energy contained in the 
vibrational mode exceeds the energy required for dissociation. 

Later researchers (see, e.g., 32, 82, 90) remedied this problem by trun- 
cating the number of vibrational levels so that the energy of the uppermost level 



and the dissociation energy were approximately equal. The thermodynamic 
relations that result for this "truncated harmonic oscillator" are also simple; 
however, the effect of truncating the number of vibrational energy levels is to 
over-compensate in the opposite sense at the higher temperatures (4, 58) in 
that the vibrational partition function and the energy function are too small. 

In this section a vibrational model is developed that includes a represen- 
tative number of vibrational energy levels and the effect of molecular anharmo- 
nicity, thus yielding more accurate values for the vibrational quantities. The 
model also retains some of the simplicity characteristic of the harmonic- 
oscillator models. The development evolves naturally from the approach taken 
in the text to study the effects of vibration-dissociation coupling (see Chapter 
II). The model still contains some of the shortcomings of other models (36, 58, 
93) in that the effects of vibration- rotation coupling are not considered. The 
lack of reliable values of the rate parameters pertinent to rotation precluded 
the addition of rotation-vibration relaxation effects into the rate equations (see 
Appendix B); consequently, it is questionable whether it is of value to include 
these effects in the thermodynamic quantities. For an explanation of some of 
the effects of vibration-rotation coupling, the reader is referred to the papers 
by Bauer and Tsang (4) and Bauer (5). 

The approach taken in developing the thermodynamic quantities for the 
vibrational energy mode is to assume that molecules can be split according to 
vibrational energy into two separate groups, a lower and an upper, in the same 
manner as described in Chapter II of the text. The separation level is denoted 
by a (see fig. A-l). The lower group of molecular energy states has relatively 
large energy spacing; in fact, it has the spacing commonly ascribed to the 
vibrational mode of harmonic-oscillator models (e.g., see 97). The upper 
group of energy levels has a reduced energy spacing. The number of vibra- 
tional energy levels is truncated so that the uppermost level has an energy 
corresponding to the dissociation limit. An additional separation b is intro- 
duced to make the thermodynamic quantities derived here identical with the 
quantities required in the text. As far as the derivations in this section are 
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concerned, this separation may be assumed at any arbitrary position within the 
upper group of vibrational levels. The set of energy levels in figure A-l illus- 
trates the placement of the levels for the model molecule. The levels are 
numbered consecutively starting with zero. The separation levels a and b 
are labeled, and the A t , A 2 , and B groupings of energy states are indicated at 
the left-hand margin. For a Morse oscillator the vibrational energy E v for 
the vth level is given by references 9, 40. 


where 


( e v) m = Hv[l - (xo) M v] v = 0,1 N m 

t v 

X °M 4D 0 

N M = 1/2 (x 0 ) M (A10a) 


The best values for E v , from spectroscopic data, are based on the relation 

(E v ) h = hjo) 0 v - 0) 0 (x 0 ) H v 2 + “o(y 0 ) H v 3 J v = 0,1 N h (AlOb) 

where the constants a; Q , co 0 (x 0 ) u ., and v 0 (y 0 ) are g iven > for example, by 

H ri 

Herzberg (40), and Npj is defined as the largest integer such that — D 0 . 

It is worthwhile to point out that the Morse-oscillator model is known to con- 
tain too many vibrational energy levels (see, e.g., 58), and the partition func- 
tions that result are therefore too large at the higher temperatures. 

The value commonly assigned as the vibrational energy-level spacing for 
harmonic-oscillator models is E^, defined by the relation^ 

e A a = ( Ei )h ~~ ( E o)h ~ ^oj^ - ( x o)mJ (AlOc) 

We will utilize this quantity as well as an additional level spacing E A , the 
upper level spacing, defined analogously by 

+ The differences in E Ai , depending on whether (x 0 ) M or (x Q ) H is the anhar- 

monicity coefficient used in equation (AlOc), are insignificant for most 
practical purposes. 
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Model 

molecule 


Anharmonic 

molecule 



Figure A-l. Diagram showing relative placement of vibrational energy levels for model as compared 

to anharmonic-oscillator molecules 



E A 2 - E B - ( E a) M - ( E a-l) M 



= Hco 0 [l - 2a(x Q ) M ] (AlOd) 

Here it is important to use the Morse oscillator harmonicity constant, ( x o) M 
that yields smaller energy-level spacing than would be obtained from (x 0 )^ in 
order to compensate partially for the fact that the levels are spaced even 
closer near the top of the potential well. Although this evaluation of 
results in only a qualitatively correct level spacing, it will be shown later to 
lead to good agreement of the thermodynamic quantities as compared with those 
obtained from equation (AlOb). 


The energies at the separation levels are given by 



E a = < a - DE Al + E A 2 

(Alla) 


E b = E a + ( b “ a ) E A 2 

(Allb) 


E n = E b + (N - b) E Az 

(Allc) 

where the constants a, b, and N are defined as the largest integers that satisfy 
the relations 


a - [< E a> M - (Ea, - E A ,)]/E Al 

(Alld) 


b ' a + N m - <*■>*]/** 

(Alle) 


N . b + [d q - (E b ) M ]/ EA! 

(Allf) 

The quantities ( E a) M > ( E b)]vr anc * D o are specified. In terms of these 
relations the energies for the various levels for the model are given by 



E v = E Aj v v ^ a - 1 

(A12a) 


= E a + E Aa (v - a) a^v^b-1 

(A12b) 


= Ejj + E a ^ (v-b) b =£ v < N - 1 

(A12c) 
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The important thing to observe is that the thermodynamic quantities for the 
vibrational degrees of freedom are defined in terms of three parameters, E^ , 
e A 2 ( or E a)> ancl N ( or E o)« T ^ e Morse oscillator contains two parameters, 
c o Q and (x G ) M (or D 0 ), and the truncated harmonic oscillator also contains two 
parameters E^ and Njjo(or D Q ). The constants u) Q and insure that the 
various models have proper low-temperature dependence, and the inclusion of 
D 0 , directly or indirectly, insures that the energy cutoff is accounted for. The 
model has one additional parameter, (or E a ), that lies between the lowest 

and highest levels of the potential well and improves the behavior at inter- 
mediate temperatures. 


(A-2a). Partial Vibrational Partition Functions 


The partition functions corresponding to the various grouping of energy 
levels can be evaluated in terms of parameters defined previously. The partial 
partition function is defined by the sum 


Qi(T) 



(A13) 


where and are the appropriate bounds corresponding to the grouping of 
energy levels indicated by the subscript i. The subscript 3 , denoting a 
specific molecular species, is omitted from what follows to simplify the 
notation. Substituting the appropriate values for the energy E v (eqs. (A12) ) 
into the preceding equation, we obtain 


a-l 


«A,f r > - S 

v =0 



1 - 


1 - 


exp 1 


exp 


(- 


aE A 

kT 

fAj 

kT 


0 

) 


(A14a) 


208 



(A14b) 




/ E a ^ b v ‘ f (v - a)E^"| 

^ ^ H- — w — J 

- exp (- h) q 'a i ( t > 

„ „ , / E b\ N v* r ( y - b > E A,i 

Q b (T) = exp(- — ) £ e*p[- S J 

v=b 

= exp (- w) q bW 


where the primed quantities are defined by 

(b - a) E a 


Q' A (T) - 


1 - exp 


[- 


kT 




1 ‘ exp (' 15 s ) 
(N- b )E Ai - 


Qg(T) = 


1 - exp 


[- 


kT 


1 - exp 


m 


(A14c) 


(A15) 


(A16) 


The primes attached to the thermodynamic quantities and Q B denote that 
the relevant partition functions are evaluated relative to their separate zeroth 
states, E a and E^, respectively, and not to the molecular ground state E Q . 


The fractional population of molecules that have vibrational energy corre- 
sponding to some specific level within a grouping of levels (say, for example, 
the vth level within the B grouping of levels) is represented by 


(-30 


exp 


B 


Qn(T) 


Q^(T) 


(A17) 


b — v — N - 1 


Furthermore, given the total vibrational partition function 


Q V (T) = Q Al (T) + Q Az (T) + Q B (T) (A18) 
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we can obtain the fractional number of molecules in any group as a ratio of the 
partial partition function and the total partition function. For example, the 
number of molecules in the A 2 -grouping of states relative to the total number 
of molecules is given by 


n A 2 QA 2 ( T ) 

n ~ Q V (T) 


These ratios are important, for example, when the initial values for the numer- 
ical integration are required. It is necessary to point out that the latter ratio 
is valid only for vibrational equilibrium, while the former ratio (eq. (A17) ) is 
always valid provided the appropriate vibrational temperature characteristic of 
the grouping of states (T 3 for the example cited) is assigned. 

In figure (A2) the molecular vibrational partition functions are compared 
for molecular oxygen, computed on the basis of four different models, plotted 
as a function of temperature. The curve labeled "Morse” was obtained from 
equation (A13) with values of E v computed from equation (AlOa). The curve 
labeled "Herzberg" was obtained in a similar manner, but with E v computed 
from equation (AlOb) and with values for the anharmonicity coefficients, x Q and 
y Q , obtained from reference 40. Equation (Ale) was used (with internal param- 
eters having values a = 9, b = 16, N = 32, 0 . = 2234, and © A = 1800) and the 
plotted result is labeled "Model." The evaluation of equation (A14a) (a = 26 and 
© Ai = 2234) yields the curve labeled "Truncated Harmonic Oscillator." 

The Morse oscillator contains too many energy levels which, furthermore, 
are spaced too close together near the top of the potential well. Consequently, 
its partition function is generally conceded (see, e.g., 58) to be too large. The 
truncated harmonic oscillator, having relatively large level spacing for the 
upper vibrational levels, has too few vibrational energy levels and a partition 
function that is therefore too small (see, e.g., 4). The model introduced in this 
paper is in reasonable agreement with the Herzberg result that excludes any 
effect of rotation (as used here). As pointed out earlier, it is perhaps un- 
warranted to come to any specific conclusion as to the validity of any one of the 
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models for flow calculations, since in all cases the effects of rotation have not 
been considered and these effects may be large, particularly at the higher 
temperatures (see, e.g., 4). 


(A- 2b). Partial Vibrational Energy Functions 


The average vibrational energy associated with molecules in the ith group- 
ing of energy states is defined by 


-1 

h 


Qi(T) = £ E v .^ 

v=&t. 


(A20) 


The quantity q^(T) is readily related to the partition functions defined in the 
previous section by the following procedure. With the substitution of the appro- 
priate ratio for n v /ru (e.g., see text preceding eq. (A17) ), the preceding 
equation can be written 


-1 _ 

h E 


Qi(T) = £ 

V=£. 


Vi 


exp 


(-fa) 

V kT / 


Qi(T) 


kT 2 BinQ^T) 
8T 


(A21) 


The above definition is consistent with that normally encountered in statistical 
mechanics (e.g., see 73) and is convenient in practice since it gives the vibra- 
tional energy directly by differentiation of the partition function. 


Substituting the expressions for the partition functions into the preceding 
equation yields 

aE Ai' 


%( T > = 


E A exp (-‘5 1 ) aE A 1 «p(— ep) 


1 la) 1 - exp (“ *r) 


(A22) 


q, (T) = E + q' (T) 


(A23) 



q B ( T ) ~ + q B (^) 


(A24) 



where 





(T) 


Ea 2 ex P 

(-— ) 
kT / 

(b - a) E Az exp 

(b - a)E A l 

kT J 

1 - exp 1 

f, E v 

V kT ' 

1 - exp 

(b - a)E Az ' 
kT 



E A2 e xpl 

('15?) (N-b)E A2 exp 

' (N - bjE^' 
kT 

1 - exp l 

e a 2 \ 

\ kT ) 1 6XP 

" (N - b) E a ” 

L kT J 



(A25) 


(A26) 


These expressions can be modified to obtain an alternate form that is also 
useful. The first terms in equations (A22), (A25), and (A26) represent the 
expressions that one would obtain for an infinite harmonic oscillator (harmonic 
oscillator containing an infinite number of levels). The second terms are 
important only at the higher temperatures and, in effect, limit the values of 
the first terms. This may be shown in the following manner, with equation 
(A22) as an example. Dividing the right-hand side by the first term and retaining 
it as a factor, we obtain^ 


q Al (T) = 


E Al exp 


1 - exp 



exp 


(a - 1 ) 
kT 


E A, 


“ ini 


/Qa^ 


(A27) 


^All quantities defined in this and the next section contain factors having the 
functional form 1 - exp (a). This form is convenient for the accurate evalua- 
tion of the thermodynamic quantities. Because it is not possible to evaluate 
1 - exp (cv) accurately for all possible arguments on an electronic computer 
having only a finite number of significant digits, the following relation is used 
when the arguments are small (a < 10 -3 ) 

1 - exp (a) * -ajl + f [l + f (l + f )]) 


When the right-hand quantity is expanded, the equation is seen to be simply a 
truncated Taylor T s-series expansion of the quantity on the left. 
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Because of their frequent occurrence in the text, it is worthwhile to assign 
symbols to the separate factors in this equation. We first define 




(A28) 


As mentioned previously, this quantity is the vibrational energy of a harmonic 
oscillator containing an infinite number of equally spaced energy levels with a 
spacing E^. The second factor in equation (A28) will be called the "truncation 
factor" for reasons that will be explained, and is symbolized as follows: 


•V T) 




(a - l)E Al 

kT 


- in 



(A29) 


In terms of these functions equation (A27) is then simply 


%< T) - «*<*> V T) 


(A30) 


Before proceeding it is worthwhile to discuss the above functions. In figure A- 3 
dimensionless equivalents of these functions (with internal parameters having 
values a = 9 and 0^ = 2234 representative of molecular oxygen) are plotted as 
functions of temperature. Also indicated on the abscissa scale are the dimen- 
sionless units T/0^. Dimensionless quantities are plotted so that all three of 
the quantities can be plotted conveniently on one graph. At low temperatures, 

T << , the infinite-harmonic-oscillator energy function increases nearly 

exponentially with temperature. At sufficiently high temperatures, T » 0^, 
the function can be approximated by q^(T)« kT (or q^(T)/kT ~ 1 as indicated 
in the figure). The function continues to increase without bound with increasing 
temperature. Because of this characteristic, nontruncated harmonic-oscillator 
models are limited in their usefulness to fluid flows that have relatively low 
temperatures, T < a0 Aj . 

As regards the truncation factor, at low temperatures this factor is unity 

(as indicated by the constant slope in the figure that is approximately 

Q, / 0 o . ) and remains nearly unity for all practical purposes until the second 
a— l 
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term becomes important, that is, when T ~ (a - 1)0^ = © a _ t (the dimensionless 
quantity increases linearly in this region because of the factor T). At higher 
temperatures, T > @ a _ l5 the truncation factor may be accurately approximated 
by 0 a _j/2T (from the figure we observe that X^fTJTA^.j The maxi- 
mum value of the product q A (T) = q°? (T) X A (T) at these extreme temperatures 
is given by k0 a _ 1 /2, a constant (this is implied in the figure as a result of the 
approximate value q A ^(T) /ke^-j 


One can obtain similar functions for each grouping of energy levels. There 
results 


£' Or) = 1 - 
■C b(t) = 1 - 




With the defining relations 


qT (t) 

a 2 


E a , exp 


1 - exp 


£ 5 ) 

(-%) 


(A31) 


(A3 2) 


(A33) 


q~( T ) = q“ (T) (A34) 

the primed values of the vibrational-energy functions take the form 

^CD = q Az (T) Z' A { T) (A35) 

Qg(T) = q“(T) Z' b (T) (A3 6) 


The corresponding quantities for the unprimed vibrational energy functions are 
given by 


V T) = "a 2 (t) £ a/ T ) < A37 > 

q B (T) = q*(T) X fi (T) 
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(A3 8) 



where 


X A = X + (A39) 

A2 * %(T) 

X B (T) = X^(T) < A40 > 

B 

When each grouping of vibrational energy states has the same vibrational 
temperature, the vibrational energy per molecule (including all groupings of 
vibrational energy states) can be written 

n Ai n A, n R 

V T ) = ~ V T) + ~ %( T) + ~ q B ( T ) 

= + «A ! ( T >5a ! < T > + %( T >"V T >] 

(A41a) 

or, in the case that only A* and A 2 states are to be combined, 

V T > - %< t > + V T) V t, 1 < A41b > 

Figure A-4 shows the results from equation (A41a) for molecular oxygen. Also 
included are the corresponding quantities computed from three other vibrational 
models. These quantities are obtained by differentiation (see eq. (A21) ) of the 
relations displayed in figure A-2. The curve labeled "Morse" has values that 
are large and the curve labeled "truncated harmonic oscillator" has values that 
are small relative to the curves labeled "Herzberg" and "Model." The general 
behavior of the various models is thus the same here as in figure A-2. The 
energy quantity agrees very closely (differing at most by about 2 percent) with 
the Herzberg result. 

To conclude this section it will be shown that, in the limit as the spacings 
E^ and E^ become equal, the total vibrational energy tends to a function 
recognized as that for a truncated harmonic oscillator with equal energy-level 
spacing. When the appropriate values of the product expressions for the energy 
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Figure A-4. Comparison of the vibrational energy function computed using 
different models 
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functions (noting qy = q Aj - q^ = qg when = E^) are substituted into 
equation (A41), we have 


q y (T) = q~(T) 


Qa/T) • c a 1 ( t ) + Qa 2 (t)-Ca 2 (t) + Qb(T)Xb(T) ~ 
Qy(T) 


(A42) 


The expression in brackets is the truncation factor for the function qy(T). Sub- 
stitution of the appropriate expressions for the truncation factor and simplifi- 
cation of the result yields 

r(N - l)E Al /Qv(T)Y1 

X V (T) = 1 - exp|~ — — - in(— — )\ (A43) 

This is identical to X A ^(T) with N substituted for the parameter a as one 
would expect. 


(A-2c), Par tia l Vibrational Speci f ic Heats 

The contribution of the vibrational internal energy mode to the specific 
heat at constant volume is found from the temperature derivative of the average 
vibrational energy. Consistent with the relations defined in the previous sec- 
tions, there will be three separate contributions, A t , A 2 , and B, depending on 
the specific grouping of vibrational- energy states considered. Additional 
quantities corresponding to the specific energy of a non-truncated harmonic 
oscillator will also be introduced and such quantities will be superscripted oo 
similarly to the analogous energy quantities. Truncation factors £ ^ will also 
be introduced. 

The temperature differentiations are readily found and the resulting 
specific-heat quantities are as follows 

c v (T) = c“ (T) £ (T) (A44a) 

v At v A t Ai 

c^CT) = c£ (T) £ a ^T) (A44b) 


t These functions also contain the functional forms 1 - exp (a); see footnote in 
previous section. 



(A44c) 



c Vb (T) = c^ 2 (T)£ b (T) 

= Cy B (T)£ B (T) 


where 


Ov” A i ( t ) . 


dqJ/T) Bo(-^)^p(-^) 


dT 


1 - exp 


(4)1 


l2 ~ Ro 


'^(T) 


LRoT J 


/%\ 

exp^~j 


(T) = 


£ Aj (T) = 1 - exp 


£ A,( t ) = 1 " exp 


£-d(T) = 1 - exp 


■SO 1 ) 

f q A (T) 1 
a 2 

2 

exp( 


dT " R ° 

_ RqT _ 


< T / 


(a - 1) 

T ~ A l 


®A- " 2 ■*“' 


^Qa/t)' 


(b - a - 1) 0 Kw 

T e A 2 2 in \b - a , 


_ fN - b - 1) 


0A 


/Qb(t)V 
2 H N^b ) 


(A44d) 


(A45a) 


(A45b) 


(A46a) 


(A46b) 


(A46c) 


We note that the primed and unprimed values of the specific heat are equal, that 
is, the specific heat is not affected by the ground-state reference energy. 

If consecutive pairs of energy groups, or all the energy states collectively, 
have the same temperature, they may be readily combined. For example, the 
specific heat of the combined Aj and A 2 states can be obtained by differentiating 
equation (A41b) with respect to temperature. There results 


°va< T > - R o{^ 


R 0 + \R 0 tJ 


Qa, 


% 


^A, 


-i 


Q a l R 0 \RoT; 


f<lA> 

,RoT, 


(A47a) 


If all the energy states have the same temperature, the total vibrational specific 
heat is given by 


220 




The behavior of the vibrational specific heat (assuming that the vibrational 
energy may be defined in terms of a single temperature, eq. (A47b)) is exhibited 
for molecular oxygen in figure A-5. This is a dimensionless plot of the temper- 
ature derivatives of the energy functions displayed in figure A-4. The curves 
for the Morse and Herzberg models are obtained from equations very similar 
to equation (Ale), the electronic contribution of the specific heat, except that 
for the vibrational model the degeneracy g is unity and in place of the elec- 
tronic excitation temperature 0^ we use the vibrational temperature E v /k, 
where as before E v is obtained from either equation (AlOa) or equation (AlOb). 
The Model results are found from equation (A4f) or equation (A47b). 

These quantities display the same relative characteristics observed in 
figures A-2 and A-4. The curves for the Morse and Herzberg results show a 
slight anomalous behavior (the slight indentation of the curves between the 
temperatures 4,000° and 7,000° K; also the model curve crosses the Herzberg 
curve at the temperature 7,000° K). The Morse and Herzberg quantities are 
computed as differences of terms obtained by summation. The anomalous 
behavior is attributed to a slight loss of numerical significance for these 
quantities resulting from the procedure used in their computation. Since they 
are not used for the computations given in the text, it was not felt worthwhile 
to make them more accurate. The model results, however, are accurately 
computed. 
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Dimensionless specific heat, C Vw (T)/R 







APPENDIX B 



VIBRATIONAL TRANSITION RATES 

In Chapter n values are required for the vibrational transition rates 
k v -i f v> or > more specifically, for the parameters k^ 0 , k^ 0 , x(T), and the level 
separation a. The latter quantities are those used for the approximate rate 
relations (11-45) through (n-48) developed in the text. Obtaining these quanti- 
ties requires the introduction of an approximation for a special quantity "the 
collisional rate of vibrational quantum transfer." The linear vibrational models 
that have had widespread utility, for example, the Landau-Teller (49) or Bethe- 
Teller vibrational relaxation model (8) described in Chapter I approximate this 
quantity, in effect, in terms of a temperature-dependent linear function of the 
vibrational-energy-level number v. Such a relation is unsatisfactory since it 
does not take into account the effect of the more rapid transitions that occur for 
vibrationally excited molecules. It will be shown that a more accurate relation 
can be obtained if at least 2 connected linear functions are used for the "rate- 
of-quantum-transfer" quantity. 

In the next section a discussion will be given of the transition rates appli- 
cable for anharmonic-oscillator models, and the results of several different 
theories will be compared. In the last section the procedure for approximating 
the transition-rate relations will be introduced. 

(B-l), Transition Rates, Analytical Form 

It will be shown that the vibrational transition rates, applicable for 
vibration- translation energy interchange involving atom-molecule collisions, 
may be represented in general by 

k v-i,v = vk o,i( T ) G ( T > v ; x o) (Bia) 

k v,v-i = exp ( E v 5 v- i/kT ) k v -i ( v (Bib) 

= vk l,o( T ) G(T,v;x 0 )exp|[E v>v - 1 - E 1>0 ]/kTj (Blc) 




where 



v = 1, 2, f N- 1 

E v,v-i = E v “ E v-i (Bid) 

and x G is the anharmonieity coefficient (xo = 0 for harmonic oscillators), v is 
the vibrational quantum number, E v is the energy of the vibrational level des- 
ignated v, k 0jl and k lj0 are the rates for transitions between the ground state 
and the first excited state of a molecule, and G is a factor that accounts for 
the effect of molecular anharmonieity on the values of the transition rates. The 
above formulations display the effect of anharmonieity explicitly through the 
function G. For example, setting G equal to unity (i.e., G(T, v; 0) = 1) yields 
the Landau-Teller transition rates (46, 49) for harmonic oscillators. 

One may speculate on the functional behavior of the anharmonieity factor 
G. Since the low-lying molecular states are nearly equally spaced and the 
effects of anharmonieity are therefore small, G should be approximately unity 
for values of v near unity. The energy-level spacing decreases for the higher 
vibrational energy states as a result of molecular anharmonic effects, and we 
may expect the rate of interchange of vibrational energy occurring during atom- 
molecule collisions to increase. We may therefore expect that G is a mono- 
tonically increasing function of v. In this section the G factor will be evalu- 
ated from the transition rates given by Keck and Carrier (48). It will be shown 
that the factor differs very little from a similar G-factor obtained from tran- 
sition rates derived by Nikitin (65) and used by Bray and Pratt (18). The rates 
given by Keck and Carrier are based on quantum- mechanical treatment of the 
collision processes, while the Nikitin formula, obtained from a perturbation 
treatment of the Landau-Teller results (49), is based on a semiclassical 
approach to the collision problem. In both cases the rates are described as 
being only qualitatively correct (see, e.g., 18); their main shortcoming is that 
they apply to nonrotating molecules; that is, they do not include the effects of 
vibration- rotation coupling. The effects of this mechanism may be important 
at the higher temperatures (see, e.g., 4, 5, 85), effectively increasing the rate 
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of energy interchange. Completely reliable values for the rate parameters 
describing the molecular kinetics, particularly the rotational processes, are 
not available. 

The rates given by Keck and Carrier (48) are for a nonrotating quantized 
Morse oscillator, for which the potential energy is given by (see 48 or 49) 

V M (r) = D e jl - exp[- /3(r - r Q )]j (B2) 

The vibrational transition rates are written 


_ . h. 

k v-i,v - *0,1 ai f ( Wl r) 


(B3) 


This formulation differs from reference 48 only in the notation and the fact 
that the factors independent of v in that reference have been lumped together 
here in the coefficient k 0jl . (The present representation is useful since the 
factor k 0?1 appears explicitly. This factor maybe found from experimental 
data as in Chapter m or its value can be estimated from the following 
relations.) The expressions given in equation (B3) are defined as follows: 


_ Z(m,n) f(o; 1 T) ai 
M “ mn 


(B3a) 


where 


Z (m,n) = dn >m c mn 
is the bimolecular collision rate, 

_ 2ju 3 2 Nv 
&v ~ ©Mi 2 (2N - v) 

is the resonance transition probability for a Morse oscillator, 

f(x) = x 2 f e~y csch 2 (xy _1 / 2 )dy 

j Q 


is the adiabaticity factor. 




kT (E v - 


(B3b) 


(B3c) 


(B3d) 


E v _ t )A = <Oo(l ~ §) 


(B3e) 



is the angular frequency of the transition, 

E v ~ © 0 v [l ~ = © 0 v[l - x Q v] (B3f) 

Hoj 0 

is the energy (ergs per particle) of the vibrational level designated v, 0 o = 
is the ground- state level spacing divided by kT, d^ m is the kinetic cross 
section for collisions between molecules n and collision partners m, and m 
and n are the concentrations of colliding partners and molecules. The con- 
stant |Ug is the reduced mass for the collision, (j, i2 reduced mass for 

the molecules, r= 7 tlV^ 3 /2 kT is the effective collision time for an exponential 
interaction potential V 3 = V D exp (- r 3 /L) , and u> 0 = /?V2D e/^ 12 is the angular 
frequency of the transition between the ground and first excited vibrational 
states. The quantity c = v8 kT / tt is the mean speed for collisions of mole- 
cules and collision partners, and N is approximately equal to the number of 
vibrational levels in the potential well. This notation has also been altered 
somewhat from that of Keck and Carrier (48) to conform with the notation of 
this thesis. For x » 1, the adiabaticity factor, equation (B3d), can be evalu- 
ated by the method of steepest descent and is found to be 

f(x) ~ fj(x) = 8 (~3') x?//3 ex P(“ 3x 2 / 3 ) (B4a) 

By numerical integration and curve fitting. Keck and Carrier (48) obtained the 
expression 

f(x) ~ f 2 (x) = ^[ 3 - ex p( _ x)J exp (“if) (B 4b ) 

which they point out is accurate (±20%) for 0 ^ x == 20, and bridges the gap 
between impulsive and adiabatic energy exchange. 

It will be shown here that equation (B4a) can be related to the rate relation 
found by Nikitin (65). On this basis, the G- function, defined by equation (B4a), 
can be written 



G,(T,v;x e ) , p[- aW - ^/»)] 

' + l( 1 -^i) + 3(‘*> T ) ,/, ( 1 - n) 7 

*n-m\ (B6) 


We obtain an approximation for G that is valid when v is small compared 

with the total number of levels, N, by the substitution of the approximations 

2/3 2 11 

in(l ±6) « ± a , (l ±6) » 1 ±5 6 , and — ~ 1 where 6 in this case is some 

small number 5 < 1. The G-function can then be represented by 


G(T ,v ; x D ) = exp { 2 ^ [(x 0 t) - l]J 

(B6a) 

= exp 1 4 x Q (v - 1) [(x 0 t ) 2/3 - l] J 

(B6b) 


where the subscript has been dropped in order to differentiate the expressions 
denoted by equations (B5) and (B6). The second of the two equations follows 
from the fact that for a Morse oscillator the number of energy levels is given 
by N = 1/2 x Q (see, e.g., 9, 40, or 48), that is, the anharmonicity coefficient 
x Q is determined from 


**>0 J_ 
X ° 4D 2N 


(B7) 


The factor co 0 t may be written in a more convenient form such that its tem- 
perature dependence is explicitly displayed. We first introduce a dimension- 
less constant defined by 


% 



(B8) 


In terms of this quantity we then obtain 



<B9) 
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The G- function may then be written 


G(T,v;x q ) = exp|4x 0 (v - l) 
= exp|a? F (v - 1)J 


<Po\ 


' ^0 

kT 


1/3 


where 


°T 


= 4x r 


<P 


f * u o 
°\ kT 


1/3 


(BIO) 


(BlOa) 


is introduced to further abbreviate the notation. 


The expression in equation (BIO) is relatively simple in outward appear- 
ance. When substituted into the equation for the transition rates, equation 
(Bla), it yields an expression very similar to the following relation obtained 
by Nikitin (65): 


^V-l,v) Nikltin = ex P 


4x 0 (v 



1/31 


vk M G Nikitin 


(Blla) 


where 


G Nikitin ~ ex P 


/ a w 0 \ 1//3 

4x 0 (v - i)<PoygfJ 


We observe that 

G Nikitin = G exp[4x 0 (v - l)] 


(Bllb) 


(B12) 


The relations thus differ only at the higher vibrational energy levels; at the 
uppermost level, where v = N (recall x Q = 1/2 N), they differ by a factor e 2 . 
Because these expressions are only qualitatively correct at these high vibra- 
tional levels, factors of order e 2 are certainly not unreasonable. 

It is interesting to compare the G- functions obtained from the adiabatic ity 
factors (eqs. (B4a) and (B4b) with the function of equation (BIO). This com- 
parison is shown in figure B-l, where the parameters have been evaluated for 
oxygen (see constants in Appendix G). The solid-line curves are plots of 
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Figure B-l. Plot showing relative deviation of the rate anharmonicity coeffi 
cients computed using different adiabaticity factors 


[in Gj - q: f (v - l)]/a-p (v - l) , and the dotted curves are plots of [in G 2 - Qf^x 
(v - l)]/a-p(v - l) ; in both figures the vibrational number v is plotted along 
the abscissa. The plots are made for three different temperatures that yield 
values of 1ico 0 /kT of 8, 1, and 0.02, and the corresponding curves are labeled 
1,2, and 3, respectively. In effect, the plots illustrate deviations of the G- 
functions from exponentials whose arguments vary linearly in v and inversely 
in the one-third power of the temperature. Such functional behavior is exhib- 
ited by horizontal lines in the figure. The solid-line curves, corresponding to 
the adiabaticity factor f t of equation (B4a), deviate significantly from such 
functional behavior, especially when the vibrational number v is large. Under 
these conditions the parameter oyr (see eq. (B3e) and discussion) is small 
(of order unity), and the quantum-mechanical and semiclassical derivations 
leading to these quantities are questionable. Keck and Carrier (48) introduced 
the adiabaticity factor f 2 of equation (B4b) to formulate more accurately the 
correct behavior of the transition rates for the smaller values of ay r. The 
dotted curves corresponding to this adiabaticity factor are nearly horizontal. 
These curves, however, exhibit a strong low- temperature dependence that is 
suspect. When v is small and the temperatures are lowest, the dotted curves 
correspond to percentage deviations from zero that are quite large (25 to 30 
percent), and we observe that the quantity ay r has its largest values. In this 
case we expect that the solid-line curves are perhaps more accurate. This 
supposition follows since at relatively low temperatures and for the smallest 
values of v the solid-line curves are equivalent to the Landau-Teller k 0>1 
transition rates for which there is good experimental agreement. Also, for 
this case, there is close agreement with the Nikitin formula (eq. (BH)). We 
shall consider that the G-factor to be used in equation (Bl) is given by equation 
(BIO) for all temperatures and for the complete spectrum of values of the 
vibrational quantum number v. 

(B-2), Approximating the Transition Rates 

From the relations of the previous section, one can see that a vibrational- 
relaxation model containing transition rates that vary linearly with the 


vibrational quantum number v cannot account for the rapid relaxation that 
occurs in the upper vibrational levels. Linearly varying rates, such as result 
from the study of collisions of atoms and harmonic-oscillator molecules 
(Landau-Teller (49) and Jackson and Mott (46) ), are very attractive, because 
the resulting relations greatly simplify the energy-moment equations derived 
from the master equations. This is exemplified by the analysis in Chapter II 
(see also, Bethe and Teller (8) and also p. 198 of Vincenti and Kruger (97)). 

To account properly for the nonlinear effects resulting from molecular anhar- 
monicity, one has the choice either of developing a relatively complex model 
based on nonlinear rate constants that are only qualitatively correct or of 
introducing approximations and thus extending the utility of the existing models. 
Such models have heretofore been remarkably successful in describing vibra- 
tional relaxation at the lower temperatures and have had considerable appeal 
in flow-field calculations because of their relative simplicity. A procedure 
will be introduced here whereby one can utilize the advantages of linearly 
varying rate relations, but still account qualitatively for the more rapid 
relaxation occurring in the upper states. 

The essential problem is one of correctly approximating Tj and T 2 in 
equations (11-49). This can be accomplished by appropriately approximating 
the quantities 

Ev,v - i ^v,v— l (B13a) 

E v+i,v ^v,v+i (B13b) 

by linear functions (note that E Y ^ 1 is defined by eq. (Bid) ). These relations 
represent the rate of transfer of vibrational quanta as a result of collisions 
(vibration- translation energy interchange). If one assumes a harmonic oscilla- 
tor then the energy differences E v V _ 1 and E v+1 v are constant and the transi- 
tion rates are given by vk i 0 and (v + l)k 0>1 , respectively. These relations 
have been basic in previous models. For a Morse oscillator the energy differ- 
ences between successive vibrational energy levels are not constant, but vary 
linearly with the vibrational number v and, furthermore, the transition rates 



depend exponentially on v (see eq. (BIO) ). A single linear relation is there- 
fore not satisfactory. We shall therefore approximate the above relations by 
two linear segments, as exemplified in figure B-2. 



Figure B-2. Approximation scheme for including anharmonic effects 

The ordinate is the product under consideration made dimensionless by the 
factor (E^ 0 ) _I and the abscissa is the vibrational energy E v ; the smooth 

curve I is the function to be approximated. We observe that this curve attains 
a maximum very near the dissociation energy D 0 . This maximum is a result 
of the Morse oscillator equations and is probably not a real effect. The curve 

I will here be approximated by the two straight-line segments II and in that 
intersect near the separation E a . The plot (and the approximations) is found 
with the vibrational energy E v as the independent variable rather than the 
vibrational number v for two reasons: First, the variation of the quantity 
E v v-i kv,v-i i s least in this case; and second, the thermodynamic energy 
quantities q^(T^) (see Appendix A) are more accurately included. The segment 

II appropriate to the lower states (with energy less than E a ) is the vibrational- 
relaxation function obtained from Landau-Teller theory (8, 49). At relatively 
low temperatures only these states are appreciably populated. Since existing 
low-temperature theories have been successful in their description of vibrational 
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relaxation, particularly behind shock waves, there is no reason to alter this W 

* 

function. The segment HE applicable to the upper states requires a different 
approach. This segment lies between a point on the smooth curve (with a 
coordinate E v ) that lies near the curve maximum and the intersection of the 
segment I at the separation E a . The value of the curve maximum depends on 
the temperature T. The value (a constant) of E a is chosen somewhat arbi- 
trarily and will be discussed later. The procedure is described in detail later, 
and compared wife others by means of figures at the end of this section. 



From the above discussion the required relations can readily be found: 
( E v,v-i k v,v-l) prox . = vk i,o E A X 0 - v =£ a (B14a) 


'approx. 

^approx. 


( E v,v-1 k V,V-l) aDDrox . - k l,o 


E 


(v - a)x R (T) + a "e^ j 


a < v < N 

(B14b) 


where 


X P (T) = 


c (2) 

^>0 

A 1 ) 

C l,0 


k v ,v- l j 

' ,v_1 k l 0 J 


- aE 


'Ai 


E, 


/m R 


- E f 


(B14c) 


is the slope of curve IE shown in the sketch and the parameter a is an integer 
evaluated in such a manner that (a - l) E^ + E^ = E a (see Appendix A) . The 
quantity v m depends on the temperature; its value is obtained as an approxi- 
mation to the actual coordinate v max at which Ey^j ky^y.! is a maximum. 
The procedure will be discussed shortly. 


The above relations are associated with rates of transfer of vibrational 
energy for transitions in the "backward" direction, that is, from the level v 
to v-1. Equivalent expressions for transitions in the "forward” direction, 
E v,v-i k v-i,v> are found by requiring that the model equations have the proper 
behavior at equilibrium. At equilibrium we have n v k VjV _ 1 = n w _ 1 k v _ ljV which 
follows as a consequence of "detailed balancing" or "microscopic reversibility" 
(e.g., see 97, 100). For the model equations, there results 
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k V-l,V / n v \ 

k v,v-i \ v_1 /equilibrium 

(B15a) 


( N 

= exp l- — J 

1 < v < a 

(B15b) 


- exp (- -4 s ) 

a < v ^ N 

(B15c) 

The required rate parameters, of the form given by equations (11-45) in the 

text, are obtained by dividing equation (B14a) by E Al and equation (B14b) by 

E A2 . We obtain 

k V,V-l - vk l,0 

1 < v ^ a 

(B16a) 


k v,v— i = (v - a) k$ 

+ 4; a _! a ^ v ^ N 

(B16b) 

where 

k i,o = k i,o 


(B16c) 


k M = k M X r (T) 


(B16d) 


, (2) ,(i) Ea i 

k a,a-! = k i,o a i7 
A 2 


(B16e) 


The reverse rates (eqs. (B15) ) are unchanged. 

To obtain the coordinate v m required in the above rate relations we need 
the derivative with respect to the energy E v . With Morse-oscillator values 
for the energy spacing (Ey^j =Ha; 0 [l - x 0 (2v - l)]) and equations (Bl) and 
(BIO) for the transition rates, the energy derivative of the product of interest 
is given by 


9E~ ^ n [ E v,v-i k v,v-tl T=i 


T=const. 


2x^ 


1 - x 0 (2v - 1) 


+ i + VT) - - -° ° 


KT 


Hw 0 [l - 2x 0 v] 


(B17) 
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To simplify notation we define 


21ia>nX n 

a B (T) = a F (T) jgr- (B18a) 

f /UofeV 1/3 I 

" 4 *l«>\ h 7 " 2H “ X J (B18b) 

The quadratic equation in the variable v max that results when the above 
derivative is set equal to zero is easily solved. The solution is given by 

v max( T ) = "2 

(B19) 


2 1 1\ //_ _2_ 31 If 4 ( X + x o) 

2x 0 2 y yy Oij^ 2 xq 2/ 2 Xq of 


The root with the negative radical is discarded because negative values of 
v max are meanin g less ' 

Since we are seeking only approximate relations, it is not worthwhile to 
use the complicated expression given above. A simpler one is obtained as 
follows. Neglecting the terms 1/2 as compared with the terms l/x Q and l/f3 
and x Q as compared with unity, we have 



The ratio 4x 0 /of^ is always less than 1 for the temperatures of interest. We 
can approximate the radical using the first two terms of a binomial expansion 
and obtain for v max the approximate quantity 


%/ T ) 





(B21) 


An alternative, less accurate formulation is obtained if terms of order 
(2x 0 /qi r ) 2 are neglected: 



2x o 

“h(T) 


1 


(B22) 




The quantities v max , v mi> and v m2 are plotted versus temperature in 
figure B-3. We see that v mi is always greater than v max and differs most 
at a temperature of about 24,000° K, where the error is about one vibrational 
level (or about 3 percent). The approximate quantity v m is always less than 
v max and differs at most from the correct value, v m£UC , by about five vibra- 
tional levels (about 13 percent error). The product E Y , Y _ X k V)V _ 1 actually in- 
creases relatively slowly with increasing v until it reaches its maximum and 
then decreases very rapidly (e.g., see fig. B-2). The most accurate value, 
v mi , yields values for k VjV -l in the rapidly decreasing region (v m > 

v max) an d the least accurate value, Vm 2 , yields values for the product quan- 
tity in the slowly increasing region (v m2 < v max ) . Consequently, whether one 
uses v mi or v m2 to represent v max actually makes little difference in the 
value of the quantity E V)V -j[ ky^y-i that is of interest. We shall therefore use 
the simpler relation v m2 when values of are required. The quantities 

( E v,v-i k v,v-i)/( E i,o k i,o)»( E v+l,v k v,v+i)/(Ei,o k i,o) and their approximated equiva- 
lents are plotted versus dimensionless vibrational energy, E v /E^, in figures 
B-4. The product relations appropriate for the M forward M transitions are 
labeled (f) and the "reverse" quantities are labeled (r). Figure B-4a is for a 
temperature ratio T/©^ of unity (9^ = 2234° K) , and figures B-4b, B-4c, and 
B-4d are for ratios of 2, 5, and 10, respectively. 


Discussing figure B-4a first, one notes that in this case the approximating 
method is extremely poor. It is obvious that the approximate rates have values 
that are too large in the region adjacent to and above the separation point 

(E v /E * ~ 9) and the forward rate is not large enough at the higher levels. It 

''■*1 

may be shown, however, (e.g., see Chapter HI and IV) that for normal-shock 
waves the upper levels are not appreciably populated at these temperatures and, 
therefore, the large differences between the approximate and Morse relations 
have little effect in the application of the approximate rates. At successively 
higher temperatures, illustrated by figures B-4b, B-4c, and B-4d (note that the 
ordinate scale changes in these figures), the approximate relations improve as 
regards their difference from the Morse relations, and the difference between 
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Figure B-3. Comparison of the vibrational energy level computed by different 
methods that yields a maximum "rate of vibrational quantum 
transfer" 
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Figure B-4. Comparison of the "rate-of-quantum-transfer quantity," (x R ) 
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Figure B-4 Continued 
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Figure B-4 Continued 
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(d) t = 10 e Aj 

Figure B-4 Continued 
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the forward rate relations also becomes less* For nozzle flow where the tem- 
peratures are high, the rates are therefore accurately approximated* Freezing 
of the vibrational mode generally occurs before the lower temperatures are 
reached. 


The value of the separation E a , at which the lower linear segment inter- 
sects the upper segment, was chosen so that the value of the approximate rates 
found for the lower segment, differ only by a factor of about 2 from the Morse 
values. This choice is somewhat arbitrary, and E a may be determined on the 
basis of other criteria when such approximations are worthwhile, for example, 
when more reliable values of the quantities being approximated are available. 
Then it may also be advisable to make the separation E a temperature depen- 
dent. It may also be worthwhile to introduce additional segments for these 
approximations (see Chapters III and IV). 


To assess the effect of the approximating procedure for the upper states, 
it is convenient to have a different set of relations from those just described. 
Such relations are obtained by following a procedure similar to that used to 
obtain equations (B14) through (B22). Bather than being based on the reverse 
rates, however, the approximating procedure is based here on the forward 
rates. It turns out that the resulting equations are very similar, except that in 
equations (B14b) and (B16d) X ^ is replaced by X F , where 


E v,v-i k v-i,v 

/% - %\ 

k 0,l 

expl T J 

v=v m _ 


aE 


Ai 


Xp(T) = 


'm F 


- E e 


(B23a) 
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The results from replacing X-^ by Xp is illustrated in figures B-5. The 
labeling has the same meaning as in the previous figures. Here we observe 
that the curves labeled (f )^ 0( j e ^ intersect near the maximum of the curves 
(f >Morse where in the Previous four figures the curves labeled ( r 
intersected near the maximum of those labeled ( r ) Morse - Temperature 
dependence is similar in that the behavior at low temperatures (see fig. B-5a) 
results in a relatively poor approximation. We also note that the differences 
between the curves subscripted "Model M and those subscripted "Morse" be- 
come less at the higher temperatures, and the approximated rate relations are 
more accurate. 

The quantities X R and Xp are plotted versus temperature in figure B-6, 
and their differences may be observed by a comparison of these curves. We 
see that the differences are greatest at the lower temperatures (about a factor 
of 4 at about 1000° K) and become less at the higher temperatures. At about 
24,000° K, X R and Xp are approximately equal. It is worth pointing out again 
that the X quantities, defined as the ratio k^g/kj|g (see eq. (B16d)), are also 
a measure of the effect of molecular anharmonicity on the rate processes. For 
a harmonic oscillator we have Xp = X^ = 1, since in this case both E v ^ v _ 1 
k v _ ljV and E v ,v-i k v,v-i are linear functions of v. It then follows that for 
those temperatures where the X(T) are large, the effect of molecular anhar- 
monicity is expected to be greatest; conversely, when the X(T) are small, the 
effects of anharmonicity are least. We thus observe in figure B-6 that 
the effects of molecular anharmonicity are greatest at the lower tempera- 
tures; although they still exist at the higher temperatures, they are greatly 
reduced. This result is not unexpected if one considers that when kT»E v>v _ 1 
(regardless of v) then the effect of the smaller level spacing is essentially 
lost in the collision process. This should not imply that the rate of energy 
transfer is less at higher temperatures. The quantity X(T) does not include 
the effect of k|*g (see, e.g., eq. B14a); k^g increases with increasing 
temperature. 
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Figure B-5 Continued 









To conclude this section, the rate parameters defined in Chapter II by 
equations (11-39), (11-56), (13-68), and (11-73) will be rewritten in terms of the 
quantities derived here. We obtain 
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where 



(B28a) 



where 



(B29a) 



(B29b) 


The temperature-dependent parameters a-p and a R are defined by equations 
(BlOa) and (B18b). 
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APPENDIX C 



PREFERENTIAL DISSOCIATION MODEL OF MARRONE AND TREANOR 


Marrone and Treanor (58) introduced a model that yields the vibrational- 
level dissociation rates in a useful and rational manner. Their model, referred 
to as the "preferential dissociation model", is relatively simple and introduces 
few parameters into the problem of vibration-dissociation coupling. Further- 
more, it is consistent with recent theoretical and experimental evidence that 
the dissociation probability is higher for higher vibrational levels. The deriva- 
tion given here is slightly different from that in the original paper, the most 
significant difference being that here the level dissociation rates N are 
assumed to be negligibly small for all vibrational levels below some level 
designated b. No such level was introduced in the original model. The final 
equations remain very similar to those of Marrone and Treanor. 

In the derivation a function p y is introduced that represents the fraction 
of dissociation events taking place from the vt h level. Thus the sum of p v over 
all levels is equal to 1. The function p y is given by 

Vn’v™ 

" ( Cl > 

k^m 


N-i 

= Z kyNiyn (C2) 

* 


where 


t'V 


k F n B m 


From collision theory (e.g., see p. 215 in 97), one obtains 


Z(n v ,m) M(D-E V ) f(n v ,m) 
Z(n B ,m) M(D-Efj) f(n B ,m) 


(C3) 
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where the various quantities are defined as follows: 



Z(m,n) 

M(e) 


bimolecular collision rate associated with a mixture of gas 


particles m and n ^n m n n d ^ m ^ ; 


activation or energy factor, that is, the fractional number of 
collisions with line-of- center kinetic energy greater than 
the activation energy e^exp^--^-^ ; 


f(n,m) steric factor, that is, the fraction of sufficiently energetic 

collisions between the gas particles m and n that actually 

result in a dissociation; N- 1 

C normalizing constant evaluated such that 7] p v = 1 . 

v=b 


If we assume that the collision diameter of a molecule is unaffected by its vibra- 
tional state , then we have Zfn^mJ/Zfn-gjin) = n y /n B . Furthermore, if the 
molecular population in the vibrational energy states is described by a 
Boltzmann distribution, the population ratio n^/n^ is given by exp(-E v /kT B )/ 
Qb(Tb), where -Qb(T b ) is the vibrational partition function associated with the 
upper states evaluated at the vibrational temperature, , associated with 
these states (see Appendix A). The activation factor (e.g., see p. 218 in (97)) 
is given by exp(-e/kT) provided only the kinetic energy of the relative motion 
along the line of centers between two colliding particles is active in promoting 
a reaction requiring an activation energy e. The key to Marrone and Treanor T s 
preferential dissociation model (58) is the hypothesis that the ratio of steric 
factors in equation (C3) increases exponentially with increasing vibrational 
quantum number; that is, 


f(n v ,m) __ /_ Dp-Ey 

f(n B ,m) (, 


(C4) 


where T^ (with dimensions of temperature), describes how rapidly the dis- 
sociation probability drops off for low values of v . 
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With the above relations substituted into equation (C3), we obtain an alter- 
native representation for p y as follows: 


Pv = C 


r m_l i jj\ 

exp L k \t b T~Tul 
Qb ( t b) ex P ( E b/ kT > 



(C5) 


The normalizing constant C is readily evaluated as 


C = exp (E b AT)exp(D 0 /k Tu ) 

and the quantity Tp is defined by 



(C6) 


(C7) 


The final form for p y is obtained by the substitution of C into equation (C5) to 
obtain 


Pv = 



Q b ( t f> 


(C8) 


This expression is very similar to that obtained by Marrone and Treanor. At 
equilibrium we note that T-p = -T^j and (P v ) e q = exp(E v /kT B )/QB(“Tu). In the 
case that T^ - °o, we have p v = exp(-E v /kT m )/Q B (T m ) where l/T m = l/T B - l/T. 
This corresponds to nonpreferential dissociation (equal probability that dis- 
sociation occurs from any B-state vibrational level; see, e.g., 58 or 93). For 
the other extreme case where T U <T, T B , we have p y = exp(E v /kT u )/ 

Qb(“Tu). We may approximate Qb(~ t u) exp(Ej/kT[j) exp[(N-b-l)EA 2 /kTu] 

(see Appendix A), and after introducing the appropriate relations for E^ and 
E v in terms of 0 a 2 from Appendix A, we have p v » exp [(v-N +1)0 a 2 /Txj] 
where N-l — v ^b. In the limit that T-q is small compared to 0^ we note 
that p v is negligibly small for all v except v = N-l. This case corresponds 
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to highly preferential dissociation— in essence, dissociation involving only the 
very uppermost bound vibrational state N-l. 

Substitution of the probability factor p v into equation (Cl) yields an expres- 
sion for the level dissociation probability k y N . This expression, however, is 
not yet useful since it introduces still another parameter into the problem, that 
is, the effective dissociation rate k F . We can relate kp to the quantity 
k;p e q.(T), the dissociation rate constant that would exist with vibrational equi- 
librium at the local translational temperature T. By definition the level dis- 
sociation rate, k v ,N> depends only on the translational temperature and, hence, 
is independent of the number of molecules in any of the B grouping of vibra- 
tional levels, that is, the vibrational temperature T-g. It follows then from 
equation (Cl), that the level transition probability may be defined by either of 
the following pair of equations: 



~ k-p P v 


Qb( t b) 



(C9b) 


Equating these relations, we obtain 


V(T b ,T) 



Q B ( t > (Pv) eq 
Qb( t b) Pv 


exp 




(ClOa) 


Qb(T)Q B (T F ) 

Qb( t b)Qb(-Tu) 


(ClOb) 


The factor V(Tg,T) thus defined is called the "vibrational coupling factor." At 
equilibrium, where T;g = T, we have V(T,T) = 1. The parameter kp e q (T) 
is, as already mentioned, the equilibrium dissociation rate. 



The level recombination rate expression k^ jV (T) is obtained by invoking 
"detailed balancing, 1 ’ that is, by assuming that k^ jV (T) is not affected by the 
degree of vibrational nonequilibrium. The appropriate expression can be found 
by requiring that the rate at which vibrationally excited molecules are formed 
as a result of atom recombination (and having energy corresponding to the vth 
quantum state) be identically equal to the rate that such molecules are dis- 
sociated (at equilibrium). There results 

^N,v( t ) = (kv,N)eq (^f) (Clla) 

\ a / eq 

■A. 

= (ky,N) eq/ n v\ _ k F( T ) e q , . 

“ ^T) \*b) eq ' k B (T) 

(Cllb) 

where k b (T) is the "partial" thermodynamic equilibrium constant associated 
with the reaction and is defined by 


*b( T ) 



(C12) 


The "partial" equilibrium constant k b (T) is not the same as the thermodynamic 
equilibrium constant normally encountered in studies of diatomic molecular dis- 
sociation, but differs by the factor n/n-g = Q(T)/Qg(T); that is, 


^b< t > 


Q V (T) - 

(C13a) 

_ Qv< T ) - / E b\ 

Qg(T) K * T * eXp VkT/ 

(C13b) 


where Qy(T) is the complete vibrational partition function associated with the 
diatomic molecular species and k (T) is the normally encountered equilibrium 
constant (e.g., see eq. (A6)). 



To recapitulate, the form of the equations representing the level dissociation 
and recombination rates will be set down. These equations, obtained by appro- 
priate substitution of equations (C9) and (CIO), are as follows: 



= k 


nB 

F Pv ny 


(C14a) 
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where p v is defined by equation (C8). 
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APPENDIX D 


A STUDY OF THE EFFECT OF THE TRUNCATION FACTOR: PROBLEM OF 

UNCOUPLED VIBRATIONAL RELAXATION 

The simplicity of the Landau-Teller vibrational relaxation equation as 
obtained by Bethe and Teller (see Chapter I and also (8)) is noteworthy. Fur- 
thermore, as a result of the theoretical investigations of Rubin, Montroll, and 
Shuler (62, 75, 76, 81) and the favorable agreement with experiment the equa- 
tion has received widespread acceptance. Some investigators, however, (e.g., 

57, 58) have altered the original formulation by the use of accurate thermody- 
namic vibrational -energy quantities with the aim of generalizing the equation 
to account (1) for the effects of anharmonicity on the vibrational energy and 
(2) for the fact that there are only a finite number of vibrational energy levels. 
Such generalizations have been heuristic in that no attempt was made to justify 
their use in a formal manner on the basis of the master equations. Such "im- 
provements" are questionable and, furthermore, may lead to erroneous conclu- 
sions (57) concerning the effects of anharmonicity. 

To exemplify this the differences obtained by the inclusion of the truncation 
factor in the vibrational relaxation terms will be investigated for the problem of 
uncoupled vibrational relaxation. The effects to be discussed are somewhat 
hypothetical in that when the differences are sufficiently large the effects of 
dissociation would also be important and should not really be ignored. We are 
concerning ourselves then only with the study of the vibrational relaxation terms. 

We shall compare vibrational relaxation as given by the previously used 
equation 


dqy (Ty) q v (T) - qy (Ty) 


(Dl) 


dt 


t(T) 



with that given by 


d<ly (Ty) 
dt 



£y(T v ) 


%( T ) 


'V 


(T v ) 


JCv(T) 


~ q y ( T y) 


t(T) 


(D2a) 


which is the relaxation equation obtained in Chapter II, except that the effects of 
anharmonicity are ignored and only the truncation factor is included (see 
eq. (11-106)). In this equation the function qy is that given by equation (A 42). 
Its value of N = 26 yields a characteristic dissociation temperature 0p that 
is approximately correct for molecular oxygen (0p = 59,368° K). It must be 
emphasized that the effect of anharmonicity is not being taken into account in 
this investigation. The comparisons that will be made would show even larger 
differences if such effects were accounted for in the manner discussed in 
Chapter II. 

Equation (Dl) can be readily integrated to obtain 



( T v) ~ ^ ( T ) 

qy(To) - Qy (t) 


(D3) 


where T 0 is the initial temperature (i.e., the temperature at t/r = 0). We can 
also solve equation (D2a) formally by quadrature as follows: 



c v (x) dx 


Z v( x ) _ 

£y(T) 


q v (x) 


(D4) 


This integral can be evaluated numerically. 
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The solutions are compared in figures D-l and D-2. Here the abscissa is 
the dimensionless time t/r and the ordinate is the dimensionless vibrational 
temperature Ty/T. The dashed -line curves are obtained from equation (D3) 
and the solid-line curves from equation (D4) with the integration done numeri- 
cally using an electronic computer. We shall presume that our problem is that 
of studying uncoupled vibrational relaxation behind a normal -shock wave (of 
course, ignoring effects of dissociation) in a dilute mixture of oxygen molecules 
in an inert diluent of argon atoms . The fluid temperature T is impulsively 
changed by the shock wave from ambient T 0 upstream of the shock to some 
large constant value T immediately downstream of the shock. For shock tem- 
peratures of T = 4000° K (not plotted), there is little perceptible difference 
between the two solutions. At the higher shock temperatures, however, the 
introduction of the truncation factor causes more rapid vibrational relaxation. 
This is illustrated by the fact that the solid-line curves are closer to the 
equilibrium asymptote (Ty/T = 1) than the dashed line curves regardless of 
the time t/r . This is perhaps better illustrated by the following table, where 
the relaxation time required for the solutions to approach to 10 percent of their 
equilibrium value is compared. 

Table D-l 

Comparison of characteristic 90-percent relaxation time 


T 

(t/r)i 

(t/r) 2 

8,000 

2.12 

2.28 

12,000 

1.94 

2.51 

16,000 

1.75 

2.75 


We see that the stronger shocks (indicated by the larger values for T) 
yield the greater differences. There is nearly a factor of 2 difference in the 
last case listed. These differences occur only because the ratio £y (Ty)/ Z y (T) 
is not unity (see discussion in Chapter II). Changing the functions qy in equa- 
tion (D3) affects primarily the driving force qy (Ty) - qy (T) and would be 
expected to have only a negligible effect on the relaxation time. Hence, it may 
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be anticipated that investigations, such as that conducted by Maillie and Hsu (57), 
will show only a minor effect of "anharmonicity." 

The present model contains a number of terms of the type given by equation 
(D2a) (e.g., see eq. H-95). The parameter a , which is equivalent to N in the 
case of the first term of equation (H-95a), has values less than 26, and the dif- 
ferences exemplified in table D2 are even greater. Figures D-2 illustrate the 
differences when a = 11 . (The quantity qy has the same functional form as 
given above except that the parameter N is denoted instead by a ; compare, 
e.g., eq. (A27) and (A42).) In this case we obtain the characteristic relaxa- 
tion times given by the following table, where the values have the same mean- 
ing as previously: 

Table D-2 

Comparison of characteristic 90-percent relaxation time. 


T (°K) 

(t/Oi 

(t/r ) : 

8,000 

1.50 

2.40 

16,000 

0.98 

3.60 


At the higher temperature the values differ by nearly a factor of four and the 
effect of the truncation factor is even greater. 
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Time, t It dimensionless 


(a) T = 8,000° K 

Figure D-l. Comparison of relaxation behavior for case of uncoupled vibrational relaxation where 
a representative number of vibrational energy states are included 








Time, t/r dimensionless 


(a) T = 8,000° K 

Figure D-2. Comparison of relaxation behavior for uncoupled vibrational relaxation where only 

a portion of the vibrational energy states are included 




APPENDIX E 



DEVELOPMENT OF THE NUMERICAL PROCEDURE FOR 
SOLVING THE MODEL AND FLUID-FLOW EQUATIONS^ 


(E-l). Introductory Comments 

The differential equations that describe the internal properties of the model 
molecule were derived in Chapter II; special solutions not requiring integration 
were obtained and discussed in Chapter in. To obtain the "complete" solutions 
described in Chapter IV necessitates a coupling of the model rate equations 
with the fluid -flow equations and requires a fairly involved numerical proce- 
dure for integrating these equations. In particular, the model equations will 
here be coupled with the one -dimensional fluid-dynamic equations to obtain a 
system of nine first-order nonlinear ordinary differential equations. A discus- 
sion will then be given of the numerical procedure used to integrate these 
equations . 

It is beyond the scope of this paper to analyze the model equations for all 
possible flow fields. One can infer the behavior along streamlines in a complex 
flow field, however, by reasoning from the results obtained for shock waves and 
nozzle flow. Furthermore, considerable data are available for such one- 
dimensional and quasi-one-dimensional flows, and results may be compared 
quantitatively both with experimental data and with other vibration-dissociation 
models. Although the time -independent one -dimensional fluid problem is one of 
the simplest to attack, such studies are not without difficulty (18, 29, 77). The 
major problem is that of integrating a system of equations containing greatly 


t The author wishes to draw attention to the many valuable discussions he has 
had with his colleagues Mr. Harvard Lomax and Dr. Harry E. Bailey, of NASA 
Ames Research Center, concerning the numerical integration of the model 
equations given in this paper. The numerical work was undertaken while they 
were completing the research for several papers that are now published (51, 
52, 53, 53a, 54). The discussion here may be considered a brief summary of 
the material relevant to this work and contained in those papers. 
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different time constants . This problem is often referred to as that of integrat- 
ing, "stiff” equations (29, 25); the reason for this designation will be apparent 
later. 


(E-2). Fluid-Flow Equations 

The model equations to be integrated are given by equations (11-110). These 
equations must be solved together with the fluid-dynamic equations to compute 
the properties of flow fields subject to vibration-dissociation coupling. We 
make the following basic assumptions: 

(1) The flow is steady, one-dimensional, inviscid, nonconductive, 
and nonradiating. 

(2) The fluid behaves as a mixture of perfect gases. 

The equations for conservation of mass, momentum, and energy, respectively, 
are then 


pvA(x) = M 

(Ela) 

vdv + ^ = 0 
P 

(Elb) 

V 2 

h + -= ht 

(Elc) 


where p and v are the mass density and speed of the fluid, and A(x) is the 
cross-sectional area of a streamtube or a channel, a known function of the 
distance x measured along the flow, h is the specific enthalpy, and M and 
h^ are constants. In derivative form with x as the independent variable the 
equations are 


pA 


— + vA ^ 
dx dx 



(E2a) 


dv dp 

pv T" + ^ = 0 

dx dx 


(E2b) 


dv dh 

v 4 

dx dx 


0 


(E2c) 


r 


These constitute three equations for the four unknowns v, p, p, and h. To 
complete the set we must introduce additional relations specifying the thermo- 
chemical state of the gas. These are given by the thermal equation of state 

p = p r 0 T Z y j 
j 

(E3a) 

= pRZT 

(E3b) 

= Z n i kT 

i 

(E3c) 

and the caloric equation of state 


h = Xrjhj 

(E4) 


j 


where H 0 is the universal gas constant (gas constant per mole), R = R 0 J] ( y j)o 

j 

is the low-temperature specific gas constant (gas constant per unit mass), T is 
the local or "translational” temperature of the gas, Z = ( y j)o * s 

j j 

compressibility factor, 7j is a measure of the concentration in units of moles 
per unit mass of gas, and the concentration variable subscripted o means 
"evaluated at low temperatures." The enthalpy per mole, hj , is defined 
according to 


hj = R 0 T + e j (T) + qj (Tj) + h Gj 

= hj (T) + ^ 
r i 


(E5a) 

(E5b) 
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where ej (T) is the internal energy per mole and h G j is the heat of formation 
per mole (see Appendix G for values). The asterisk (*) is applied only when 
the subscript j denotes a molecular species and only when the vibrational 



energy is -considered separately in the internal energy or enthalpy. In this case 
the variable qj , already defined, represents the vibrational energy per mole, 
and the variable ej the vibrational energy per unit mass. The functions ej (T) 
and their asterisk equivalents ej (T) are defined in Appendix A. 

In general, we can consider the nonequilibrium thermochemical state to 
be specified by the two state variables p and T and the nonequilibrium vari- 
ables and . In terms of these variables the equations of state have the 
functional form 


P = P(P.T, 7j) (E6a) 

h = h(T,r j> e j ) (E6b) 

With the state equations we can eliminate p and h from equations (E2). By 
introducing the rate equations (11-114) and (11-110), we then obtain the follow- 
ing complete set of differential equations: 


» dv .dp 
p A — + vA — 
dx dx 


= "PvA x 


dv 9pdp Bp dT Bp d7A Bp d7B Bp d7a 

^ V dx Bp dx BTdx By^ dx Byg dx 8y a dx 
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dx 


= Q 


y A 


dyg 

dx 


= Q 


y B 


dg A 

dx 


= Q 


e A 


d€ B _*B 

~ - W 

dx 



(E7) 


The Q's, with appropriate superscripts, denote the right-hand side of the rate 
equations (31-110), divided by the velocity v . These equations are written for 
one -dimensional flows in general. The results for flow behind a normal shock 
are obtained by setting the cross-sectional area constant, so that dA/dx = A x = 0. 

Equations (E7) can be written in matrix -vector form if the vector w* is 
introduced such that its transpose is w* T = (v,p, T, 7^, 7 a ), the 

vector c* such that c * T =(-p vA x ,0,0, A ? Q e B Q y a) 5 anc j the 

matrix [B*] having elements that are the coefficients of the derivatives 
dw*/dx. The result is 


[B*] ^ = c * <E8) 

In the general case the matrix [B*] will become singular when the velocity 
v has the value corresponding to the local speed of sound af (e.g., see Appen- 
dix F) . For this reason it is advantageous to introduce a new independent 
variable s defined such that 

dx 

— = x' = a' det ([B*]) (E9) 

there cr is an arbitrary scaling constant that is equal to 1.0 x 10 -2 ° for normal- 
shock flow and equal to -1.0 x 10“ 2 ° for nozzle flow. If we multiply both sides of 
equation (E8) by the product det ([B*]) [B*] 1 and use equation (E9), 

^ = w*' = F* = det ([B*]) [B*]" 1 c * (E10) 

QS 

where det ([B*]) [B*] 1 is the adjoint of [B*]. Finally, we define the new 
vectors w and F with one element more than their starred counterparts 
according to 


w T = (w*T, x ) 


(Ella) 



that is, 


W T = (v,p, T, y A , y B , e A , e B , y a , x) (Ellb) 

and, further, 

F T =(F* T , det(lBl)) (Ellc) 

The explicit representations for the elements of these vectors are given in 
Appendix F, along with a minor alteration that facilitates the integration of 
the equations. With the abbreviated notation we then have a set of nine simul- 
taneous equations represented in vector form by 

s w' = F(w) (E12) 

In this there is no explicit dependence of F on the independent variable s ; 
that is, the equations may be described as being autonomous. These equations 
and their dependent variables are analogous to the set of equations and variables 
investigated by Lomax and Bailey (53, 53a). The equations are first-order, non- 
linear, ordinary differential equations and, ingeneral, are not easily solved. In 
the next section the difficulties associated with these equations will be dis- 
cussed, and the numerical procedure required for their solution will be 
described. 

(E-3). Review of the Numerical Integration Methods 

The equations described in the previous section, equations (E12), have 
features not unlike those of similar equations encountered in the study of 
effects of chemical nonequilibrium in fluid flow. Considerable research 
and computer time has been spent on the numerical integration of such 
systems of equations, and many different methods have been developed for 
this purpose (63, 68, 69, 95). These methods fall into two principal categories. 



In one, the nonlinear differential equations are reduced to nonlinear difference 
equations by the substitution of difference -differential expressions for the de- 
rivatives, and the equations that result are then solved numerically (27, 28, 

29, 30, 31, 32, 95). In the other, the differential equations are first linearized 
locally, and the resulting linearized form is solved either exactly (63, 68, 69, 
72), approximately (26), or by finite -difference methods (96). On reviewing 
these methods one observes that three basic considerations occur regarding 
a comparison of the methods: accuracy, stability, and efficiency. The precise 
definitions of these terms are given in references 53, 53a, and 54. Rather than 
involve the present discussion with the intricacies of these definitions, it suf- 
fices to state that (1) a numerical method is accurate if it produces a solution 
that would agree with an analytical solution (if one could be obtained) within a 
controllable error (i.e., an error that may presumably be made small by reduc- 
ing the size of the difference interval); (2) a method is stable if its solution does 
not diverge; and (3) a method is more efficient than other numerical methods if 
it requires less computer time to obtain the same numerical solution, numerical 
error being roughly the same. 

From the experience gained with the various methods, instability appears 
to be the most troublesome problem in nonequilibrium studies in which conven- 
tional integration methods are used. An unstable method can often be made 
stable by reducing the integration step size (29). This, however, often decreases 
the efficiency of the method (many steps are required while the dependent vari- 
ables change by only very small amounts). One also finds that methods that are 
always stable may not always be the most efficient or even the most accurate. 
The interplay of stability, accuracy, and efficiency is thus extremely important 
in numerical integration methods. Such considerations form the subject of 
three recently published papers by Lomax and by Lomax and Bailey (53, 53a, 

54). Their recommended procedures will be followed here. Since an under- 
standing of a few of their findings is essential for the proper utilization of the 
numerical procedure to be described, it is worthwhile to digress and discuss 
their theories briefly. 
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Lomax and Bailey have concentrated their investigations primarily on 
coupled linear ordinary differential equations. They also discuss, however, 
how their results extend to certain nonlinear systems; specifically, those that 
may be characterized as quasi-linear (i.e., in which the highest order deriva- 
tive is linear). The equations of the previous section are of this type. To apply 
their theories it is advantageous first to obtain a linearized form of equa- 
tions (E12). Consider the set of autonomous equations 

= w' = f(w) (E13) 

where the elements, say F* , have no explicit dependence on the independent 
variable s . If the right-hand side of the equation, Fi , is expanded about a 
local point, referenced as n where s = nh, we obtain 

w i = ( F i)n + E ( w j -w jn) (ff:) n + 0[(w-w n ) 2 ] (E14) 

We define 





(E15) 


as the elements of the matrix [A]^ , where the subscript n signifies that the 
quantities are evaluated at the discrete reference point denoted by n. If 
w = Wjj+j we note that the higher-order terms in w - w n can be written 
h 2 [w n+ i - WnJYh 2 or h 2 [w^] 2 plus terms of order 0(h 3 ) . Thus in vector- 
matrix form we have 


W' = [A] n w + (F) n - [A] n w n + 0 (h 2 ) (E16) 

The above equation is, after higher-order terms are neglected, a locally lin- 
earized form of the original equations. Further, since the original equations 
were autonomous, the linearized equations have constant coefficients. 


Equations (E16) are readily solved (e.g., see 54) and the homogeneous solu- 
tion may be written ^ 


w i = E C t j exp (Xj s) (E17) 

j 

where the C^j are constants depending on the initial conditions and that Xj are 
the eigenvalues associated with the system of equations; that is, they are the 
solutions of the characteristic equations [A^j - Xj 6^] = 0 (<5-^ is the Kronecker 
delta function such that 6^ = 1 or 0 for j = i or j f i) . It is typical of equa- 
tions for nonequilibrium flow that in certain regions of a flow field the eigen- 
values Xj are large, negative, real numbers and some are much smaller in 
magnitude than others. (These characteristics are also found in the differential 
equations that describe the dynamics of "stiff” springs. For this reason 
systems of equations having such properties are often referred to as "stiff” 
equations (25)). 

To facilitate explanation of the problem, consider a system of two differen- 
tial equations with eigenvalues represented, for example, by ( x i) n = -1000 /i 
and (x 2 ) n = -M , respectively. In the integration of equations with such eigen- 
values, two cases can occur. (1) If the effect of (x t ) n over the small region 
where exp(-1000£is) is significant and is to be resolved, then we must perform 
our calculations at points spaced very close together. (2) If the value of s is 
large enough that exp (-1000 jus) is negligible compared with exp(-/i*s), then we 
need use only the much coarser spacing that resolves exp (-jits). However, if 
the integration is carried out numerically with many of the conventional integra- 
tion schemes (e.g., "explicit" methods such as the fourth -order Runge-Kutta 
method) , violent instabilities occur for the coarse spacing. Lomax and Bailey 

tA different relation is required when the eigenvalues Xj are not all distinct. 
This point is considered in reference 54 but will be ignored here as not par- 
ticularly essential either to the discussion or to the analysis. The eigenvalues 
are not actually required in the numerical methods, but are only monitored 
occasionally as a check on the integration procedure. This point will be made 
clear in the discussion. 



(51, 53, 53a) refer to those like (x 1 ) n as "parasitic” eigenvalues and those like 
(\ 2 ) n as "driving" eigenvalues. 

Just how the eigenvalues affect the various numerical methods canbe under- 
stood from the paper by Lomax (54). He shows in a rigorous manner that if 
each equation of the set of differential equations is operated on by identical 
linear difference -differential operators (we define such operations as "the 
numerical method") then the numerical method "detects" the eigenvalues of 
the differential equations and the success or failure of the method is measured 
by 

(1) Its accuracy in resolving the eigenvalue for which it is most accurate. 

(2) Its stability with respect to the eigenvalue for which it is most unstable. 

Briefly, we may say that linear numerical methods, in effect, "decouple" the 
system of equations, and one needs only to determine whether a method has the 
ability to resolve the effects produced by the separate eigenvalues. The accu- 
racy and stability criteria of a number of numerical methods are cataloged by 
Lomax and by Lomax and Bailey in references 53 and 54. The efficiency of a 
few of the more useful methods in the integration of stiff equations, is discussed 
in reference 54. Lomax and Bailey conclude that at least two types of difference- 
differential operators are required for rapid integration of the equations of non- 
equilibrium flow. Two satisfactory methods are the explicit and implicit 
methods described in the following discussion. 

The explicit method used by Lomax and Bailey (53) with considerable suc- 
cess,^ and employed to obtain the solutions of Chapter IV, is referred to as 
either an "Euler predictor with a modified Euler corrector" or as a 


^ Lomax and Bailey (51, 52, 53, 53a, 54) do not extol any particular method or 
combination of methods as being universal techniques. On the contrary, they 
have investigated and reported on the merits of a number of methods. 


"second-order Bunge-Kutta method". The method is defined by the difference- 
differential equations 


4+i = w n + h F(wn) 
h — ► (l) 

w n +i = w n + -[F(w^!) + F(w n )] 


(E18) 


The quantity h is the interval size and the superscript in parantheses denotes 
the iteration number. (The superscripts are omitted from the symbols that 
represent the final values of the dependent difference variables. Note the omis- 
sion of (2) on w n +i in the second equation above.) Equations (E18) have a 
truncation error (xh) 3 / 6 when applied to linear equations and a real stability 
boundary of -2.0 (see (54) for the definitions associated with the terminology 
and for a discussion on these given values). If xh > -2.0 then numerical insta- 
bility will not occur. We may exemplify this further by returning to the "two- 
eigenvalue" problem discussed earlier. We see that stability will be assured 
provided the integration interval size is the smaller of the two values 
(h}) n < 2 x 10““ 3 /p and (h 2 ) n < 2/p. The smaller value, of course, is (h 1 ) n and 
is the interval based on the "parasitic" eigenvalue. Also the effect of trunca- 
tion error on the dependent variables is, to a great extent, determined by the 
"driving" eigenvalue. Although the solution values are little affected by the 
parasitic eigenvalue, stability depends on this eigenvalue, and if the effect of 
the parasitic eigenvalue could be ignored, a three-order-of-magnitude increase 
in interval size would be possible, in principle. Implicit methods have such 
desired characteristics. 

The implicit method used by Lomax and Bailey (53) is the "modified Euler 
method" defined by the difference -differential equation 

w n +! = w n + - [F (w n+1 ) + F (wn)] 


(E19) 




This equation contains derivatives that must be evaluated in terms of the un- 
known dependent variables w n +i . Application of this equation is facilitated by 

_ dw n +i 

the following procedure. If we substitute Wn+i = — given by equation (E16) 

for F (w n +!) , the equation can also be written 

W n+1 = W n + | {[A] n w n+1 + (F) n - [A n ] w n } = 0 (h 3 ) (E20) 

or we can combine terms and obtain the set of simultaneous equations given by 

[[I] -|[A] n ](w n+1 -w n ) = hF n + 0 (h 3 ) (E21) 

where [I] is the unit matrix (matrix with diagonal elements of unity and off 
diagonal elements of zero). These equations may be readily solved for the 
unknown quantities w n -H , with the result^ 

W n+1 = W n + h jjl] - | [Ain]' 1 (E22) 

Although, the equation may now be considered explicit since the unknown quan- 
tities are now related to quantities that are readily evaluated, we shall continue 
to refer to this method, and as well as methods requiring the evaluation of the 
Jacobian matrix A, as implicit. It can be shown (53,54) that the above method 
is unconditionally stable for all real negative values of kh . (The truncation 
error is (\h) 3 /l2 when applied to linear equations. See 54.) 



^In practice the matrix inversion that yields equation (E22) may cause machine 

"overflows" if the matrix j^[I] - — [A] n J is not first properly scaled. This is 

easily done. If we let [S] be a diagonal matrix of which the elements are the 
reciprocals of the maximum -valued elements of the row vectors of the matrix 

to be scaled, then Sjj jll] - “ [A] n j^ — 1 and inversion difficulties are not 


encountered. Thus, in practice, the inversion procedure given by 


w n+1 = w n + h [S] 


[m - 1 Mn]f 


[S] F n 


(E23) 


is to be preferred. Equations (E22) and (E23) are exactly equivalent. 
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With the use of the implicit method we have gained 1 'stability M when the 
parasitic eigenvalues are important, but at the expense of additional computa- 
tional complexity. Specifically, the elements of the Jacobian matrix, a^, must 
be computed. Although additional computer time is required to advance a step 
with the implicit method, considerable saving in the computational time required 
to obtain the entire solution may be obtained since large step sizes are possible. 
In particular, the intervals are large relative to those required for explicit 
methods that are dominated by stability criteria. 


The application of the implicit method is made considerably easier if one 
obtains the matrix elements by numerical rather than analytical differentiation. 
The procedure to use is that suggested by Lomax and Bailey (53) and is indi- 
cated by the formula^ 


a 


ij 


3F i AF i F i (l.01xw i )-F i (0.99xw j ) 
3wj Awj 0.02 x wj 


(E24) 


T Some of the individual elements, aij , obtained by numerical differentiation, may 
not agree with the similar elements obtained by evaluation of the relations found 
by analytical differentiation. This occurs when there are large -magnitude dif- 
ferences between the individual terms. For example, if 


where 





a i,k w k| 10 9 a i s k-i w k-i 


(E25) 


one obtains "numerically" a^ = 0 and "analytically a^= a^ ^ using an elec- 
tronic computer with an eight digit mantissa. The full implication of this dif- 
ference is not yet understood, and the disagreement is pointed out only as a 
potential source of difficulty. Actually, neither the author nor Lomax and 
Bailey (53) encountered any difficulties from equation (E24) (except when try- 
ing to compare matrix elements). Considerable time was expended and diffi- 
culty experienced, however, in trying to evaluate and use analytically obtained 
matrix elements. The Jacobian matrix contains 81 nontrivial elements for the 
case studied in Chapter IV, and only a few are identical zero. 
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One additional numerical parameter requires some discussion. This is the 
interval size "h” contained in both integrating methods. During the course of 
a numerical integration the interval size was allowed to vary in a manner simi- 
lar to that discussed by Lomax and Bailey, but according to the following 
specific criteria. 

(1) If for all variables Wj we have [(wj) n “ ( w j) n « 1 ]/( w j) n _ 1 — 0.025 e , 
then the interval for the next point was doubled, that is, h n +! - 2h n , 

(2) If for any variable Wj, say w k , we have [(w k ) n - (w k ) n _J/(w k ) n ^ 1 > 0.1 e, 
then the interval size h n was halved and the (w) n were recomputed. 

(3) If neither condition (1) nor condition (2) was satisfied, then the vari- 
ables (w) n+1 for the next point were computed with unchanged 
interval, that is, h n+1 = h n . 

For most calculations e = 0.5 was used. With this scheme it turns out that 
any reasonable "initial value" h t may be used for starting the integration 
procedure. The above criteria were more stringent than the conditions used 
by Lomax and Bailey (53) where e = 1.0 . These criteria were required to 
assure sufficient accuracy for the cases investigated. 

The importance of using a combination of explicit and implicit methods can 
be illustrated by plots of the eigenvalues and of certain key parameters based 
on these eigenvalues for a few of the cases discussed in Chapter IV. In particu- 
lar, we shall exhibit the effects that occur when the stability criteria are not 
satisfied. Figures E-l show the eigenvalues and the parameters |^ m i n h| and 
| x max k | corresponding to the solutions given by figure IV-1 of Chapter IV. 

The similar figures E-2 correspond to the solution of figure IV-6. The abscissa 
scale here is the same as that in the figures in Chapter IV. The eigenvalues 
were computed for every third point of the solutions. In figure E-l the transient, 
quasi-steady, and final relaxation regions are again labeled by the numerials 1, 
2, and 3, respectively. In obtaining the solutions the first 125 points were com- 
puted with the explicit method, equation (E18), and the remaining points (roughly, 
between 200 and 250 points) were computed with the implicit method, equa- 
tion (E22). The "switch point," at which the integration method was changed, 


is indicated on these figures by a dashed vertical line. One observes that the 
equations are quite "stiff" throughout the entire flow region, since all the large- 
magnitude eigenvalues are negative and are very large compared with the driv- 
ing eigenvalues. The driving eigenvalues are the very small quantities having 
values of order 10“ 5 to 10~ 6 in figure E-l and 10“ G to 10“ 7 in figure E-2. The 
precise values of these small eigenvalues are not to be trusted, since there is 
some doubt as to whether the computer program can accurately compute eigen- 
values having values seven orders of magnitude less than the largest eigenvalue. 
One may be sure, however, that the small eigenvalues are no greater than those 
shown. (Note also that the small eigenvalues are distributed less randomly as 
the largest eigenvalue becomes smaller.) 

The fact that the equations are "stiff" appears to be a characteristic of 
vibration-dissociation coupling and may explain the computational difficulty 
experienced in references 18 and 77. Such behavior is usually not observed 
for shock -wave flow governed only by chemical relaxation (note the solutions 
exhibited in (53)), and explicit methods are then quite satisfactory. Actually 
the implicit method could well have been used throughout the entire flow region 
for the shock-wave cases that were investigated for this thesis; however, only a 
very small amount of computer time was expended in using the explicit method. 
Although, it may appear on the bases of the figures, that a large fraction of the 
entire flow region was computed with the explicit method, the reader is to be 
reminded that the variables are plotted with a logarithmic abscissa scale. 

The figures show that the eigenvalues vary primarily in the transient 
region and remain relatively constant in the other regions. Another feature 
to note is that figures E-l and E-2 exemplify strong and weak shock waves, 
respectively, and there is roughly an order-of-magnitude difference in the 
eigenvalues. The strong shock wave has parasitic eigenvalues with the larger 
values. The general character of the locus of eigenvalue points remains rela- 
tively unchanged independently of the strength of the shock wave. 

The effect of exceeding the stability criteria may be observed in figures E-lb 
and E-2b, particularly the latter figure. A comparison of these figures with 



figures IV-1 and IV-6 shows that whenever the quantity h| exceeds the 

stability criterion (i.e., whenever jx min h| > 2), then unstable behavior is 
exhibited in the solutions. This is particularly evident in a comparison of 
figures E-2b and IV-6, where e = 0.75 rather than e = 0.50 was used to 
find the interval size. The most dramatic effect, which exemplifies the value 
of implicit methods for overcoming numerical instability, is the almost imme- 
diate damping of unstable behavior directly after switching methods. Since 
the stability criterion no longer applies, there is then a rapid increase in the 
size of the integration interval. This is shown in figures E-lb and E-2b by the 
rapid increase in the value of |^ m i n h| . 

Figures E-lc and E-2c are plots of the quantity |x max h| , which is a mea- 
sure of the the influence of truncation error on the solutions (for a detailed 
description of this quantity, see (53)). So long as this parameter remains less 
than 0.01 we may be reasonably confident that the truncation error has little 
effect on the final values of the solutions. This condition is slightly violated, 
but not to the extent one would expect an error greater than 0.01 percent. 

Figure E-3b shows the parameters for nozzle flow. Here the sign of the 
scaling constant cr (see eq. E9) is changed, and the sign of the eigenvalues is 
changed accordingly. The parasitic eigenvalues are now positive, and the role 
of the parameters |^ m i n h| and |x max hj is interchanged. The implicit method 
was used exclusively for obtaining the nozzle-flow solutions; as a result, the 
stability criterion, j^max hj , is always greater than 2. The f, error n crite- 
rion, |\£nin h| > reaches a maximum value of 0.12 about 1 -centimeter down- 
stream of the throat, indicating that the maximum error due to numerical 
truncation effects is slightly greater than 0.01 percent. Everywhere else in 
the flow the error is much less. It should be mentioned again that the values 
of the stability and error parameters depend on the value of e in the procedure 
described for determining the step size. 



Eigenvalue 




(a) Eigenvalues, \ 

Figure E-l. Parameters characterizing the numerical aspects of shock wave 
solutions, versus adjusted shock passage time (associated model 
parameters are listed along first row of table IV-1) 
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(b) |^ min h| 
Figure E-l Continued 
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( c ) | ‘'•max h| 
Figure E-l Continued 
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(b) j^ m in hj 
Figure E-2 Continued 
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Adjusted shock passage time, sec 


< c > | x max h | 
Figure E-2 Continued 
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X, distance measured from throat, cm 


(b) j ''•max h| 
Figure E-3 Continued 
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( c ) i^min h ) 
Figure E-3 Continued 
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APPENDIX F 


EXPLICIT FORMULATION OF THE FLOW EQUATIONS 


The procedure for obtaining the differential equations of one-dimensional 
flow was described in Appendix E , and the required equations were given im- 
plicitly. Here the explicit formulation of these equations is set down. The 
singular behavior of the coefficient matrix [B*] (see Appendix E) is also 
discussed. 

The set of equations required to describe one-dimensional flow with 
vibration-dissociation coupling are given by equation (E8), where the coefficient 
matrix [B*] is defined by 


[B*] = 
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(FI) 


In general, chemically reacting systems with vibrational equilibrium also ex- 
hibit similar coefficient matrices in that most of the diagonal elements are 
unity (see, e.g., 53, 53a). 


The matrix of minors and the determinant of such matrix systems are 
readily evaluated. (This is the same as saying that such systems of equations 



are easily solved.) That this is true can be demonstrated as follows: We may 
partition the matrix [B*] into square submatrices according to 


[B*] = 


[A] 

[ 0 ] 


[B] 

[I] 


(F2) 


The determinant of such a system can be written 


det [B*] = det [ A] x det [I] - det [B] x det [0] (F3) 


However, the determinant of the submatrix fo], the null matrix, is zero and the 
determinant of the matrix [I], the identity or unit matrix, is unity. Hence, the 
problem of obtaining the determinant of the high-order coefficient matrix is re- 
duced to a simpler problem of obtaining the determinant of the lower order sub- 
matrix [A]. The result is given by 


det [B*] = det [A] = det 


pA pA 0 


pp 


9p 8p 
dp 8T 

0 ^ 

8T 


= A 


9h 


(pi P^h + u2 8E. _ pu2 
dp d T dT p dt 


) 


(F4) 


The singular behavior of [B*] is particulaidy important, since it is the 
source of integration difficulties. It will only be demonstrated here that such 
behavior exists when the flow velocity v is equal to the "frozen" value of the 
speed of sound. We first introduce the following relations: the specific heat 
at constant pressure 


c 


P 


-alL- 

0T p 


the isentropic exponent 


7 = 



(F5a) 


(F5b) 
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the frozen speed of sound 


'i 


a f =VrRZT 

(F5c) 

and the identity 


Cp — Cy — RZ 

(F5d) 

We can evaluate the determinant of the coefficient matrix (F4), 
these quantities. The result is 

in terms of 

de,[B*] - 4KI-U’) . | H] 

(F6a) 

= pA[ Cp (RZT -u 2 ) + RZu 2 ] 

(F6b) 

” pA - RZ > Op - RZ - 

! (F6c) 

= pA c y [YRZT - u 2 ] 

(F6d) 

= pA Cy [a | - u 2 ] 

(F6e) 


We observe that the determinant [B*] is zero when the flow velocity is equal to 
the frozen speed of sound a^ (see, e.g., 97). This singular point may be identi- 
fied as a saddle-point singularity and is the source of considerable computation 
difficulty near the throat of a nozzle (see, e.g., (53a)). 

The adjoint matrix is required before an explicit form of the flow equations 
can be given (e.g., note eq. (E10)). An expression for this matrix is readily 
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obtained once the matrix of minors is given. With the "partitioning scheme" 
described earlier, the matrix of minors is readily evaluated. We obtain 


Minor [B*] = 
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where 
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(F8a) 
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The adjoint matrix [B*] 1 det [B*] is obtained from the matrix of minors 
in the following straightforward manner: First, multiply each element of the 
matrix of minors by the factor (-l) 1+ j, which changes the sign of alternate 
elements. (The resulting matrix is called the cofactor of [B*] .) Second, inter- 
change rows and columns; that is, obtain the transpose of the cofactor. The 
resulting matrix is the adjoint matrix of [B*] and is given by 


Adj [B*] = [B*] 1 det [ B * ] = 


ap ah 
a p aT 

-uA^ 

aT 

at 



- uA ^EJlL 


-"‘*5 

' u p (fr “ It) 

P a|* 
H aT 

-pA|^ 

l>A 

pA d{T,y B ) 

oA 3 ,(P- h ) . 
pA a(T,e A ) 


pA 9(T^) 

ap 
" u ap 

u 2 A 

M 

-Ap 

r A 

~ A \ 

-A/3 € 

A 


~APy 

r a 

0 

0 

0 

|B*| 

0 

0 

0 

0 

0 

0 

0 

0 

IBM 

0 

0 

0 

0 

0 

0 

0 

0 

IBM 

0 

0 

0 

0 

0 

0 

0 

0 

IB* 1 

0 

0 

0 

0 

0 

0 

0 

0 

IBM 


(F9) 


The derivative vector F , equation (E10), is obtained from the vector- 
matrix product adj [B*] cr The elements of the resulting vector are the 
derivatives that are required. Explicitly they are as follows: 


F? = “puA x — - uA Y ~(Pi h ) 

1 p x 0T U a(T,i)^ 


F 2 = u 2 p 2 A : 




(FlOa) 


(FlOb) 


F 3 * = u 2 pA x - aYj 


(FlOc) 


= IB* |Q A 


(FlOd) 



(FlOe) 



F 5 = IB* I Q YB 
F e * = IB* | Q 6a 
F 7 = |B* | 

Fg = IB* I Q 7a 


(FlOf) 

(FlOg) 

(FlOh) 


where £ = 1, 2, . . 5 denotes the nonequilibrium variables 7^, y B> e A ’ 

6r>, and 7 , respectively. This formulation is useful in that its symmetry sug- 
- D a 

gests the generalization required for more complex systems, for example, 
chemical systems with additional reacting species. In such cases there will be 
additional rate equations (Q quantities); one additional equation is required for 
each additional dependent species-concentration variable. 

The complete system of equations, denoted by the vector F (see eq. 

— ► * 

(Ellc) in the text), is obtained by adding one more element to the set F 
to obtain 


F T = (f * T , det [b*]) (Fll) 

An alternative representation of the first three rate equations, equations 
(F10 a, b, and c) is useful. We first introduce the identities 

/ y y y\ 

Hi = R 0 T V3 A + Q B + Q a / (F12a) 

H 2 = h* y / A + h* rB Q 7B + Q 6a + Q B + hy a Q r a (F12b) 
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In terms of these relations the summation quantities in the first three equations 
may be- written 


5^«*‘ > ( rzh2 _ ™ Hi ) <Fi3a> 

Yj Q £ = p[(RZT-u 2)H 2 - uZHj] (F13b) 

£ 

Substitution of these quantities into the relevant equations (F10) yields 

F t = Ff=-puA[||r(RZT^|^- - H t ) + RZH 2 ] (F14a) 

F 2 = F 2 = P 2 a[u 2 (|^ - + RZH 2 - f^rH t ] (Fl4b) 

F 3 = F 3 * = Pa[u 2 RZT - (RZT-u 2 ) - (F14c) 


These equations rather than equations (F10 a, b, and c) were used to obtain the 
results reported in Chapter IV. 

One additional change was made. In practice it was found that equations 
(F14), the associated rate equations, and the choice of dependent variables indi- 
cated by equation (El lb) of Appendix E are not always easily solved on an 
electronic computer that is limited to arithmetic operations involving numbers 
between certain fixed bounds. These bounds are of order 10" 38 and 10 38 for an 
IBM 7094 electronic computer when the single-precision arithmetic option is 
employed. The reason for the difficulties that occur can be illustrated as 
follows: In the flow field immediately behind a shock wave the dependent vari- 
able Y b is proportional to exp(-E^/kT B ) (T^ = 300° K immediately behind 
the shock) , and for values of E^/k greater than 26,000° K, 7-g is less than 
10“ 38 and its value is rounded to zero by the computer. Since this variable 
occurs as a product in the rate equation Q a , equation (FlOg) (see also eq. 



(Il-llOd) in the text), considerable difficulty is then encountered until the tem- 
perature Tg is sufficiently large that exp(-E^/kTg) > 10~ 38 . Similar 
problems are experienced with the dependent variable 7 a - In addition, these 
variables also change rapidly in the transient region behind the shock, increas- 
ing in an exponential manner more than 33 orders of magnitude (see, e.g., fig. 
IV-1). The integration methods actually involve low-degree polominal "fitting." 
Consequently, many points and a corresponding large amount of computer time 
were required to advance through the transient region. Another problem was 
that of determining the temperature T in the transient region. If 7 is 
identically zero then so is €g and T B is then indeterminant. All these 
problems were alleviated by introducing the dependent variables 7^, 7^, T^, 
T b , y a in place of ? A , ? B> « B , y a where 

\ = inr A (F15a) 

r B = ^nr B (FI 5b) 

r a = in r a (Fl5c) 


The appropriate rate equations for these more convenient variables can be 
obtained by inspection. They are as follows: 


Fj = - P^[lr ( RZT " H i) + RZH 2 
F 2 = p 2 a|u 2 (gf - Rz) + RZH 2 - fp Hi 

F 3 = pA^u 2 RZT“~^ - (RZT - u 2 ) H 2 - u 2 Hj 

f 4 = |b*|q 7a 
f 5 = |b*|q Tb 


(F16a) 

(F16b) 

(F16c) 

(F16d) 

(Fl6e) 


(F16f) 




Fg = | B*| Q A 

f 7 = |b*|q Tb 

F. = |B*|Q Y a 


(F 16g) 
(F16h) 
(F16i) 

The quantities and H 2 are evaluated according to equations (F12), where 


F 9 = | B*| 


A = 


'q^A c 


q B = r 


c Va (T a >Q T a ] 

(F17a) 

o Vb (T b )Q T b] 

(F17b) 


B 


The rate quantities , Q^B , Q^A ) Q^B , and Q^ a are obtained from 
equations (11-113) (note eq. (11-114)). The rate equations represented by F 
were integrated to yield the results discussed in Chapter IV. 
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APPENDIX G 


PHYSICAL CONSTANTS AND OTHER REQUIRED PARAMETERS 


The fundamental physical constants required for the analysis are listed 
here for the reader's convenience. Also given are the assigned values of other 
constants used throughout the work. The reader is referred to the text for an 
explanation of the listed values. We have used c.g.s units throughout. 

(1) Physical Constants 


Avagadro's number 

N 0 = 6.02486X10 23 

(g-mole)" 1 

Planck's constant 

h = 6.62517xl0 -27 

erg-sec 


•R = 1.05443xl0“ 27 

erg-sec 

Boltzmann constant 

k = 1.38044xl0" 16 

erg °K~ 1 

Universal gas constant 

R 0 = 8.31662X10 7 

erg (g-mole) 


(2) Other required parameters 
(a) Molecular oxygen 


Anharmonicity separation 

a = 9 

(see Appendix A and B) 

Temperature separation 

b = 16 

(see Chapters III and IV) 

Number of energy levels 

N = 32 

(see Appendix A) 

Lower states energy 
level separation 

©* = 2234. 3° K 

A i 

(see Appendix A) 

Upper states energy 
level separation 

©^ = © B = 1800.2° K 

(see Appendix A) 


© a = 19,675°K 

(see Appendix A) 


©jj = 32,276°K 

(see Appendix A) 

Model dissociation 
temperature 

©N_i = 59,280°K 

(see Appendix A) 



Characteristic 

dissociation temperature 0^ = 59, 368° K 


Characteristic rotational 
temperature 
where cr (symmetry 
factor) is 2 


= 2.07cr°K 

K,o 2 

(Eo) M = 2256°K 


(X 0 )M = 00962 


Oxygen molecular mass m^ = 32 g/mole 
Heat of formation h^ = 0 

Electronic excitation parameters 


J g EJU,o 2 e Eza,o 2 { ° K) 
13 0 


(fo | o 2 —o 2 ~ 5.2587 
<Po |o 2 -o = 4.5939 

<Po |o 2 -Ar = 5.4454 
(b) Atomic oxygen 

Mass m Q =16 g/mole 

Heat of formation ho/R 0 = 29690°K 

Electronic excitation parameters 


i 

g Ei,i,o 


1 

5 

0 

2 

3 

228.0 

3 

1 

325.8 

4 

5 

22825 

5 

1 

48609 

6 

5 

106112 

7 

3 

110166 

8 

19 

126051 

9 

8 

137878 

10 

84 

142875 


1 

(see Appendix A) 

(see eq. (A2d)) 

(see eq. (AlOa)) 

(see eq. (AlOa)) 

(see eq. (E5a)) 

(Eq. (Aid)) 

(Eq. (B8)) 


(Eq. (E5a)) 
(Eq. (Aid)) 
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(c) Argon 


Mass m Av = 3 9 - 944 

Heat of formation h^ r = 0 
Electronic excitation parameters 


g/mole 


£ 

1 

2 

3 


s E£,£,Ar 

1 

5 

3 


e EiU,Ar ( K) 
0 

93144 

93751 


(Eq. (E5a)) 
(Eq. (Aid)) 


(3) The following table was often used for converting various energy units in 
references and is reproduced here because of its usefulness. 


Table G-l 


Energy Unit Conversion Factors 


Unit 

cm -1 

erg/molecule 

cal/mole 

ev 

°K 

erg/ mole 

1 cm' 1 

1 

(hc)=1.98618xlO' Ifi 

2.863 

1.23976X10'* 

(hc/k)=1.43880 

(hcN)=l. 19665x10® 

1 erg/molecule 

(he) -1 =5. 03479x1 Q 1S 

1 

1.441X10 16 

6.24196X10 11 

(k 1 ) =7. 244 07x10* 5 

(N)=6.02486xl0‘ 23 

1 cal/mole 

.3493 

6.938X10' 17 

1 

4.331X10 -5 

.5026 

4.18X10 7 

1 ev 

8.06608X10 3 

1.60206X10" 12 

2.309x10* 

1 

1.16054x10* 

j 9.65219X10 11 

1 °K 

(k/hc) = .69502 

1 

(k)=l .38044x1 0 -16 

1.990 

8.6164X10' 5 

1 

(kN=R o )=8.31696xl0 7 

1 erg/mole 

(hcN) _1 =8.35666xl0" 9 

(N' 1 )=1.65979X10' 24 

,2.392xl0' 8 

1.03603X10' 12 

(kN) -1 =l. 20236x10'® 

1 


N Avagadro's number 
c Speed of light 
k Boltzmann constant 
h Planck 1 s constant 
R 0 Universal gas constant 
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